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ABSTRACT 

We present the most extensive combined photometric and spectroscopic study to date of the enor- 
mous globular cluster (GC) system around M87, the central giant elliptical galaxy in the nearby Virgo 
cluster. Using observations from DEIMOS and LRIS at Keck, and Hectospec on the MMT, we derive 
new, precise radial velocities for 451 GCs around M87, with projected radii from ~ 5 to 185 kpc. We 
combine these measurements with literature data for a total sample of 737 objects, which we use for a 
re-examination of the kinematics of the GC system of M87. The velocities are analyzed in the context 
of archival wide-held photometry and a novel Hubble Space Telescope catalog of half-light radii, which 
includes sizes for 344 spectroscopically confirmed clusters. We use this unique catalog to identify 18 
new candidate ultra-compact dwarfs, and to help clarify the relationship between these objects and 
true GCs. 

We find much lower values for the outer velocity dispersion and rotation of the GC system than in 
earlier papers, and also differ from previous work in seeing no evidence for a transition in the inner 
halo to a potential dominated by the Virgo cluster, nor for a truncation of the stellar halo. We find 
little kinematical evidence for an intergalactic GC population. Aided by the precision of the new 
velocity measurements, we see significant evidence for kinematical substructure over a wide range of 
radii, indicating that M87 is in active assembly. A simple, scale-free analysis finds less dark matter 
within ~ 85 kpc than in other recent work, reducing the tension between X-ray and optical results. 
In general, out to a projected radius of ~ 150 kpc, our data are consistent with the notion that M87 
is not dynamically coupled to the Virgo cluster; the core of Virgo may be in the earliest stages of 
assembly. 

Subject headings: globular clusters: general — galaxies: star clusters 



1. INTRODUCTION 

Giant elliptical galaxies at the centers of galaxy clus- 
ters, often called brightest cluster galaxies (BCGs) or 
cluster-central galaxies, are extreme systems that pro- 
vide stringent tests for theories of galaxy formation. Cur- 
rent models of BCG formation predict dualistic histories 
where most of the stellar mass is formed very early in the 
highest-cr peaks of the dark matter distribution, while 
much of the assembly of mass occurs later through the 
merging of smaller galax ies that are already old (e.g., 
iDe Lucia fc Blaizotl 120071 ). 

Observations support and inde ed motivate the ol d 
stellar ages in this scenario (e.g., lLoubser et aLll2009h . 
However, the assembly situation is less clear. There 
are observational reports of high-z BCGs having similar 
sizes and stel l ar masses to those in t he local Universe 
(|Collins et al.l l2009t iStott et all 120111 ). Such findings 
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could contradict the theoretical expectations for late- 
epoch growjth J _a]idiojagh the evidence is still in dispute 
(e.g., lAscaso et all 1201 If ). 

The late assembly of BCGs can be considered an 
extreme example of the general theoretical picture of 
two-phase assembly for massive galaxies, where a seed 
galaxy grows gradually through the accretion of an ex- 
tend ed outer envelope, primarily through minor merg- 
iZolotov et all 120091: lOser et all 
iDommguez-Tenreiro et al.H 2011: 



ers ( Abadi e t al. 
2010HWang eTail 



ng . 

Font et all 1201 1[) T Direct evidence for this process has 
been found via number density, chemical, and kinemat- 
ical transitions in the halos of various nearby galaxies 
dSearle fe Zinnl 119781 iCarollo et al.l 120071: iCoccato et all 
1 20091 12010b iProctor et al.l 120091 lAlves-Brito et all 120091: 
Forbes et all 120111: iHuxor et all 120111: lArnold et all 
120111) . but BCGs provide arguably the clearest opportu- 
nity to study such transitions — especially in cases with 
massive, extended cD envelopes. 

A few cD galaxies are known to show dramatically 
increasing stellar velocity dispersions with radius, im- 
plying cD envelopes that are associated with the sur- 
rounding cluster and m ay originate in disrupted or tidally 
stripped galaxies f e.g.. Ipresslerlll979l: iCarter et ailll999l : 
IVentimigha et al.ll2010HNewman et al.1120111 ). Such tran- 
sitions can be mapped out to larger radii using the kine- 
matic s of halo globular clusters (GCs; e.g.. lRichtler et all 
120111) . The remarkable abundance of GCs around BCGs, 
of up to 50,000 per system, not only makes kinematic 
studies feasible but has also long implied an unusual for- 
mational pathway for these galaxies' outer regions (e.g., 
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2011 
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iForte et al.lll982l) . compared to the halos of most other 
galaxies which have relatively fewer GCs. 

M87 (= NGC 4486) is the nearest example of a mas- 
sive elliptical galaxy at the center of a cluster (Virgo) . It 
hosts a population of ~ 10 4 GCs along with a mild cD 
envelope, and is a natural target for exploring the build- 
up of BCG halos in detail. Deep imaging of the core 
regions of Virgo has revealed a host of faint tidal stream- 
ers around M87 tha t presumably originate from ongo- 
ing a ccretion events ()Mihos et al.ll2005t Uanowiecki et al.l 
|2010() . However, despite this system's proximity, its kine- 
matics and dynamics are not well constrained. 

M87 was in fact the first galaxy outside the Local 
Group to be successfully studied wit h GC kinematics, 
followe d by further e x tensiv e studie s (IHuchra &: Brodlel 
iMould et all 



1987; 

19971 ICohenl 
Hanes et alj 120011: 



19871 Il990t iCohen k Ryzhoyi 



2000 : iKissler-Patig fc Gcbhardt 1998; 
Cote et al.l I2001D . However, in 
recent years the observational focus has turned to 
other systems, leaving the current data set for M87 
relatively obsolete. The typical velocity precision is 
much worse than that possible with modern instrumen- 
tation, and there are minimal data outside of ~ 35 kpc, 
which is just where th e cD envelope begins to emerge 
(jKormendv et al.ll2009l ). 

One of the most notable results from these GC studies 
was the discovery of rapid rotation in the outer regions of 
M87, suggesting either spin-up from a major merger or 
shear from large-scale motions of the Virgo core. Another 
finding was a rapid outward increase in velocity disper- 
sion that implied a transition to dynamics governed by 
the large-scale cluster potential. More recently, a hand- 
ful of planetary nebula (PN) spectra have been used to 
argue for just the opposite : a truncated stellar envelope 
and little sign of rotation (Doh ertv et al.ll2~009f) . 

The goal of the present study is to revisit the kine- 
matics of M87's GCs by exploiting the new generation of 
spectrographs on large telescopes. The combined depth, 
accuracy, and field of view of the observations presented 
here are the best yet obtained for any galaxy. We aimed 
to map the kinematics over a large baseline in radius, 
covering the possible transition region between galaxy, 
cD envelope, and surrounding cluster. As we will dis- 
cuss, many of the previous findings about the M87 halo 
are overturned in the light of this extensive data set. 

Given the richness of the data, there are numerous 
topics of interest that could be examined but are not 
covered in this paper. In particular, the themes of GC 
sizes, metallicities, and phase-space s ubstructure are con- 
sidered in depth in separate papers dBrodie et alJl2011t 
Strader et al., in preparation: IRomanowskv et al.ll2011[) . 

The distance we adopt for M87 is 16.5 Mpc, which 
is consi stent with analyses o f surface brig htness fluctu- 
ation s (jBlakeslee et al.l 120091 ) . GC sizes dMadrid et al 



2009ft, and the tip of the red giant branch (|Bird et a! 
20101 ). One arcmin of angular distance corresponds to 
4.80 kpc. 

The effective radius R e and luminosity of a galaxy 
like M87 are highly un certain and perhaps il l-posed 
quantities. For reference, iKormendv et al.l ([2009) found 
R e ~ 3.2' (16 kpc) and M v ~ -22.86. 

This paper is outlined as follows. $2] presents the new 
GC observations, while <|3] discusses the spectroscopic 



data reduction. <21 compares and combines old and new 
data, and provides careful classification of GCs and "con- 
taminants" such as foreground stars. The photometric 
and kinematical properties of the GC system are ana- 
lyzed in $5] and $6j respectively. <J7] provides a brief dy- 
namical analysis. $8] discusses several implications of our 
results, including the classification of GCs and the for- 
mation of M87 and BCGs in general. Sj9] summarizes the 
details of our findings. 

2. DATA 

Here we outline the acquisition of the new GC data 
sets presented in this paper. Spectroscopic data from 
Keck/DEIMOS, MMT/Hectospec, and Keck/LRIS are 
covered by $2.11 12. 2[ and 12.31 respectively. A log for all 
spectroscopic observing runs can be found in Table [TJ 
Those interested in an overview of the spatial positions 
of the GC spectra can glance forward to Figure [TTJ The 
archival CFHT imaging is discussed in §2.41 

This paper does not incorporate the small data set 
of velocities obtained with DEIMO S around the outer 
"strea m A" , which is presented in IRomanowskv et all 

(poll . 

Note that we will use the usual term "radial velocity" , 
referring to the observer's line-of-sight, which should not 
be confused with other quantities involving the internal 
(projected) radius of M87. 

2.1. Keck/DEIMOS Data 

We selected GC can didates from th e Subaru/Suprime- 
Cam imaging study of T amura et al.l (|2006al ). who iden- 
tified potential clusters using BVI imaging. Preference 
was given to objects with / < 23 and to blue candidates, 
since most of the clusters are large radii are expected to 
be metal-poor. 

Spatially, the iTamura et al.1 ()2006al ) data extend di- 
rectly east from M87 to an approximate projected galac- 
tocentric radius of ~ 2 deg, encompassing also the ellip- 
tical galaxy NGC 4552 at a radius of ~ 1.2 deg. We de- 
signed four Keck/DEIMOS slitmasks, each covering an 
area of ~ 16' x 5'. These were arranged with minimal 
overlap in a 2 x 2 grid covering a total area of ~ 32' x 10', 
starting ~ 8' from the center of M87. 

The two inner masks had 70-80 candidates each and 
the outer masks 50-60 GC candidates each, for a total 
of 254 possible clusters (two objects were duplicated on 
another mask, giving 256 slitlets). Not all of these were 
high-quality candidates; objects outside of the nominal 
range of color expected for true M87 GCs, and objects 
likely too faint to produce a usable spectrum, were used 
to fill unused mask area. All slits were cut with a width 
of 1" and a minimum length of 5". 

The four masks were observed using DEIMOS on 2007 
March 20-21, with an exposure time of 3-3.5 hours per 
mask, divided into individual exposures of 30 min. The 
weather was clear and seeing ranged from 0.8-0.9". We 
used the 1200 1/mm grating centered at 7500 A, yield- 
ing a resolution of ~ 1.5 A. This setup gives coverage of 
the full Ca II triplet (CaT, with rest frame wavelengths 
of 8498, 8542, and 8662 A) for most clusters with radial 
velocities typical of objects associated with M87. At the 
blue end, we also cover Ha for many clusters, which can 
serve as a valuable velocity check for low signal-to-noise 
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(S/N) spectra. This is our standard setup for GC obser- 
vation s with DEIMOS, us ed al so in | Romanowskv et all 
(|2009[) . lFoster et alJpOll . and l Arnold et all (|2011[ ). but 
here with longer exposure times to compensate for the 
fainter median magnitude of our targets due to the low 
densities of GCs in the far outer halo of M87. 

2.1.1. Additional DEIMOS Data 

As a check on some of the unusually high radial veloci- 
ties reported in the literature, we observed one additional 
mask in March 2010, centered ~ 7' south of M87 at an 
approximate position angle of 120 degrees. Candidates 
were split among those with known velocities and previ- 
ously unobserved objects. The total exposure time was 
36 min, split into two exposures, with approximately 1.0" 
seeing. The setup was identical to that for the earlier 
DEIMOS data. 

2.2. MMT/Hectospec Data 

The initial results of our DEIMOS observations mo- 
tivated us to obtain follow-up data with complete az- 
imuthal coverage. Hectospec is a fiber-fed instrument 
on M MT with 300 fibers c overing a 1 degree field of 
view (|Fabricant et al.l [2005) . GC candidates were se- 
lected from the CFHT imaging described in £12.41 Pref- 
erence was given to candidates with g < 22, but with no 
color restrictions beyond the basic GC candidate cuts. 
Objects were primarily selected at large radii, since lo- 
cal sky subtraction must be performed with individual 
sky fibers and sky subtraction in the inner regions of the 
galaxy (well interior to any of our fiber locations) is more 
challenging. 244 fibers were assigned to GC candidates, 
with the remainder to guide stars or sky. Each fiber sub- 
tends a diameter of 1.5". 

A single Hectospec field, centered on M87, was ob- 
served for 2 hours (divided into 20 min exposures) in 
February 2010. The conditions were clear and the see- 
ing was 0.6-0.7". Observations were made with the 270 
1/mm grating, which yields a spectral coverage of 3700- 
9100 A and a resolution of ~ 5 A. 

2.3. Keck/LRIS Data 

Most spectroscopic studies of GCs in M87 have tar- 
geted very luminous GCs, with typical i-band luminosi- 
ties of ~ 1.2 x 10 6 Lq, surpassing all but the brightest 
clusters in the Milky Way. With a goal of determining 
the metallicity distribution of somewhat fainter GCs in 
M87, in April 2010 we observed four slitmasks in the cen- 
tral region of the galaxy (~ 1.5-5.5') with Keck/LRIS. 
The stellar populations of these GCs will be discussed 
in a separate paper; here we consider only the radial ve- 
locities. As with the Hectospec observations, candidates 
were selected from CFHT imaging ( EI2.4jl . The median 
targeted object had g — 22.5, and no preference was 
given to objects of a certain color (beyond the basic GC 
selection criteria), so the objects observed should be an 
accurate reflection of the intrinsic color distribution of 
GCs at the galactocentric distances surveyed. 

Each mask had ~ 50 — 55 candidates, with 209 tar- 
gets total. All slits had a width of 1" and a minimum 
length of 4". At the time of the observations, the red 
side of LRIS was malfunctioning, so observations were 
made only on the blue side. All data were taken binned 



2x2, yielding a scale of 0.27" per pixel. We used the 600 
1/mm grating, which gives ~ 2600 A of spectral coverage. 
The exact wavelength coverage for each object depends 
on the spatial position in the mask, but a typical range 
was from ~ 3300 — 5600 A with the upper limit set by 
the dichroic cutoff. With a 1" slit, the resolution was 
~ 3.6-4.1 A depending on wavelength, with a binned 
dispersion of ~ 1.26 A per pixel. 

Total exposure times for each mask ranged from 2.5- 
3.5 hr, with individual exposures of 30 min. Seeing was 
variable between 0.6-0.9". 

2.4. Photometry 

The iTamura et all (|2006aT > imaging extends to ex- 
tremely large radii, but with non-uniform azimuthal cov- 
erage. For this reason, we chose to use for our default 
GC photometry the new gri CFHT/Megacam imaging 
of M87 published in lHarrisI ([200 9aH These data cover 
~ 1 square degree and so offer complete radial coverage 
to a projected radius of ~ 30' (144 kpc). 

We downloaded the stacked CFHT images from the 
CADC archive and reduced and analyze d the data in a 
manner very similar to that described in [Harrisl ([2009a); 
all of our measurements are consistent with that study, 
although we differ in the interpretation of a few points 
discussed below. The gri magnitudes were roughly cali- 
brated to the SDSS system by matching objects with the 
SDSS DR7 catalog and calculating a constant offset for 
each filter. The matches were done to the DR7 psfMag 
values (since the objects are unresolved) and over the 
magnitude range 18-22 for each filter. We cannot rule 
out the possibility of small variations in these zeropoints 
over the Megacam field of view, but the resulting accu- 
racy of photometric calibration is sufficient for nearly all 
of our scientific goals. The one exception is our investiga- 
tion of possible GC color gradients (see further discussion 
in g53). 

We corrected the photometry for foregro und red- 
dening based on iPeek fc Graved (|2010l ). which 

uses passively evolving g alaxies to correct the 

ISchlegel. Finkbeiner. fc Davisl (|1998f ) reddening 
map over the SDSS DR7 spectroscopic area. Note 
that we have applied spat i al va riations from the 
ISchlegel. Finkbeiner. fc Davisl |l998) map (which has 
a resolution of 6'), rather than taking a single mean 
reddening value. This exercise is important when 
studying the data on large angular scales, as there are 
highly non-uniform patches of Galactic cirrus in the 
Virgo direction, and we found variations in the E(g — i) 
correction of up to 0.02 mag. 

3. REDUCTION AND RADIAL VELOCITIES 

In this section we discuss the data reduction and 
determination of radial velocities for all spectra. 
Keck/DEIMOS, MMT/Hectospec, and Keck/LRIS data 
are covered in §3.1, 3.2, and 3.3 respectively. §3.4 com- 
pares velocities for common objects among the three in- 
struments as a consistency check. 

3.1. Original DEIMOS Data 

9 As Harris mentions, these data will soon be superseded by 
the Next Generation Virgo Cluster Survey (NGVS), which should 
allow superior selection of GC candidates. 
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The spectra were reduced in a standard manner by us- 
ing the IDL spec2d pipeline^- Two-dimensional images 
were flatfielded and wavelength calibrated using inter- 
nal lamps, followed by sky subtraction and an optimal 
extraction of one-dimensional spectra. 

We derived helioc e ntric radial velocities as described in 
iRomanowskv et ail (|2009t ). In brief, GC spectra around 
the CaT were cross-correlated with a library of DEIMOS 
stellar templates from mid-F to late-K. Confirmation of 
the velocity required visual identification of at least two 
lines; these were usually CaT lines, but in several cases 
a single Ca line and Ha were used. 

Radial velocity uncertainties were estimated from the 
width of the cross-correlation peak and by the variation 
in velocit y obtained from dif f erent templates. As dis- 
cussed in IRomanowskv et al.l ((2009) . cross-checks from 
multiple DEIMOS observations of the same GCs in pre- 
vious work and from different instruments all suggest 
that our DEIMOS uncertainty estimates are accurate in 
a random sense. However, these uncertainties do not re- 
flect the possibility of catastrophic misidcntifications of 
the principal cross-correlation peak. This is a particu- 
lar problem for low signal-to-noise spectra in the CaT 
region, which is prone to sky-line residuals. Hence we re- 
quire visual confirmation of multiple lines; nonetheless, 
such catastrophic failures have evidently occurred in the 
literature ( £j4.2[) . and we cannot exclude their existence 
from the present data. 

We classify our 254 candidate GCs as follows, based 
on the spectroscopic results. 63 appear to be bona fide 
M87 GCs, of which 60 are newly confirmed. A further 44 
are Galactic stars (see the i)4.4.2l for details of star-GC 
separation). 26 are background galaxies, mostly objects 
with multiple narrow emission lines. For the remaining 
122 candidates, no reliable velocity could be obtained, 
generally because of low S/N, A subset of these objects 
had cross-correlation velocities consistent with M87 GCs, 
but had at most one identifiable line; many could be true 
clusters. 

These numbers suggest a nominal success rate of 
~ 50% for identifying GCs, but the true percentage is 
higher, since many marginal objects were included as 
mask filler. This suggests that three-band imaging in 
decent conditions is an efficient method for finding GCs, 
even in the far outskirts of galaxies where the surface 
density of clusters is low. 

The basic data for all of these observed DEIMOS ob- 
jects are given in Table [3J The median velocity uncer- 
tainty for confirmed GCs is 9 km s _1 . 

3.1.1. Supplementary DEIMOS Data 

As discussed above, we observed one additional mask, 
closer to the center of M87, for 36 min in March 2010. 
Data reduction and analysis were identical to that for 
the earlier DEIMOS data. Owing to the short exposure 
time, reliable velocities could only be determined for 29 
objects of the 112 initial targets. 27 of these are GCs 
and two are stars. Data for these 29 objects are given in 
Table [3] 

3.2. MMT/Hectospec 
10 http:/ /astro. berkeley.edu/~cooper/deep/spcc2d/ 



The Hectospec data were pipeline reduced in a stan- 
dard manner as described in lMink et al.l (|2007[) . While 
this pipeline also produces radial velocities through cross- 
correlation with a template library, we chose to manually 
derive velocities for our objects. Heliocentric radial ve- 
locities for the Hectospec data were estimated through 
cross-correlation with a normalized Hectospec template 
of an M31 globular cluster over the wavelength range 
4000-5450 A. Re gions with significant sky emission were 
excluded. 

We estimated velocity uncertainties by Monte Carlo 
simulations using the error spectra. To these values, we 
add an error of 8 km s _1 in quadrature due to uncer- 
tainties in the wavelength calibration. The median total 
uncertainty for all of the spectra is 14 km s _1 . 

Of the 244 candidates, 240 produced an identifiable 
radial velocity. Of these, 172 appear to be true GCs 
and the other 68 foreground stars (see ^4.4|) . Since this 
sample of candidates is brighter than in the DEIMOS 
observations, a higher fraction of foreground stars and 
a lower fraction of background galaxies is expected, as 
observed. 

The basic data for all observed Hectospec objects are 
given in Table [4] 

3.3. Keck/LRIS 

The bulk of the LRIS data reduction was done with 
the XIDL/LowRedux packag^ZI- This pipeline performs 
the standard steps of bias subtraction, flat fielding, il- 
lumination correction, wavelength calibration, and sky 
subtraction before an optimal extraction. As the LRIS 
calibration lamps cannot be used for flat fields in the 
blue, the pipeline uses pixel flats taken with the grism in 
place during twilight. 

For approximately 15% of the slits, the pipeline ex- 
traction did not produce an optimal spectrum; in most 
cases this was due to very short slits for which the sky 
subtraction failed. These slits were reduced manually in 
IRAF using standard tasks. 

Unlike Keck/DEIMOS and MMT/Hectospec, LRIS 
has significant flexure — shifts of 1-2 pixels over two hours 
are quite common. Such shifts correspond to radial ve- 
locity shifts > 100 km s _1 and so must be addressed P^l 
Each individual exposure was corrected for flexure by 
cross-correlation of the sky spectrum with a template 
before coaddition of the exposures for each object. The 
proper location of the bright O I line at 5577.34 A indi- 
cates the overall fidelity of the process. 

Heliocentric radial velocities of the coadded spectra 
were derived through cross-correlation with a library of 
appropriate templates over the wavelength range 3900- 
5300 A. The intrinsic velocity uncertainties were esti- 
mated with Monte Carlo simulations using the error 
spectra. Total velocity uncertainties were computed by 
adding these values in quadrature with the wavelength 
calibration uncertainty and that due to the flexure cor- 
rection, which varied depending on the mask. The me- 
dian uncertainty for all of the spectra is 26 km s _1 . 
While we believe that these estimated uncertainties are 

11 http:/ /www. ucolick.org/~xavier/LowRedux/ 

12 LRIS flexure issues may be responsib le for the apparent veloc- 
ity zero-point offs et between the studies ofjCohcn & Ryzhov( 119971 ) 
and lCohenl ||2000|') . 
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reasonable, few of our objects have existing velocities in 
the literature because of their faintness (see §3J, and so 
only minimal external checks can be performed. Given 
the potential flexure issues, future study of a subset of 
these GCs would be valuable. 

The four masks contained 209 designed slits, and the 
data reduction yielded 204 usable spectra. 200 are prob- 
able M87 GCs, with only 2 stars, one background galaxy, 
and one ambiguous object (see ^4.4.3|) . for an excellent 
targeting efficiency of > 95% in the abundant inner re- 
gions of M87. 

The basic data for all observed LRIS objects are given 
in Table [3 

3.4. Self- Consistency of the New Radial Velocities 

The new radial velocities were derived from spectra 
taken with three different instruments on four separate 
observing runs, and so checking the consistency of these 
velocities is desirable. As these new measurements all 
have much smaller uncertainties than in previous studies, 
we first consider only new data; in §3]we compare our new 
velocities to those in the literature. 

Each of our data sets has differing coverage in galac- 
tocentric radius, azimuth, and target magnitude range. 
As such, the sample overlap is limited. Those objects ob- 
served with LRIS have no counterparts in our other data. 
There are five GCs in common between Hectospec and 
the original DEIMOS sample, and a further six between 
Hectospec and the supplementary DEIMOS data. 

The median difference between the combined set of 
DEIMOS radial velocities and those from Hectospec is 
15±5 km s^ 1 (all such errors quoted are standard errors). 
These are If the first and second DEIMOS runs are taken 
individually, the respective differences are 4± 10 km s _1 
and 18 ± 6 km s _1 . The distribution of the observed 
differences between the MMT and DEIMOS radial ve- 
locities is consistent with our estimated velocity uncer- 
tainties once these offsets are applied (see also Figure [H 
below). Given the small values of the offsets and the 
limited number of comparison GCs, we do not apply a 
common offset to the final set of velocities. However, 
we cannot exclude the possibility of systematic offsets of 
order ~ 10 km s _1 among these data. 

In Figure [T] we show sample spectra, over a restricted 
spectral range, from each of the instruments employed. 
We deliberately show spectra for each dataset that have 
io magnitudes near that of the median for that sample, 
so the plotted spectra can be considered representative 
of our data. 

4. A PHOTOMETRIC AND SPECTROSCOPIC 
COMPILATION 

This section describes our comprehensive catalog of 
photometry and radial velocities for candidate star clus- 
ters associated with M87. §4.1 compiles velocities from 
the literature. §4.2 compares our new velocities to those 
from the literature to establish a common velocity scale. 
In §4.3 we derive new half-light radii for many clusters 
using archival Hubble Space Telescope imaging. §4.4 dis- 
cusses the classification of objects, including those from 
the literature. §4.5 summarizes the catalog and §4.6 dis- 
cusses the separation of M87 GCs into metallicity sub- 
populations. 
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Fig. 1. — Sample spectra from each of the three instruments 
employed. Each GC was chosen to have a median io magnitude 
close to that of all confirmed GCs for that dataset. Representative 
spectral regions are displayed: 3900-4900 A for MMT and LRIS, 
and the Ca II triplet region for DEIMOS. Prominent spectral lines 
are marked. The spectra have each been shifted to zero redshift 
and smoothed with a 3-pixel boxcar for display. 

4.1. Literature Data 

To form a complete dataset for dynamical analysis, 
we combine our 451 GC velocities with data from the 
literature. Not all of these velocities are unique — 
there are 24 duplicates with the literature (with 5 ad- 
ditional stellar repeats), as well as some internal dupli- 
cates. As our starting point, we take the collecte d cat- 
alog of radial velocitie s from | Hanes et~aLl pOOll here- 
after H+01) (see also iCote et al.l 120011 ). Other than 
the new CFHT/MOS velocities published in that pa- 
per, the previous data in t h e compilation were take n 



fromlHuchra fc Brodiel (IT 9871) : IMould eTal] (Il987l fl99f]l) 



ICohen fc Rvzhovl (|1997D : and ICohenl (|2000D . 

To this compilation, we have added radial velocities 
for some objects observed as part of studies of "ultra- 
compact dwarfs" (UCDs) or "dwarf-globul ar transition 
object s" (DGTOs) i n the Virgo cluster. IJones et all 
(2006), IFirth et all (120081) . and IPaudel et all (120101 ) 
present low-resolution spectreP^I; we also include higher 
resolution studies of a s mall number of obje c ts, using 
Keck/ ESI, u ndertaken bvlHasegan et all (j2005[ ). IHaseganl 
(|2007l) . and lEvstigneeva et all (|2007l) . Together, these 
papers report radial velocities for 42 unique objects as- 
sociated with M87 or the Virgo cluster itself. In <JO]wc 
address the classification of these objects. 

H+01 used both linear fits and simple offsets to cor- 
rect all of the radial velocities to a common system. We 



13 There are also 6 velocities reported in Firth et al. (2009). 
However, these are all duplicate measurements from the literature 
and from our new data, and inter-comparisons of the velocities 
suggest that these authors' measurement uncertainties may be un- 
derestimated. Therefore we do not use them in this paper, pending 
further investigation. 
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adopt slightly different corrections than did H+01, pre- 
ferring constant offsets if possible. Th ere are 48 objects 
in com mon between CFHT/MOS and iCohen fc Rvzhovl 
(1997) with reliable radial velocities. The median differ- 
ence is 33 ± 19 km s _1 , which we adopt as a common 
offset. While there are only eight GCs in the overlap be - 
tween CFHT/MOS and the later study of lCohenl (pOOh . 
the evidence for an offset is even stronger (median differ- 
e nce 77 ± 33 km s -1 ), and we use this value. 

iMould et all (|1990t ) published original velocities as well 
as previously published one s from IMould et al.l ()1987[ ) 
and lHuchra fc Brodid p~987). The rep orted velocities for 
the fiv e clusters in common between iHuchra fc Brodiel 
(|1987|) and their work are straight averages. Considering 
the sixteen objects in common between the Mould et al. 
studies and the CFHT/MOS data, there is no statis- 
tically significant evidence for an offset in the velocity 
scale. There is marginal evidence for an offset between 
the velocities in IHuchra fc Brodiel (|1987l ) and those of 
both H+01 and Mould et al. However, the sign of the 
offset differs between the comparison samples, and so we 
choose to make no correction. 

Given the substantial set of new radial velocities and 
the large velocity dispersion of the GC system, these 
offsets actually have very little effect on the first-order 
dynamical analysis. However, they are relevant for the 
examination of outlying velocities and higher order mo- 
ments of the velocity distribution, as well as overall con- 
sistency. 

In addition to these zero-point co rrections, we fol- 
low iRomanowskv fc Kochanekl (|2001l ) in renormalizing 
the published uncertainties of the studies prior to 
H+01 by small factors to enforce consistency among 
the uncertainties of different papers (or, in some cases, 
providing uncertainties where none were published). 
In particular, we adopt a common vel ocity uncer- 
tainty of 10 6 km s" 1 for all objects from IMould et al] 
(fl987l [l990h : IHuchra fc Brodiel (Il987l) ; ICohen fc Rv^ovl 
(fl997T ) ; and a value of 77 km s" 1 for lCohenl (pOOOf ) . When 
combining velocities the weighted average was used. A 
few of the oldest (1990 or previous) velocities were ex- 
cluded from combined values if they deviated signifi- 
cantly from more than one recent measurement for the 
same object. 

There are two exceedingly popular objects in the vicin- 
ity of M87, S3143 and S1280, each of which has five 
published radial velocities. Regrettably, the latter is a 
foreground star. 

4.2. Merging the Datasets 

In E i3.4l we established — at least to the degree possible 
with the current observations — that our new radial ve- 
locities appear to be self-consistent and have reasonable 
uncertainty estimates. Now we discuss the integration of 
the current results with velocities from the literature. 

There are eighteen objects in common between the 
CFHT/MOS radial velocities reported in H+01 and the 

14 We refer to objects from the H+01 compilation with an 
"S" prefix s i nce t hey are on the revised numbering system of 
I Strom c t al. (1981). Rather t han ex t end tha t system further, all 
new GCs with imaging from lHarrisI H2009al) have IDs beginning 
with "H" ; those with imaging from Tamura ct al. ( 2006a) are given 
a prefix of "T", even though a subset of these objects may be in 
the Strom et al. catalog. 



present study. These matches are not all happenstance; 
the supplementary DEIMOS mask was designed to re- 
target a number of existing GCs with unusual reported 
velocities. Of the eighteen matches, fifteen show good 
agreement. The median difference between the new and 
old velocities is 7 ± 23 km s _1 and the uncertainty esti- 
mates appear to be accurate. 

However, three of the objects are "catastrophic" out- 
liers, such that the new velocity differs by > 5<r from the 
lite rature velocity (for other examples i n the literature, 
see lLee et al.ll2008l and lPark et aLll2010H . These objects 
are listed, together with the other matches, in Table HI 
The area around each of the GCs was checked for close 
companions to rule out target misidentification, but no 
viable interlopers were identified. The likely culprit is 
the incorrect choice of the primary peak in the velocity 
cross-correlation in the previous work. A general method 
to identify further outliers is unclear. While one of these 
objects (S923) is among the faintest in the H+01 MOS 
sample and is also peculiar as we shall see later, the other 
two (S878 and S1074) do not have unusual magnitudes 
or colors. It should be noted that we specifically targeted 
extreme velocities in our supplementary DEIMOS mask; 
our new observations do not provide strong evidence that 
the fraction of unreliable velocities in the literature is 
high (see also 34.4.3j) . 

T he other principal source s for the H+ 1 com pilation 
are ICohen fc Rv^movl ([1991 and ICohenl (pOOOl) . Since 
(at least nominally) these two papers have already been 
normalized to a common scale, we consider them to- 
gether. For thirteen objects in common between the 
present study and these two papers, the median differ- 
ence is 8 ± 24 km s _1 , and the distribution of the nor- 
malized differences supports the accuracy of the adopted 
uncertainty estimates. 

In Figure [2] we summarize these results. The top 
panel plots the radial velocity difference between our 
new data and the domin ant literature studies (H+01 and 
ICohen fc Rv"zliovl [19971) as a function of io- The scatter 
is consistent with the reported error bars, with the ex- 
ception of the three catastrophic outliers. The bottom 
panel, with identical axis limits, plots these same quanti- 
ties but for the repeat measurements between our MMT 
and DEIMOS spectra (§3.4). The higher precision of the 
new radial velocities is evident. 

4.3. Cluster Sizes 

For this study, there are two reasons why sizet(3 are 
useful. First, they can be used to distinguish unresolved 
foreground stars from actual GCs when the radial veloc- 
ity is ambiguous. In addition, a subset of our objects 
are probable UCDs with potentially different kinematics 
from the GCs; these objects can be partially identified 
through their large sizes. 

For as many GCs as possible, we collected sizes from 
the literature . The relevant papers included extens ive 
ACS studies ([Jordan et al.ll2009t iMadrid et alJl2009D as 
well a s the older WFPC2-based paper from lLarsen et al.l 
(|200lT ). In add ition, we obtained sizes for UCDs (or can- 
didates) from lEvstigneeva et all (|2008|) , Hascga n et al.l 

15 This term should be taken as a synonym for projected half- 
light radii for the remainder of the paper, unless otherwise 
stated. 
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Fig. 2. — The top panel shows the difference between our new 
radial velocities and the literature values as a f unction of ip mag- 
nitude , restricted to H+01 (filled circles) and Cohen & Ryzhov 
(1997;, open circles) for simplicity. The three "catastrophic" out- 
liers discussed in §4.2 are apparent. The lower panel shows the 
same quantities, but comparing our MMT and DEIMOS repeats, 
discussed in 53.4. 



(2005), and |Hasegan| (|2007l but see below). We took 
sizes for a total of 147 objects from these papers, but 
this is only a fraction of the total spectroscopic catalog. 
All sizes were scaled to a common distance of 16.5 Mpc. 

For the remainder of the objects, we searched the Hub- 
ble Space Telescope (HST) archive for images suitable for 
estimating GC sizes. We considered images taken with 
any of ACS, WFPC2, or STIS, and in any broadband 
optical filter. Owing to the large number of cycles in 
which WFPC2 was the default parallel instrument, the 
majority of useful images are from this camera. 

Our image analysis varied somewhat depending on the 
instrument used. For ACS data, we used the default 
drizzled images from the Hubble Legacy Archive, and 
estimated the point spread function (PSF) directly from 
a number of bright stars on the drizzled image. For 
WFPC2 data, we instead downloaded the basic cali- 
brated images from the MA ST archive. After co smic ray 
rejection with L.A. Cosmic ()van Dokkum 2001), the in- 
dividual chip images were averaged together. In order to 
preserve the PSF in severely undersampled WFPC2 im- 
ages, we only coadded images taken at the same position; 
we made no image shifts. TinyTinC3 was used to calcu- 
late position and filter-dependent PSFs for each WFPC2 
cluster. These theoretical PSFs were subsequently con- 
volved with the standard diffusion kernel that accounts 
for charge diffusion in the CCD. Charge diffusion is likely 
to be wavelength dependent and is one of the dominant 
sources of error in estimating sizes for barely resolved ob- 
jects in WFPC2 images. Finally, the few GCs with STIS 
data were analyzed in a similar manner to the WFPC2 

16 http:/ /www. stsci.edu/software/tinytim/tinytim. html 



images. 

Using these images and PSFs, we e stimated clus ter ef- 
fective (half-light) radii using ishape ([Larse 319991) . This 
program convolves models of star clusters with a supplied 
PSF and determines best-fit parameters by maximizing 
the likelihood function. Given the moderate signal-to- 
noise of most of our clusters on the images, we fit ellip- 
tical King models with a fixed concentration (r t /ro) of 
30 to each of the GCs. For nearly all objects, a fitting- 
radius of 5 pixels was used; we adopted a larger radius 
for a few unusually large clusters. The effective radii 
reported are circularized and assume a distance of 16.5 
Mpc. Besides the new size measurements, we have also 
remeasured the sizes of a handful of the fainter objects 
from the Hasegan catalogs, now using fixed concentra- 
tions (see lBrodie et al.ll2011l for further details). 

Determining uncertainties for cluster sizes is challeng- 
ing even with homogeneous data; the dominant uncer- 
tainties are systematic. With our range of instruments, 
filters, and exposure times, supplying an accurate er- 
ror estimate for each cluster is implausible. To assess 
the typical uncertainty, we analyzed WFPC2 images 
in F555W and F606W of the central regions of M87 
th at overlap wi t h the published ACS GC size catalog 
of I Jordan etaU (|2009D . F555W is the filter for which 
the charge diffusion is best characterized; F606W is the 
dominant filter in our catalog. These WFPC2 images 
have exposure times close to those of the ACS images, 
and the comparison objects have a distribution of mag- 
nitudes similar to those of our catalog . 

Considering the Uordan et al.l (|2009() sizes themselves: 
there are measurements in two different bands (F475W 
and F850LP), with a systematic difference of ~ 0.25 pc 
betw een the niters (F850LP sizes are larger; [Jordan et al.l 
I2005T ) . The typical reported uncertainty among the 
bright ACS GCs is 0.15-0.20 pc. We use a straight av- 
erage of the size estimates in both filters fo r comparison 
to the WFPC2 sizes. Unlike for our sizes, Uordan et al.l 
(2009) do not assume a fixed King model concentration, 
but (i) moderate variations in the concentration affect 
the half-light radius very little, and (ii) the S/N of the 
majority of GCs is insufficient to constrain this parame- 
ter in any case. 

The median difference between our F555W sizes and 
those in Jordan et al. is —0.70 pc, with a standard 
deviation of ~ 0.5 pc. In F606W, the median differ- 
ence is +0.17 pc, but with a nearly identical dispersion. 
Given the small uncertainties of the ACS sizes them- 
selves, this suggests that the typical random errors for 
our WFPC2 sizes are 0.4- 0.5 pc (about 15% for a normal 
GCV lLarsen et all (|2001| ) suggested approximate uncer- 
tainties of 1 pc for WFPC2 data of GCs in Virgo cluster 
galaxies, but the bulk of their objects were fainter than 
ours. 

More troubling is the large systematic difference be- 
tween the WFPC2 F555W and F606W half-light radii. 
While the ACS sizes are not free from error — and indeed 
the filter-to-filter size variatio ns indicate the pr esence 
of systematic effects (see also iSpitler et al.l 120061 ) — the 
smaller pixels and better-sampled PSF suggest naively 
that ACS sizes are probably more accurate than those 
from WFPC2. Thus it is surprising that the difference 
from "truth" is largest in the F555W filter, for which 
the effects of charge diffusion are expected to be best un- 
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derstoocFI. These results suggest that WFPC2 charge 
diffusion may be poorly characterized, even in F555W. 
To our knowledge no other direct comparisons between 
WFPC2 and ACS sizes exist. In any case, the bulk of 
our sizes are measured in F606W, which seems to agree 
best with the ACS effective radii, and we do not make 
any corrections to our measured sizes. 

Table [7] lists our size estimates, along with the instru- 
ment and filter, for the interested reader. In total, we 
estimated half-light radii for 197 objects that are spectro- 
scopic members of M870- Thus 344 clusters — nearly half 
of the spectroscopic catalog — have sizes. We have also 
verified that an additional 39 objects with low radial ve- 
locities are unresolved foreground stars; these are marked 
as such in Table [8] This is by far the largest publicly- 
available data set of GCs around any galaxy with both 
measured sizes and spectroscopic confirmation, and en- 
ables us to clarify the relation s between UCDs and GCs 
(gXU and IBrodie et aLllMTl . 

4.4. Object Classification 

One of the key ingredients in our kinematical and dy- 
namical analyses is to have a clean sample of very-high- 
probability GCs. To derive this sample, we use a combi- 
nation of colors, luminosities, sizes, and velocities to dis- 
tinguish GCs from other "contaminating" objects. The 
two obvious classes of interlopers are background galaxies 
and Galactic foreground stars; a third, more subtle group 
are the UCDs, which inhabit a similar color-magnitude 
space to the GCs. 

4.4.1. Ultra-compact dwarfs 

Considering the UCDs first: the largest objects, with 
sizes approaching ~ 100 pc, are almost surely a dif- 
ferent class of objects than "normal" GCs, and can be 
considered bona fide galaxies or UCDs. However, there 
is a substantial gray area of smaller objects that could 
include both the most compact UCDs and the most ex- 
tended GCs. True UCDs are unlikely to share the kine- 
matical properties of normal GCs and so should be ex- 
cluded from the pure GC sample. 

Distinguishing between UCDs and GCs is not a prob- 
lem we will solve in this paper, but we will try to at 
least provide some clear boundaries for the case of M87. 
In Figure [3] we show a diagnostic diagram of size ver- 
sus magnitude. Considering both the sample of objects 
imaged with HST/ ACS in the central regions of M87 
(Jordan et al.) and our spectroscopic sample with sizes 
over a wider field ( 34.3[) . there is a clear ridge-line of 
GCs centered at ~ 2.5 pc over most of the luminosity 
range, with a gradual upturn to ~ 3.0-3.5 pc beginning 
at magnitudes brighter than iq ~ 21. At these bright 
magnitudes, there is also a broad tail of objects extend- 
ing to very large s i zes. Th ese general observations paral- 
lel those of lHarrisI (I2009bfl for GCs in a combined sample 
of HST/ ACS BCGs, although there appear to be fewer 
large objects in the Harris sample than in M87; whether 

17 Convolution with an older, more diffuse kernel, which was at 
one time recommended by the WFPC2 Data Handbook, produces 
an even larger discrepancy with the ACS sizes. 

18 A single object, H38352, is unresolved in its WFPC2 images. 
Its radial velocity of ~ 1100 km s _1 confirms M87 membership; 
we estimate an upper limit of ~ 1 pc for its effective radius. 
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Fig. 3. — Half-light radius versus i-band magnitude for com- 
pact stellar systems. Small gray points show GC candidates in the 
central HST/ ACS field Uordan et alj l2009). while larger symbols 
show objects in the spectroscopic sample. Dotted lines show our 
diagnostic boundaries for GCs, UCDs, and uncertain/transition 
objects, with accordingly different symbols for the spectroscopic 
objects in these regions. Spectroscopically-confirmed objects with 
no size measurements are plotted as tick-marks near = 0. The 
vertical dashed line shows the approximate GCLF turnover mag- 
nitude. A histogram on the right-hand side shows the number of 
objects in size bins. There are another 20 UCDs with > 20 pc 
that are not visible on this plot. 

this difference is real (perhaps because such objects are 
preferentially found at large radii) or due to selection 
effects is unclear. 

Following the envelope of GC sizes from faint to bright 
magnitudes, we identify ~ 5.25 pc as an approx- 
imate upper boundary to the "normal" GC distribu- 
tion. To classify objects as UCDs, we use a boundary of 
rh > 9.5 pc, which is somewhat arbitrary but consistent 
with findings in other environments from t he literature 
(e.g.. lForbeslTaII[2fJol iTavlor et alJl2010h . Using this 
criterion, there are 34 UCDs that are confirmed mem- 
bers of M87, including 18 that we have newly idenfitied. 
Objects in the intermediate size range of rh ~ 5-10 pc 
may include a mix of GCs and UCDs, and for now we 
designate them as transition objectJ^l. 

We illustrate our adopted size boundaries in Figure [3l 
We emphasize that although many previous studies of 
UCDs have included a minimum luminosity as a classi- 
fication criterion (e.g., Mj < —11), the low luminosities 
of UCDs that we have discovered here demonstrate that 
such a bou ndary is an artifac t of observational selection 
effects (see lBrodie et al.ll2011l for a fuller discussion). At 
very faint magnitudes, Mj > —9, M87 may harbor a pop- 
ulation of extended clusters as found in the Milky Way. 
This could include some bona fide GCs that have under- 

19 A notable Local Group object in the transition region is u Cen 
(rh ~ 7 ~8 P c )> while M54 and Gl (with r h ~ 3-4 pc) would be 
considered "normal" GCs. NGC 2419 would also be classified as a 
UCD with our criteria. 
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Fig. 4. — Color-magnitude diagram of objects around M87. Gray points are all point sources in the CFHT field-of-view, while the 
other symbols are spectroscopic objects identified as GCs, UCDs, transition objects, foreground stars, background galaxies, and ambiguous 
objects (in the ~ 150—350 km s — 1 range), as identified in the legend at upper left. All objects with spectroscopic data, whether old or 
new, are included here. Note that a few of the GCs were outside the CFHT field-of-view, so their photometry shown here is based on 
approximate conversions from the Subaru BVI photometry; such conversions for background galaxies are less accurate, and we only show 
cases with direct gri photometry. Typical photometric uncertainties are shown by error bars near the left-hand axis. 



gone relaxation-driven expansion, but our spectroscopic 
sample does not go faint enough for this complication to 
be a worry. 

Even the "pure" samples may include some misclassi- 
fied objects (because there may not be clean boundaries 
in principle, and in practice our sizes are only accurate 
at the ~ 10-15% level anyway), but these should be very 
few and would not compromise our analyses in general. 
Slightly more problematic is the large sample of spectro- 
scopic objects without measured sizes, as these will in- 
clude a substantial population of of unrecognized UCDs 
and transition objects. We estimate this "contamina- 
tion" fraction will be ~ 50% at the bright end (iq < 20), 
and ~ 10% for fainter objects, and so we will generally 
excluded the brighter objects from our kinematical anal- 
yses. 

Color and magnitude diagnostic diagrams are shown 
in Figures 2] and [5j An important result from these di- 
agrams is the color distribution of the UCDs and tran- 
sition objects. These are mostly drawn from the blue 
side of the color distributions, and they appear to fol- 



low a coherent color-magnitude relation that is tilted 
and distinct from the g eneral blue GC subpopulation. 
lEvstigneeva et all (|2008l ) show that the nuclei of Virgo 
dwarf elliptical galaxies follow a similar offset relation, 
with the brightest Virgo UCDs falling at the upper end 
of said relation, as possible evidence of the "threshing" 
hypothesis for UCDs. Our observation that UCDs and 
transition objects have a similar color-luminosity rela- 
tion provides evidence that most of the transition objects 
are in fact an extension of the UCD population rather 
than being extended GCs, although this conclusion is 
necessarily preliminary. 

A more in-depth discussion of the M87 U CDs and tran- 
sition objects is covered in another paper (|Brodie et al.1 
1201 lh . One additional comment here is that Fig- 
ure [3] suggests that the size-luminosity correlation for 
UCDs is not as strong as conventionally thought (e.g., 
lEvstigneeva et al.|[2008h . 

4.4.2. Foreground Stars 
ICote et all (|2001f ) used a combination of C — T\ color 
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Fig. 5. — Color-color diagnostic diagram of objects around M87. 
Symbols are as in Figure[4] Dashed lines show GC color boundaries 
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and radial velocity to separate Galactic foreground stars 
from GCs. All objects with v > 200 km s _1 and a 
color in the range 0.8 < C — T\ < 2.35 were classified 
as GCs, while objects with v < 200 km s _1 (regardless 
of color) were classified as stars. In this paper we also use 
a combination of photometry and radial velocity to select 
GCs, with the addition of high-resolution HST imaging 
where available. We s upplement the Subaru /Suprime- 
Cam photometry from iTamura et al.l (l2006af) wi th gri 



CFHT/Megacam photo metry (see iHarrisi l2009at 



and gz HST/ ACS data (j Jordan et al.l I2009D when pos- 
sible; the original H+01 CT\ photometry is also used. 

Since M87 has a systemic velocity of 1307 ± 8 km s _1 
(|Smith et al.l 12000ft and the GCs have a large velocity 
dispersion (~ 350-400 km s" 1 ; ICote et alj|2001l) it is 
likely that some true M87 GCs will have radial velocities 
comparable to foreground stars. If we assume that the 
Galactic halo is purely pressure supported, the stellar 
velocity distribution should peak at a heliocentric radial 
velocity of 40 km s _1 (assuming Ulsr = 244 km s _1 ; 
I Bow et al.1 l2009f ) with an approximate dispersion of 
120 km s _1 . Therefore, foreground stars with velocities 
larger than ~ 300 km s" 1 should be uncommon, and 
those > 500 km s _1 nearly nonexistent. Expectations 
for M87 are less clear: the velocity dispersion of the GC 
system is large and varies with radius, substantial rota- 
tion could be present, and, as we discuss below, there 
is compelling evidence for kinematical substructure. Es- 
pecially for the latter reason, we do not perform a full 
likelihood analysis to assess M87 membership probability 
that would rely on the assumption of a Gaussian velocity 
distribution. Nonetheless, objects with very low radial 
velocities (< 150 km s _1 ) are much more likely to be 
foreground stars than GCs. 

We initially consider all objects with v > 350 km s _1 
to be GCs and all objects with v < 150 km s _1 to be 



stars, unless there is strong evidence to the contrary. Un- 
fortunately this (unavoidable) selection makes it more 
difficult to identify some of the most obvious intergalac- 
tic GCs, as we will discuss later. There are three objects 
(from the older spectroscopic data sets) apparently below 
this lower boundary that are HST confirmed GCs: S7017 
(v = 94±101 kms- 1 ), S7006 (v = 104±58 kms" 1 ), and 
S7007 (v = 138 ± 94 km s" 1 ). HST imaging is among 
the most trustworthy methods for separating stars from 
GCs, since nearly all luminous M87 GCs are reliably re- 
solved with WFPC2 or ACS images of typical duration. 
A single .ff5T-confirmcd foreground star (S161) lies above 
the upper velocity limit of 350 km s , but it has a very 
large velocity uncertainty (v — 446 ± 246 km s _1 ). 

There are 50 objects in the combined catalog in the 
most ambiguous range, 150 < v < 350 km s _1 . Of these, 
11 have certain classifications from HST imaging: three 
are true GCs, and the other eight stars. Of the remain- 
ing 39 objects, the highest radial velocity is 289 km s _1 
(even though the sample is bounded at 350 km s _1 ), 
suggesting that the bulk of these objects are indeed fore- 
ground stars. The color-color diagram (Figure [5]) gives 
supporting evidence. The locus of spectroscopically con- 
firmed stars begins blueward of the GC sequence, crosses 
this sequence in the region of metal-poor GCs, and ex- 
tends far to the red. Many ambiguous objects which lie 
far from the GC sequence can be conclusively identified 
as foreground stars (and are marked as such in the di- 
agrams). However, it is evident that halo dwarfs in the 
approximate color range 0.5 < (g — r)o ;$ 0.6 cannot be 
cleanly separated from GCs with the existing photom- 
etry. For this paper, we classify all of these uncertain 
objects (marked as open purple pentagons in Figures [4] 
and [5]) as stars, but note that a subset may be true GCs. 

Two objects in the velocity range 350 < v < 
500 km s -1 do not have HST imaging. S1644 (v = 
496 ±68 km s _1 ) is in the region in gri color-color space 
occupied by metal-rich GCs and we classify it as a GC. 
H58443 (v = 441 ±24 km s _1 ) has colors consistent with 
either a very metal-poor GC or a foreground star and is 
likely to be a GC. 

The data show some unusual features outside of the 
ambiguous velocity range discussed above. Two of the 
reddest objects in the Hectospec sample (H66419 and 
H20573) also have low radial velocities for M87 mem- 
bers: 596±11 kms -1 and 683±10 kms" 1 , respectively. 
Both objects have high signal-to-noise spectra and the 
velocities are not in doubt. However, these radial veloci- 
ties are close to or exceeding Galactic escape velocity at 
their likely (but unmeasured) distances. Therefore, it is 
unlikely that they are foreground stars. 

As an aid to dynamical analyses and future observa- 
tions, in Table |5] we have compiled a list of the probable 
spectroscopic Galactic for eground stars. Thi s table does 
not include o bjects from iCohen fc Rvzhovl (|1997l ) and 
iCohenl (|2000[ ) that were identified as stars but without 
published radial velocities (S178, S221, S246, S354, S681, 
S996, S1610). These objects all have colors very dissimi- 
lar from those of GCs and are unlikely to be re-observed 
m the future 



20 Some other surveys, such as Firth et al. (2008), also unfor- 
tunately did not publish their lists of foreground stars and so are 
likewise excluded. 
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4.4.3. Cleaning the Catalog 

A small number of objects from H+01 have velocities 
consistent with being GCs but with photometry or other 
information that contradicts this classification. We dis- 
cuss these objects in turn. 

The object S923 (mentioned in ^4.2[) was previously ob- 
served by H+01 with a reported v — 1943 ± 134 km s _1 . 
Our new LRIS velocity is catastrophically different at 
2777±26 km s _1 , which would place it among the highest 
reported for an M87 GC. The spectrum is of high quality 
and the object is clearly very metal-poor, consistent with 
its blue color (g — r — 0.49). An archival WFPC2 image 
suggests an unusual, asymmetric morphology, with a re- 
solved core (of size ~ 2 pc) and two small protuberances. 
Given the morphology and the outlying velocity, we ex- 
clude this object from our analysis; better HST imaging 
is needed for this very interesting candidate. 

S7023 presents an unusual case. The object closest 
to the position reported by H+01 is clearly resolved 
in the CFHT image; it is also present in an archival 
ACS image and has a tight core with a diffuse enve- 
lope. This morphology and the outlying radial veloc- 
ity (415+ 134 km s _1 ) are generally consistent with the 
idea that it is an ultra-compact dwarf galaxy in the Virgo 
cluster with a half-light radius of tens of pc. However, its 
g — r and g — i colors are too red for a normal old stellar 
population — and, in addition, are inconsistent with the 
reported C — T\ color. The nearest potential interloper is 
located 4" from the listed position and is a background 
galaxy. We exclude this object from our analysis but 
further investigation is warranted. Both S923 and S7023 
are listed at the end of Table [§] 

S1602 has v ~ 1010 km s" 1 , but a C - T x color of 
3.42. The CFHT 5 -band image of S1602 shows that it 
is a close pair; both components have colors too red for 
a GC, and morphologies consistent with those of back- 
ground galaxies, so we exclude it. In addition, S7008 
and S7012, which have H+01 velocities nominally con- 
sistent with M87 membership, are also apparent back- 
ground galaxies. We do not know if their identifications 
are incorrect or if these are additional catastrophic radial 
velocity failures; in any case, they are removed from the 
sample of GCs. S508, with a large reported velocity in 
H+01 (v ~ 5800 km s" 1 ), is another clear background 
galaxy. 

Finally, we have resolved a subset of the ambiguous 
object identifications in H+01. For example, in the close 
pair S8052, one of the objects is a background galaxy, so 
the targeted GC is not uncertain. In the pair S8051, only 
one of the objects has gri colors consistent with that of 
a GC. S9001 and S9002 are both close pairs for which 
one (ambiguous) member has a published radial veloc- 
ity. We have obtained new velocities for one member 
of each pair; each disagrees with the published velocity. 
In the absence of other information, we assume that the 
new observations cover novel objects, and so we can now 
correctly identify S9001 and S9002. 

4.4.4. GCs in Dwarf Galaxies? 

In our quest for a pure catalog of M87 GCs, we also flag 
any objects that have a reasonable chance of being bound 
to one of the dwarf galaxies near to M87 (NGC 4476, 
NGC 4478, NGC 4486A, NGC 4486B, IC 3443). To do 



this, we follow a similar approach to iFirth et al.l ([2007D 
and define a tidal radius i?t for each dwarf: 

where R p is its distance from M87, is its mass, and 
M r is the mass of M87 enclosed within r = R p . Both R t 
and R p should in principle be 3D radii, and the percen- 
ter rather than instantaneous radius should be used, but 
we will assume the observed projected radii work well 
enough. If we do identify an object within a dwarf tidal 
radius, we examine its velocity relative to the dwarf, and 
compare this to the escape velocity, which we define as 
v CBC ~ ^/2v c , where v c is the local circular velocity of the 
dwarf. If the object is still consistent with being dwarf- 
bound, we remove it from the overall kinematics sample. 

To estimate the values, we construct simplified 
models for the dwarfs based on spectroscopic and photo- 
metric data from the literature , and for M r , we use the 
dynamical results for M87 from lMurphv et al.l (|201lD . In 
the end, we identify one object which may be bound to 
NGC 4478: H35970. 

4.5. The Final Catalog 

Accounting for these classifications, and combining the 
new velocities with those from the literature, gives a to- 
tal sample of 737 objects associated with M87. 707 fall 
into our pure GC or transition samples, while the re- 
mainder are UCDs or other extended objects. The final 
catalog, including all potential M87 members, is given in 
Table [9] Where relevant, the listed radial velocity is a 
weighted average of the available values. The gri pho- 
tometry is listed in this table; all other photometry has 
been included in a supplementary table (Table [TU]), with 
the exception of the removal of some CT\ measurements 
that strongly disagreed with other data. Except for the 
H+01 CT\ data, the photometry is corrected for fore- 
ground reddening, as described here or in the relevant 
papers PI. 

Cluster astrometry in this catalog was recalibrated us- 
ing the CFHT g-band image and USNO coordinates for 
all objects within the field of view. These coordinates 
have proven sufficiently accurate for reliable multi-object 
spectroscopy. A few sources from the literature could not 
be uniquely identified or fell on image defects and so re- 
tain their previo us positions. I n addit ion, a number of 
objects from the iTamura et alj (|2006aD catalog lie out- 
side of the field of view of the CFHT image. 

The profusion of papers on GCs and UCDs in M87 has 
led to a corresponding proliferation in object IDs. We 
have attempted to keep our naming transparent while 
generally retaining principal IDs used in refereed publi- 
cations withspectra. All objects within the field of view 
of the lHarrisI (|2009aD CFHT imaging have been given 
"H" IDs and these are listed as alternate IDs where ap- 
plicable. 

21 F or the reader who wishes to make the correction, Co te et al.l 
(2001) report a mean reddening of E(C — Ti) = 0.045 mag in the 
direction of M87. Note also that there were some complications 
in the photometric ca libration of the Suprime-Cam BVI data (see 
Tamura et al. 2006a). Comparing to the CFHT gri data, we find 
chip-level zeropoint variations of typically 0.1 mag for a subset of 
the Suprime-Cam chips; one should keep this caveat in mind when 
making use of the BVI data. 
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4.6. Globular Cluster Subpopulation Classification 

Although one of the goals if this paper is to examine 
afresh the classification of M87 GCs into independent 
subpopulations, as a starting point we first adopt a clas- 
sical bimodal separation into two components: a blue 
(metal-poor) subpopulation and a red (metal-rich) sub- 
population. A detailed analysis of stellar populations 
from a subset of our spectra will be the subject of a 
follow-up paper, but in an initial inspection of the spec- 
tra, we see no sign of GCs younger than ~ 2 Gyr. This 
further motivates the interpretation of the color varia- 
tions as tracing metallicity differences for old popula- 
tions. 

We classify objects into subpopulations using as many 
independent photometric estimates as possible. For C — 
Ti and g — z a simple color cut is used: C — T\ = 1.42 
and g — z = 1.20, respectively. For GCs with BVI and 
gri, we empirically define the GC sequence in color-color 
space and separate clusters using a line orthogonal to 
the GC sequence and passing through V — I = 1.10 and 
g — i = 0.93, respectively. These boundary values are 
suggested by mixture modeling of the color distr i butions 
and are identical t o those used in iTamura et al.1 ([2006b) 
and Harris ( 2009a), respectively. 

For 92% of the GCs, the colors agree uniformly on the 
classification, or there is only a single color. Another 6% 
(typically GCs close to the middle of the color distribu- 
tion) have at least three independent sets of photometry 
that do not all agree; we choose the classification favored 
by two or more sources. Fourteen GCs of the total sam- 
ple of 707 have equivocal evidence for classification. In 
these cases we use, in order of priority, C — T±, g — z, gri, 
and BVI. 

Using this scheme, 486 (69%) of the objects are classi- 
fied as blue metal-poor GCs, with the remainder as red 
metal-rich clusters. For completeness, the UCD candi- 
dates have also been classified, but are not included in 
the above statistics. 

5. PHOTOMETRIC PROPERTIES OF THE GLOBULAR 
CLUSTER SYSTEM 



Although Harris! (|2009aj) has already analyzed the 
CFHT imaging data in detail, for consistency we repeat 
some of his analyses, and include various extensions of 
relevance to our overall photometric and spectroscopic 
study. In §5.1 we fit the radial surface density profiles 
of GCs using Sersic functions. §5.2 discusses contamina- 
tion, and §5.3 the ellipticity of the GC system. In §5.4, 
aided by spectroscopy, we examine trends of the GC col- 
ors and subpopulations. 

5.1. Radial Density Profiles 

The number density distribution of GCs with radius 
is important for a number of reasons, including as in- 
put to dynamical modeling of the GC system. The 
three-dimensional distribution is inferred from deprojec- 
tion of the surface density, which for an extended system 
like M87's demands wide-field imaging. The metal-poor 
and metal-rich GC subsystems can be analyzed indepen- 
dently, and require derivations of their density profiles 
separately. 

Several previous studies have derived surface density 
profiles for the separate GC subpopulations in M87, 



including I Cote et al.l (120011) ITamura et al.l (|2006bl) , 
iForte et all (|2007f) . and HarrisI (|2009aD . 

Nearly 700 spectroscopically confirmed M87 GCs are 
present in the CFHT photometric catalog. As shown in 
Figure El these empirically define the GC sequence in 
color-color space. The GCs scatter around a mean locus 
extending approximately from (g — r, r — i) = (0.48, 0.18) 
to (0.78, 0.43). We experimented with using rather strin- 
gent criteria in color-color space to select GCs, but found 
that the slight advantage in rejection of contaminants did 
not overcome the loss of many true clusters. Instead, we 
used the relati vely loose color cuts 'm g — r and g — i from 
Harris (2009a). Objects were selected over the magni- 
tude range 19 < i < 22.5, with the faint limit chosen to 
avoid incompleteness issues relatively close to the center 
of the galaxy. Extended galaxies were excluded using a 
cut on <?-band FWHM. Metal-poor and metal-rich GCs 
were separated using a line perpendicular to the mean 
color locus and passing through g — i = 0.93. 

We chose to fit the GC surface density profiles with 
iSersid (|1968f ) functions for several reasons: (i) these have 
deprojected analytic forms suitable for use in dynamical 
modeling; (ii) they extrapolate to finite total GC num- 
bers; and (iii) both the observed stellar components of 
galaxies and their simulated dark matter halos are repro- 
duced well by Sersic-like functions, so it is only natural 
to suppose that the GC systems are as well. 

The GC counts were binned in projected radius over 
the range l'-29'. Corrections were made in the binned 
profile for the chip gaps present in the CFHT images. We 
assumed Poissonian errors both from the GC counts and, 
to be conservative, from a constant contaminant level of 
0.5 arcmin -2 . We also assumed a systematic uncertainty 
in the contaminant level, set to 0.15 arcmin -2 for the full 
sample and 0.10 arcmin -2 for the individual blue and red 
subpopulations. 

We fit a Sersic function of the form: 
No exp[— (R/R s y/ m ], with normalization No, scale 
radius R s , and index m. The contamination level was 
left to vary freely. Note that "contaminant" refers to any 
objects not part of the GC system of M87 — these may 
be foreground stars, background galaxies, or potentially 
even intergalactic GCs with a much more extended 
spatial distribution than that of M87. Fits were done 
for the full sample and blue and red subpopulations 
separately. 

The best-fit contaminant levels are 0.51 ± 0.16, 0.18 ± 
0.17, and 0.36 ±0.09 arcmin -2 for the full, blue, and red 
samples re spectiv e ly. The se are consistent with previous 
estimates: lHarrisI (|2009af ) found background densities of 
0.26 and 0.36 a rcmin -2 for th e blue a nd red objects, re- 
spectively, and ITamura et al.l (|2006bf ) reported a value 
of ~ 0.2 arcmin -2 for blue "intergalactic" GCs (IGCs). 
Note, however, that the very different luminosity func- 
tions of GCs and contaminants imply that the inferred 
density of interlopers will strongly depend on the selec- 
tion function of candidates. 

The Sersic fits are plotted, together with the surface 
density estimates, in Figure El As is t ypical for GC sys - 
tems (and well-established for M87; ICote et al.l 120011 ). 
the metal-poor GCs are more extended than the metal- 
rich GCs, with the crossover point at a projected radius 
of ~ 8-10 kpc. Out to the edge of our photometric data 
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Fig. 6. — Surface density profiles for GCs in M87. The full sam- 
ple of GCs is given by the black circles, with the blue and red 
subpopulations plotted as blue squares and red triangles respec- 
tively. The Sersic fits are overplotted in the appropriate colors 
as solid, long-dashed, and short-dashed lines (f or full, blue, and 
red sa mples). The integrated light profile from Kormcndy ct al. 
(2009) is plotted as a dot-dashed line, and is a good match to the 
red GC profile over an extended radial range. The GC profiles are 
constructed in circular bins, and the galaxy light is plotted versus 
circular-equivalent radius: rgMA(l — e) 1//2 . 

at ~ 30', we see no evidence for either a sharp trunca- 
tion of the outer density profiles, nor for a flattening that 
could mark the transition to a n intergalactic or a n ob- 
viously accreted component (cf lHuxor et al.ll2~01lD . 38.41 
contains further discussion of this issue. 

Figure [5] also plots t he galaxy light distribution from 
iKormendy et all (|2009t ) . While the normalization of this 
M87 profile is arbitrary, it is clearly a good match to 
the surface density o f the metal-rich GCs over a wide 
radial range (see also lTamura et aT]|2006b| ). The devia- 
tions occur in the central few kpc, where the GCs have 
a flat core-like distribution, and in the very outer halo 
(> 100 kpc). This lends credence to the generic idea 
that the metal-rich GCs formed along with the bulk of 
the field s tar population in early- type galaxies (see the 
review in Brodic & Strader 2006). The normalization 
factor is such that the inferred mass in metal-rich GCs 
is ~ 0.28% of M87 over the matching range in radius. 

The derived parameters of the surface density fits are 
given in Table [TTJ (except for the background levels, 
which are given in Table [T2"|) . Table [TT] also contains the 
results of fits for which the background level is fixed to 
the spectroscopic estimate determined in 35.21 and listed 
in Table Q2] It is challenging to assess the uncertainties 
on these parameters, since they are strongly coupled; the 
difference between the fixed and varying background fits 
gives a general sense of the constraints. Nonetheless, 
a range of values would likely provide reasonable fits. 
For example, modulo norma l ization facto rs, the R 1 ^ 4 fits 
presented in iTamura et ail (|2006bD and Harris! (|2009aft 
(equivalent to m — 4 Sersic functions) generally match 



the data well. 

For reference, the implied globular cluster system 
(GCS) parameters for the fits with free background lev- 
els are as follows. The GCS effective radius is ~ 6.8', 
22', 23' (~ 32, 106, 110 kpc) for the red subpopulation, 
blue subpopulation, and total system, respectively, with 
corresponding Sersic indices of ~ 5.3, 3.7, and 5.5. The 
respective total number of GCs is ~ 2000, 5500, 8900, 
although nearly half of the total GC population inferred 
from these fits is beyond the edge of our photometric field 
of view, and these values exclude the faint end of the GC 
luminosity function (with i > 22.5). 

5.2. Spectroscopic Estimates of the Contamination 

Our Sersic fits left the contaminant level as a free 
parameter. Here we combine the original, large radius 
DEIMOS spectroscopy and CFHT imaging to give an 
independent estimate of the density of interlopers. 

The radial overlap between the DEIMOS and CFHT 
data begins at ~ 8' from the center of M87 and extends 
to the edge of the CFHT field at ~ 30'. Of the objects 
classified with DEIMOS and with CFHT photometry, 
there are 54 confirmed GCs and 44 confirmed contami- 
nants. We then applied the selection function used for 
the surface density estimates in tj5.1l to both the GC and 
contaminant samples, yielding 46 GCs and 20 contami- 
nants. These samples were binned in radius as for the 
surface density estimates. This procedure allows multi- 
ple independent estimates of the contamination fraction, 
albeit with large individual uncertainties because of the 
small numbers of objects. The final surface density of 
contaminants in each radial bin is calculated by multiply- 
ing the local contaminant fraction by the overall surface 
density in that bin. The same method can be used on 
the blue and red GC subpopulations separately, though 
with necessarily larger uncertainties. 

Table Q2] compares the spectroscopic and photometric 
estimates of the contaminant level for both the full sam- 
ple of GCs and for the blue and red GCs. The listed 
spectroscopic values are a weighted average of the in- 
dependent radial bins. There is excellent agreement be- 
tween the photometric and spectroscopic estimates (and, 
indeed, even better agreement between th e spectroscopi c 
estimates and the photometric values of lHarrisI 12009a) . 
These comparisons suggest that our values are accurate 
estimates of the background level for the assumed selec- 
tion function. 

5.3. Ellipticity 

Most kinematical and dynamical analyses using dis- 
crete tracers assume spherical symmetry for simplicity. 
This assumption could lead to systematic errors in the 
mass modeling of real early-type galaxies, which are 
rarely spherical. M87 is notably aspherical, with its stel- 
lar isophotal ellipticity increasing f r om ~ 0.2 at 3' to 
~ 0.45 at 13' (jKormendv et al.ll2009ft . iMcLaughlin et al.l 
(|1994| ) estimated the ellipticity of the GC system in the 
inner regions of M87, finding it roughly consistent with 
that of the galaxy light. 

In this subsection we estimate the ellipticity of the 
GC system of M87, both as a whole and for the 
blue and red subpopulations. Such a subpopulation 
analysis has been done only rarely in the literature 
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Fig. 7. — Ellipticity profiles for the M87 GC system. Symbols are 
as in Figure [6] includin g the integrated light profile of M87 from 
IKorme ndv ct al. (2009). Circular-equivalent radii are plotted. 

on early- type galaxies (e.g., Kissler-Patig et all 119971) . 
A number of methods have been used in the litera- 
ture to determine ellipticities for sparsely sampled two- 
dimensional data; here we use the method of mom ents 
(|Carter fc Metcalfe] fl980l: McLaughlin et alJ [l99l . In 
general, maximum likelihood is a preferable approach, 
but we found th at maximum likelihood e stimates using 
the formalism of [McLaughli n et al.l (|1995f ) , when applied 
to radially binned data, were frequently unstable. 

We bin the data into four unequal radial bins in the 
range l'-13', with the outer limit set by incomplete spa- 
tial coverage due to the chip gaps. Both the ellipticity 
and position angle are allowed to vary. The derived el- 
lipticity values are plotted in Figure [71 along with th e 
circularized galaxy light from IKormendv et al.l (|2009f) . 
The ellipticity profile of the GC system as a whole ap- 
pears somewhat more constant than that of the galaxy 
light. The red GCs have a steeper ellipticity profile, out 
to the last measurement where their surface density be- 
comes comparable to that of the background and thus 
very uncertain. The ellipticity profile for red GCs has a 
similar slope to that of the galaxy light. By contrast, the 
blue GCs have a nearly constant ellipticity profile, con- 
sistent with e ~ 0.29 at all radii. The ellipticity values 
for the full GC s ample in the inner region s agree with 
those obtained by IMcLaughlin et al.l (| 19941) . 

The ellipticity and position angle estimates are given 
in Table US 

5.4. Beyond Bimodality 

The preceding analyses have operated under the sim- 
plifying assumption that the M87 GCs fall into two stan- 
dard subpopulations (red and blue) whose identities are 
invariant with respect to GC luminosity and galactocen- 
tric distance. However, the richness of our data set al- 
lows us to explore the extent to which this simple sce- 
nario is valid. Below we examine color trends with radius 



5.4.1. Color Gradients 

We first examine the color distribution of GCs versus 
distance R p in Figure [5J which shows the (g — i)o col- 
ors for the full photometric sample of GCs and for the 
spectroscopic objects which pass the photometric selec- 
tion criteria (95% of the total sample). As the bulk of 
the spectroscopic objects were selected without regard to 
color (especially our new LRIS and Hectospec samples), 
the bivariate distribution of spectroscopic objects should 
be a reasonable representation of the full set of GCs, as 
it generally appears to be. The reader can look ahead 
to Figure [5] to sec density plots of (g — i)o for both pho- 
tometric and spectroscopic GC samples broken down by 
projected radius and magnitude. 

The plot of the full radial range shows several remark- 
able features. Although the stated division between blue 
and red GCs appears reasonable in the inner region, in 
the radial range ~ 7-14' (~ 35-75 kpc) the distribution 
of red GC colors appears to change: the peak eithe r shifts 
bluewa rd (as assumed in the power-law models of Harris 
l2009al — see his section 6) or is replaced by a population of 
intermediate-color objects. Intermediate-color GCs were 
in fact reported in t he central regi ons using HST/ A CS 
(| Strader et alJ 2006; also visible in iPeng et all 120061 as 
residuals from their bimodal model), and the new com- 
ponent at larger radius may be an extension of this inner 
"subpopulation" . An inspection of the inner radial plot 
shows that the overall red GC color transition may start 
at a projected radius of ~ 4' (~ 20 kpc), outside of which 
the principal red (g — i)o edge of the metal-rich GCs be- 
gins to shift. 

At large radii, the blue GCs continue to have a well- 
defined and roughly con stant-color pea k, but with no- 
ticeable radial variations. [Harris (2009a) reported a color 
gradient in the blue GCs within a radius of ~ 12', and 
with a slope of —0.025 versus log radius. We can ap- 
proximately reproduce this result through a linear least- 
squares fit to the photometric sample of blue GCs. How- 
ever, we can now see that the fit is clearly affected by 
the presence of intermediate-color objects at small radii, 
and the evidence for a slope at radii > 3' is marginal. 

To more accurately trace the mode of the blue GC dis- 
tribution, we defined a sample of fainter GCs (21 < ig < 
23) to eliminate a portion of the intermediate-color pop- 
ulation (which we find below to be a stronger component 
at brighter magnitudes), while retaining good statistics. 
Density estimates with a fixed Gaussian kernel of 0.03 
matjP^l were constructed in 1' radial bins. 

The blue mode of each bin is plotted in Figure [8j out 
to a radius of ~ 15' after which contamination of the 
photometric sample becomes significant. The "bump" in 
the color at ~ 8' appears real in that it reflects the actual 
color distribution of GCs (rather than just small number 
statistics). However, outside of a radius of ~ 3-5' (<~ 15- 
25 kpc), there is no evidence for a monotonic radial gra- 
dient among the blue GCs. This transition from an inner 
radial gradient to an essentially constant color h as also 
been seen in the GC systems of the Milky Way {Harris 

22 This value is close to the optimal Silverman (1986) kernel 
bandwidth for typical values of a and N for bins in the blue peak. 
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Fig. 8. — M87 GC colors versus galactocentric radius. Small black points show the photometric sample, while colored points show the 
spectroscopically-confirmcd sample, with blue and red colors denoting our primary color classification [using a division at (g — i)o = 0.93]. 
The cyan triangles show our estimates of the blue GCs peak location in radial bins using a photometric sample (see main text for details). 
The lower panel shows nearly the full radial range of confirmed GCs (out to 30'); the upper panel zooms in on the inner part of the cluster 
system within 9'. 



2001 1 ) and of the giant elliptical NGC 1407 (jForbes et al.l 
20111) . with implications that we will discuss in <J5J 

We also notice another color "bump" in the blue GCs 
at R p ~ 21-27' (~ 100-130 kpc), which is confirmed 
by the spectroscopic sample. These objects are spread 
out nearly uniformly in azimuth, so it is unlikely (but 
still possible) that Galactic cirrus on a smaller r esolu- 
tion than the iSchlegel. Finkbeiner. fc Davisl £l998) dust 
maps is contributing to this feature. Internal reddening 
gradients could also play a role in the color changes in 
principle, but would require invoking ad hoc dust distri- 
butiontQ 

As previously discussed in we emphasize that 

possible photometric gradients are dependent on accu- 
rate large scale photometric calibration, which we cannot 
guarantee. Therefore, these conclusions are still tenta- 

23 Reddening from dust (whether internal or external) would 
produce identical color shifts for both blue and red subpopula- 
tions, but this test is difficult to apply because the red GCs do 
not maintain a clear, contiguous color peak over a large range in 
radius. 



tive, pending the forthcoming photometry from NGVS. 
5.4.2. Color Subpopulations and Substructure 

An impression garnered from both panels is the pres- 
ence of substructure in the GC color distribution, es- 
pecially in the red GCs. No firm conclusions can be 
drawn from the spectroscopic data alone, but one has 
the impression of clumping in the metal-rich subpopula- 
tion, both at small radii and in the outer regions. 

To explore the red GC variations further, we con- 
structed density plots of the (g — i)q color of all photo- 
metric (not just spectroscopic) GCs. A fixed kernel size 
of 0.03 mag was used. Since our focus was on the red 
GCs, we binned in ellipses, starting with a semi-major 
axis of 1' and with the ellipticity equal to that of the 
galaxy light at that radius. The results for a range of 
1-5' in semi-major axis are shown in the left panel of 
Figure H] The blue GC subpopulation looks similar in 
each bin, with some evidence for a mild negative color 
gradient as already discussed. However, the distribution 
of red GCs varies significantly among the bins. There 
is no clear monotonic trend, and some of the bins (es- 
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Fig. 9. — Distributions of colors for M87 GCs, with a small amount of smoothing applied, and arbitrary normalizations. The left panel 
shows the central photometric sample in several radial bins as indicated in the panel legend. The middle panel shows the GCs with 
spectroscopic confirmation, out to larger distances, in several radial bins. The right panel shows the central photometric GCs binned by 
bright and faint magnitudes. The uncertainties in individual color measurements are generally < 0.02 mag. The gray shaded histogram 
in the left panel shows RGB data from Bird ct al. (2010), converting from metallicity to color using a relation from Sinnott ct al. (2010); 
note that the RGB data are not reliable for (g - i) > 1.12 ([Fe/H] > - 0.3). 



pecially the innermost one) show evidence for multiple 
subpopulations within the "red" peak. 

Small number statistics are very unlikely to be the 
cause of these variations. There are more than 100 red 
GCs in each bin (excepting the last one, which has 92). 
If we assume conservatively that the red GCs have an 
intrinsic color distribution with a = 0.1 in (g — i)o, then 
the standard error of the mean is < 0.01 mag per bin. 
However, the scaled median absolute deviatioro of the 
red peaks in these five bins is 0.048 mag. By contrast, 
the blue GCs show much more uniformity, and the me- 
dian absolute deviation of the blue peak color is 0.012 
mag. 

To extend this analysis to larger radii, we consider next 
the spectroscopic sample, which is virtually free of con- 
tamination. Note that later in this paper we will elect 
to use only our new data set because of the large impact 
on the kinematical and dynamical analyses of occasional 
catastrophic radial velocity errors in the older data, but 
such errors do not affect our photometric analyses and 
here we use the full combined sample. 

The middle panel of Figure |H] shows the color distribu- 
tion of the spectroscopic sample in several radial bins. At 
the small radii corresponding roughly to the ACS field-of- 
view (R p < 2.2' or < 10 kpc), the classic blue/red color 
bimodality is evident. At larger radii, the red subpopu- 
lation becomes less dominant as expected based o n th e 
relative radial distributions of blues and reds (see £)5.1[) . 
However, the average colors of the red distribution also 
become dramatically bluer, and the data are more sug- 
gestive of a shift between separate red and intermediate- 
color peaks than of a simple gradient of the red peak. 
This impression is supported by kinematic evidence in 
^gOlandlOl 

24 We use this robust estimator of the scale (which asymptoti- 
cally approaches a for a normal distribution) since the underlying 
distribution of the red mode is not necessarily Gaussian. 



This curious intermediate-color peak does not seem to 
be a manifestation of a very extended rcdwards tail of the 
blue sub population; if so, i t would be strong at all radii 
(see also iForte eFall 120091) . The third "peak" clearly 
has a spatial distribution that is more extended than the 
"classic" redder peak but less extended than the blue 
peak. For the sake of examining the third peak for un- 
usual kinematical behavior later in this paper, we carry 
out a trimodal fit to the photometric data over the radial 
range R ~ 1.5'— 4', using alternatively circular or ellipti- 
cal annuli. The three peaks inferred make roughly equal 
contributions to the GC numbers and we adopt color 
boundaries of (g - i) Q ~ 0.86 and - 1.010 

5.4.3. Bright GCs 

We next consider luminosity dependencies of the GC 
subpopulations. Separating out the photometric GCs 
with magnitudes in the range io=19-20.5, we show the 
color distribution for the small-to- intermediate-radius re- 
gions (0.5'-9' or 2-40 kpc) in the right panel of Figure^ 
along with the fainter GCs for comparison. The colors of 
both the blue and the red subpopulations are clearly dif- 
ferent between the bright and faint GCs. The blue shift 
is the well- known blue-tilt . The red shift was previously 
noticed by iHarrisi ()2009aD as a tentative and surprising 
"negative mass-metallicity relation", but we regard the 
trend as secure, and indicative of a third, intermediate- 
color subpopulation rather than of a red-tilt. 

Given the small number statistics for the bright ob- 
jects, we are not able to determine where exactly any 
magnitude boundary is between the bright and faint ob- 
jects (we have adopted io — 20.5 as a rough value). Nor 
are we able to determine much about the spatial dis- 
tributions of the bright objects except to note that the 
intermediate-color objects are somewhat more extended 

25 For example, for the circular fit the modes of the subpopula- 
tions are (g - i) = 0.76, 0.94, and 1.08. 
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in radius than the redder objects. We do see some pe- 
culiar indications that the locations of the color peaks 
for both bright and faint objects shift dramatically in- 
side radii of ~ 2' (~ 10 kpc). Our present data set is too 
limited to pursue this issue in detail, but we will find in 
£16.31 that the GC kinematics are also unusual inside this 
radius. 

For additional insights, we turn to the spectroscopic 
sample, where we find that the intermediate-color ob- 
jects, with (g — i)o ~ 0.85-0.95, include a larger fraction 
of bright (i < 19.5) and large (?*h > 10 pc) objects than 
are found at other colors, particularly at large radii. This 
result supports the idea that there is a third subpopula- 
tion with an unusual luminosity distribution, and further 
implies that this subpopulation is more strongly "con- 
taminated" with UCDs — although this fraction is still 
low (< 10%), with the fainter and smaller intermediate- 
color objects showing no unusual photometric properties. 

The possible convergence of the GC color-magnitude 
distribution to a unimodal distribution at the bright 
en d has been previously d iscuss e d in other galaxie s 
bvlOstrov. Forte, fc Geislerl (fl998h . iDirsch et al.l jl2003h . 
Harris et al.l (120061) . iRomanowskv et aLl (|2009). and 
Forbes et alJ (|2011[ ). Also relevant is the growing num- 
ber of early-type galaxies with centrally-concentrated, 
intermediate-color GC subpopulations without strongly 
unusu al luminosity functions. Exampl es include NGC 
4365 (L arsen. Brodie. fc Straden 120051 and references 
therein) and NGC 5846 (|Chies-Santos et al J 12053 ). 

We confirm this pattern for M87, but in addition, the 
high-quality, spatially extended and spectroscopically- 
confirmcd photometric data set allows us to find more 
detailed evidence for substructures in color and magni- 
tude space. Our next step will be to consider the new 
dimension of kinematics to disentangle the subpopula- 
tions. 

5.4.4. The Connection to Field Stars 

Before continuing with the kinematics analysis, we con- 
clude this section with a look at the results for resolved 
red-giant branch (RGB) stars from deep ACS imaging 
of the cen t ral re gions of M87 (~ 2' or ~ 10 kpc) by 
iBird et aTl (|2010f) . These authors inferred a stellar metal- 
licity distribution which we superimpose on our cen- 
tral GC color distri butions (Figur e |9l le ft panel) after 
using the results of iSinnott et al.1 (|2010) to transform 
(g — i)o color to metallicity. Although such compar- 
isons are ultimately critical for establishing links between 
GC and field star subpopulations, unfortunately it is 
not yet possible to reach strong conclusions. The high- 
metallicity RGB peak that nominally seems to match up 
well with the far-red GC peak is actually truncated above 
(g — i)o ~ 1.12 ([Fe/H] ~ —0.3) because of observational 
limitations. We already know that the field stars in M87 
are r edder overall than the correspo nding red GC peak 
fe.g.. lForte eTall[2007l ISpitlerll2010h . 

The apparent low- and intermediate-metallicity RGB 
peaks may be artifacts of the metallicity modeling func- 
tions as well. Still, it remains a remarkable fact that 
the rich complement of metal-poor GCs is accompanied 
by only an extremely meager metal-poor field star pop- 
ulation (as is already known in other cases such as the 
Milky Way and M31, and discussed in detail by Harris 
l200lf) . The origin of this high "specific frequency" is a 



great mystery, but it does allow metal-poor stellar halos 
in distant systems to be probed indirectly using GCs. 

6. KINEMATICS 

Here we provide various characterizations of the kine- 
matics of the M87 GC system, including its rotation, 
velocity dispersion, and velocity kurtosis. A key point 
throughout is to compare the results using the previ- 
ous versus the new velocity data. Also, unless otherwise 
stated, our "GC" samples used for kinematical and dy- 
namical analysis have been pruned of UCD-candidates 
(any with sizes of r^ > 5 pc or magnitudes io < 20), 
owing to the likelihood that these have systematically 
different orbits from the bulk of the GC system. 

One complication i s that as discussed in 
IRomanowskv et ail (|2011l ). analysis of position- velocity 
phase space with M87 GCs reveals a large unrelaxed 
substructure in the ~ 50-100 kpc region. For the 
moment will we ignore this fact and proceed with 
conventional approaches to kinematical and dynamical 
analysis that assume a well mixed and relaxed system. 
In most galaxies studied to date with GC kinematics, 
the velocity data were either not extensive or high 
resolution enough to detect substructures, while even in 
M87 there may well be additional lurking substructures 
that we have not yet clearly identified. Therefore the 
conventional approaches may have widespread liabilities 
which are beyond the scope o f this paper to consider in 
detail fcf. lYencho et al.ll2006t) . 

We begin in £16.11 with an overview of the individual 
GC positions and velocities, then present our "kinemet- 
ric" techniques along with some initial analyses in £16.21 
We carry out more detailed kinematical analyses of sub- 
populations in £j ^6.3l and 16.41 Looking further forward, 
£18.41 will examine potential halo transitions, with £18.51 
discussing possibilities for a recent merger. 

6.1. General Overview 

A succinct way to illustrate the novel regions of GC 
parameter space in M87 covered by the new data is the 
galactocentric radius versus magnitude (Figure [TU| . The 
radial range is now extended from ~ 8' (~ 40 kpc) pre- 
viously to almost 40' (~ 200 kpc), and the magnitude 
limit deepened from io ~ 21 to io ~ 22.5. The one area 
where the new data do not surpass the previous is a nar- 
row radial range around ~ 6'(~ 30 kpc). It should also 
be noted that the central regions where the old and new 
datasets have the most overlap are completely disjoint in 
the magnitude ranges probed. Therefore dependencies 
of the kinematics on magnitude should be kept in mind 
as possible explanations for any discrepancies found be- 
tween the old and new results. 

The two-dimensional positions of the velocity data 
are illustrated by Figure [TTJ The new LRIS and Hec- 
tospec data provide fairly uniform azimuthal coverage at 
small and large radii, respectively. The median veloc- 
ity of the combined GC data from these two regions is 
1300 km s^ 1 , with the medians from the two data sub- 
samples consistent within the uncertainties . We adopt 
v svs = 1307 km s _1 as the systemic velocity (jSmith et al.l 
120001) . The DEIMOS masks provide additional data 
points in the south and east regions, including an east- 
wards extension to very large radii. As we found above, 
the GCS appears to be significantly flattened, and it may 
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Fig. 10. — Magnitude versus galactocentric radius for spectro- 
scopic objects associated with M87. The io-band magnitude is 
shown on the left axis, and absolute magnitude on the right, with 
a dashed line indicating the approximate GCLF turnover magni- 
tude. Gray X sym bols show the previous data compilation from 
IHanes c t al. (2001), while circles show the new data, with open 
blue and filled red circles indicating "blue" and "red" GCs, respec- 
tively. The new LRIS data are found at small radii between i = 21 
to 22; the DEIMOS data are at large radii down to i ~ 23, and the 
Hectospec data are at a range of radii down to i ~ 21. Filled and 
open triangles show UCDs and transition objects (with gray and 
brown colors for old and new data, respectively). 



be appropriate to consider GC kinematics stratified on 
ellipses. Since the red GCS is consistent with the flatten- 
ing q and position angle of the stellar isophotes, we will 
adopt these better-determined quantities as our default 
for the elliptical circular-equivalent radius R m = R ay /q. 

We next plot the GC velocities versus radius in Fig- 
ure 1121 including both old and new data highlighted by 
different color schemes. A number of interesting features 
are visible by eye, where in particular, the scatter in ve- 
locities around v sys is related to the velocity dispersion. 
In the R ~ 1.5'-5' range, the old and new data behave 
similarly, but at ~ 5'-8', the new data appear to have a 
much lower dispersion. At even larger radii (> 8'), the 
data show strong "clumping" near v sys . The latter fea- 
ture is probably related to intrinsic unrelaxed substruc- 
ture in the outer GCS, w hich we analyze in detail in 
iRomanowskv et all (|2011[ ). In the meantime, we caution 
that our results outside ~ 8' (~ 40 kpc), in the region 
of suspected substructure, should generally be viewed as 
provisional because we will assume a well-mixed steady- 
state system which may be a poor description at large 
radii. 

The discrepancy between the old and new data at in- 
termediate radii (~ 25-40 kpc) is puzzling. A strongly 
increasing dispersion at the radial limit of the old data 
was a key result relied on by ensuing dynamical analyses 
(which we will discuss further in <JT]). We have seen in 
£|4.2l that some of the high relative velocity measurements 
that drove the high dispersion were due to catastrophic 



observational errors, which we suspect may also account 
for some of the remaining high-velocity measurements. 
However, even excluding the most extreme velocities, the 
older dispersion estimate in this region appears higher 
than ours. 

There do not appear to be any obvious differences in 
the GC properties or azimuthal positions to explain this 
discrepancy (which we will examine in more detail in 
^6.2.3|) . but intriguingly, there are four dEs located at 
projected distances of 7'-9' from M87. It is possible 
that one or more of these galaxies is somehow "stirring 
the pot" or else shedding stripped GCs in such a way 
as to cause localized spikes in the velocity dispersion. 
Such features might also be related to the aforementioned 
substructures that begin around these radii. Systematic 
study of this region around M87 would clearly be of in- 
terest. 

There are two low-velocity confirmed GCs (and one 
transition object with ~ 8 pc) around ~ 300- 
400 km s , which as shown by the dotted curves in 
Figure [T^] reside around the 3-tr boundary of the over- 
all GC velocity distribution. For ~ 460 velocities drawn 
from a Gaussian distribution, it is expected that one ob- 
ject on average is found in the 3-<r tail, and more than 
one if the distribution is leptokurtotic, i.e., with extended 
tails as generally predicted for a radially-biased distribu- 
tion of orbits. In addition, there are several high- velocity 
objects just inside the upper 3-er boundary, so we will by 
default consider the low-velocity duo to be part of the 
"normal" GC population bound to M87 (rather than, 
e.g., part of an unbound IGC population). 

6.2. Kinemetry 

We next begin to put the kinematics information into a 
simplified quantifiable format using the general approach 
of "kinemetry" , which extends the standard techniques 
used in ellipse-based galaxy surface phot ometry to quan- 
tify spatial variations of the kinematics (|Krainovic et al.l 
120061 [20081 ). This approach is designed to compress two 
dimensions of kinematics information into one dimen- 
sion by assuming that the relevant quantities (rotation 
in particular) are approximately stratified on ellipses. In 
a circular system, the kinemetry reduces trivially to sine- 
curve fitting, which is a well worn technique in the con- 
text of GC kinematics. However, the strong ellipticity 
of the outer regions of M87 suggests the more general 
approach — even though we will find in this case that the 
kinematics results are not strongly influenced by the el- 
lipticity modification, and we will end up using the cir- 
cular model in order to allow for strong kinematic twists. 

Below we present the general techniques in ^6.2.11 
which we then apply to the stellar kinematics in ^6.2.21 
and to the overall GC system in ^6.2.31 

6.2.1. Methods 

Following the methodology introduced in 
iProctor et al.l (|2009f) for application to sparsely sampled 
data, our basic kinemetric model for fitting a rotation 
field as a function of azimuth 9 in an elliptical annulus 
is: 

^mod(fl) = v sys ± — ]= == =, (2) 
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Fig. 11. — Two-dimensional distrib utio n of GCs with velocity measurements; t he left panel shows a n overall view and the right panel is 
zoomed in. Symbols are as in Figure 1 101 with g r een sq uares showing PNe from Dohcrty ct al. (2009). Dashed blue ellipses show sample 
idealized stellar isophotes from Kormcndy ct al. (2009), at R m ~ 1', 4', 10', 17', 31'. Larger blue circles show five dE galaxies, with the 
sizes marking their estimated tidal radii. 
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Fig. 12. — Radial velocities versus the galactoccntric radius. Symbols are as in Figure [TT1 with light blue star symbols showing objects 
identified as stars. One of the dEs is hard to see at ~ 9', 1350 km s — 1 . Error bars show GC velocity measurement uncertainties in the case 
of the old data; for the new data, these are not shown since they would be in almost all cases smaller than the data points. The orange 
dotted curves show the ± 3 a envelopes for the new GC data set. 
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where vq is the maximum rotation amplitude, 9q is the 
direction of maximum rotation (the kinematical major 
axis JHU q is the "axis ratio" of the velocity field, and the 
+ and — signs are for a position angle inside and outside 
the range (9 — 9 ) — [—tt/2, +tt/2], respectively. 

In the case of integrated-light stellar kinematics, the 
mean velocity is directly measured at each data point, 
(vi), which can be matched to the model v mo d(6) using 
X 2 fitting. A model for the dispersion <7 Pimo d is then 
fitted separately. We also augment the velocity measure- 
ment uncertainties by ~ 5 km s" 1 in quadrature, to keep 
outliers and systematic departures from the simple kine- 
mctric model from unduly influencing the results. 

For discrete velocities, the kinematical model involves 
both rotation and dispersion simultaneously, and uses a 
likelihood function to evaluate the probability of a mea- 
surement (vi ± Aitj) being drawn from a Gaussian distri- 
bution of model velocities. The equivalent x 2 function 
to be minimized is then: 

X 2 = E / l Z T7gL + ^ [< mod + (A,,) 2 ] 0) 

2 p.mod ' \ V 

(cf equation A. 2 oflBergond e t al. 2006|and equation 4 of 
iGnedin et al.ll2010h . The same appr oach is also u s ed for 
kine metry with d i screte velocities in lFoster et al.l (|2011[ ) 
and I Arnold et~aT1 ll20lTI FL 

With integrated light kinematics, a kinemetric model 
may be refined in great detail, including fitting for sepa- 
rate flattenings of both rotation and dispersion. With 
discrete velocities, more restrictive assumptions must 
usually be adopted, as at least ~ 1000 velocities are 
needed to constrain rotational flattening. We try out two 
model alternatives with the GC data: one that is circular 
and allows the position angle to change (<7kin = 1, #o free), 
and another where the kinematical ellipticity follows the 
stellar isophotcs including their variations with radius 
(<7kin = <?phot,0o,kin = #o, phot)- The galactocentric radii 
for the GC measurements are modified accordingly for 
each scheme. 

The uncertainties in the fitted parameters are esti- 
mated by constructing a Monte Carlo series of 3000 mock 
data sets, using the best-fit kinemetric model as a start- 
ing point. Then at each GC location, a random velocity 
is drawn from an underlying Gaussian distribution de- 
fined by the intrinsic dispersion and the spectroscopic 
measurement uncertainties. We fit a kinemetric model 
to each mock data set and use the 68% range of result- 
ing model fits to define the 1 a uncertainties. The error 
bars estimated in this way are probably somewhat un- 
derestimated because the azimuthal sampling effects are 
difficult to fully take into account. 

We also use the Monte Carlo simulations to correct for 
bias in the kinemetric parameters. In particular, when 
9 is a free parameter, v rot will be biased high. E.g., a 
non-rotating system that is sampled with a finite num- 
ber of velocity measurements Ni will always appear to 

26 This convention for 6>o differs from Romanowskv et al. (2009), 
where we used the direction of the angular momentum vector. 

27 In the case of equal uncertainties Avi, equation \3\ reduces 
to standard least-squares fitting. On the other hand, standard x 2 
fitting as used in some previous studies is mathematically incorrect 
and can yield very wrong answers when Avi -C a p and when v TO t <C 

Up. 



have non-zero rotation if the direction of rotation is freely 
chosen. In general, the bias becomes important when 
the rotation is relatively low, with the rule of thumb 
that an observed free- P.A. ro tation estimate is insecure 
if (wrot/ffp) £ 0.55 x ^/20/N t . 

As far as we know, this effect has not been explicitly 
considered in previous studies of rotation in GC and PN 
systems. 

For the simple dynamical models that we will con- 
struct in this paper, we need only an azimuthally aver- 
aged second-order velocity moment, v rms (R). We could 
use the maximum likelihood techniques as in the kineme- 
try formalism, assuming Gaussian velocity distributions 
and forcing the rotation to be zero. However, to avoid 
the Gaussian assumption as much as possible, we instead 
use a classical velocity dispersion estimator: 

^ i=JV 

"L = ]y - v *ys) 2 ~ ( Aw 2 - (4) 

The uncertainty in this quantity Av rms is 
then estimated through f ormu lae provided by 
iDanese. de Zotti. fc di Tulliol (|1980tl which do have 
to assume Gaussian velocity distributions. 

It is also important to estimate the deviations of the 
velocity distributions from Gaussianity. To this end, we 
adopt the classical fourth-order moment, the projected 
velocity kurtosis: 




where v rms is calculated before correction for the observa- 
tional uncertainties Avi, and in practice we modify the 
simple express ion above for k„ bas ed on bias-corrected 
formulae from lJoanes fc GTH ([19981 ). 

6.2.2. Stellar Kinematics 

Before continuing to the GC kinemetry, we present 
results on the stellar kinematics, where the kinemetry 
is conceptually more straightforward. The two main 
datasets we will consider are from tw o-dimensional in- 
tegra l field spectrographs: SAURON (Cappcllar i et al.l 
120071) out to ~ 0.6' (~ 3 kpc), and VIRUS-P 
(|Murphv et"al~ll201lD out to - 4' (~ 19 kpc). Note that 
the SAURON data are a revised version using an im- 
proved stellar template library, which does make a sig- 
nificant difference in the case of M87. 

The kinemetric results in three sample bins in radius 
are shown in Figure 1131 In the innermost regions (inside 
- 0.6' or - 3 kpc), both the SAURON and VIRUS-P 
datasets imply very low rotation, at the ~ 5 km s _1 
level. Differences between these datasets in the direction 
of the rotation may not be significant owing to incom- 
plete azimuthal sampling. At larger radii (~ l'-4' or 
~ 5-19 kpc), the VIRUS-P data suggest the rotation in- 
creases but to merely ~ 20 km s _1 . 

We next show a summary of kinemetric parameters 
plotted versus galactocentric radius in Figure 1141 Some 
of the features of this plot will be discussed later, but 
first we highlight the rotational behavior. It appears 
that the rotation is aligned with the major axis inside 
~ 0.3' (~ 1.5 kpc), and then twists fairly sharply by 
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Fig. 13. — Mean velocity versus position angle in three radial 
bins (as labeled in the panels). The purple stars and blue circles 
are SAURON and VIRUS-P data, respectively. Kinemetric fits to 
the data are shown as solid curves with dashed curves showing the 
uncertainties on the rotation amplitudes. 

90° to alignment with the minor axis. Such a kinemati- 
cal twist does not necessarily imply a radical dynamical 
transition but may reflect the rearrangement of a small 
fraction of the stellar orbits, as found through dynamical 
modeling of t he famous "kinematically d ecoupled core" 
in NGC 4365 (|van den Bosch et al.ll2008h . 

The next thing we notice is that the SAURON and 
VIRUS-P data differ in their velocity dispersion and hi 
(fourth-order Gauss-Hermite moment) profiles, in some 
areas at very high confidence levels given the appar- 
ent error bars on the data. This is a demonstration of 
the sometimes severe system atic uncertaint i es in a nalyz- 
ing stellar kinematics data. iMurphv et al.l (|2011l) et al. 
found for the dispersion results at least that the differ- 
ence could be attributed to template mismatch involving 
the Mg b region, which is a more dominant component 
in the SAURON data. 

6.2.3. Globular Cluster Kinematics: Old Versus New 

Moving on to kinemetry of the M87 GCs, we begin with 
the simple circular model, combining all of the data irre- 
spective of color (i.e., blues and reds together) but sep- 
arating the analysis into the "old" and "new" measure- 
ments, where the year 2003 is used as the boundary of 
obsolescence. For the purposes of this comparison, we do 
not remove old velocity measurements with known catas- 
trophic errors (e.g., S878, S1074), but we do reclassify 
the objects in some cases based on new multi-color pho- 
tometry and HST size measurements. Some "GCs" are 
now identified as foreground stars or background galaxies 
(e.g., S1472, S7008, S7012), while some "contaminants" 
are reclassified as GCs (e.g., S5002, S7017). This is dis- 
cussed in i j4.4l We also remove UCD candidates and 
transition objects (r n > 5.25 pc; ^4.4.1 1) . as well as bright 
objects (io < 20) because of their potentially disparate 
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Fig. 14. — Kinematics of M87 versus galactocentric radius. The 
panels from top to bottom show rotation position angle and ampli- 
tude, velocity dispersion, and fourth-order Gauss-Hermite moment 
/14. Different data sources are indicated by symbols as given in the 
legend of the second panel. The kurtoses from the discrete velocity 
data (PNe and RGCs) are transformed to I14 by the approxima- 
tion k p ~ 8V6h,4 Ivan der Marel fc Frarmlll993f) . The dashed blue 
curv es in the top p anel show the V-band stellar photometric major- 
axis (Kormcndy ct al. 20(H). 

kinematics (to be examined in ^6.4|) . These procedures 
allow us to focus our comparison of old and new data 
sets on the spectroscopic aspects. 

Before showing our results, w e will descr i be bu t not 
show explicitly the results from iCote et ail (|2001l also 
using a circular model). They found a rotation of 
Vrot ~ 150 km s" 1 along the minor axis inside 3.5' 
(~ 18 kpc), particularly for the blue GCs. At larger 
distances they found major-axis rotation rising with ra- 
dius, up to u ro t ~ 400 km s -1 at ~ 7' (~ 35 kpc). They 
also found velocity dispersions of cr p ~ 350 km s _1 in the 
inner regions, rising suddenly outside ~ 6' (~ 27 kpc) to 
cr p ~ 500 km s _1 . The P.A. twist and high outer rota- 
tion and dispersion were also noted from earlier data by 
iCohen k Rvzhovl (|1997t ) and iKissler-Patig k Gebhardl 
(|1998D . 

Now binning the data by galactocentric radius, and 
adopting a fixed v sys = 1307 km s , we show our kine- 
metric outcomes in Figure 1151 Using the old data, we 
generally reproduce the previo us results except t hat we 
find even stronger rotation than lCote et al.l (|2001[) did at 
the outermost radii (v I0 t ~ 500 km s _1 ). These differ- 
ences are presumably due to having a "purer" GC cata- 
log. 

The results from our new data show some further re- 
markable differences, as already suggested in §6.11 In- 
side ~ 5' (~ 25 kpc) we do again see signs of kinematical 
twisting, but the rotation amplitude is lower, and the 
twist orientations are different from in the old data. We 
suspect these inner P.A. fluctuations are not significant, 
since with v Iot /cr p ~ 0.15, it is difficult to determine the 
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Fig. 15. — Kinemetric fits to M87 GC data, versus galactocen- 
tric radius. The fitted parameters are the direction of maximum 
rotation (which wraps around past 360°), the rotation amplitude, 
and the velocity dispersion (panels from top to bottom). The filled 
orange circles show our new data, and the open gray circles show 
the older literature data. Dashed blue curves in the top panel 
show the photometric major axis of the M87 stellar light from 
IKormendv et al.l (2009, with ±180° ambiguity). 

P.A. robustly. The new rotation is in most locations 
consistent with major axis rotation, although there is 
a hint of twisting to minor axis rotation outside ~ 11' 
(~ 55 kpc). 

In general, the rotation amplitude is lower than in the 
old data, with v Iot ~ 50 km s _1 typically. In particular, 
the drastic outer rotational increase has vanished. The 
old and new dispersion profiles are generally consistent, 
except outside ~ 5' (~ 25 kpc): the new tr p declines 
where the old data showed a rise (as discussed above). 

In order to understand the reasons for these differences, 
we consider the results at ~ 8' (~ 40 kpc) in more de- 
tail. In Figure [16] we plot the individual GC velocities 
versus position angle for the old and new data in the 
same radial bin. These two data sets share the same 
systemic velocity and general direction of rotation, but 
otherwise look as though they were drawn from a differ- 
ent galaxy. The clear high-rotation pattern in the old 
data is strongly inconsistent with the new data, and the 
new velocity dispersion is significantly lower. Some of 
the highest old velocities in this bin have now been iden- 
tified as catastrophic measurement errors, and it is plau- 
sible that some of the other high velocities are caused by 
similar errors. However, the differences are not simply 
due to this effect, as the velocity distribution at ~ 180° 
is peaked near 1700 km s _1 (i.e., 400 km s _1 higher than 
v sys ) and some low velocities near 0° are found as would 
be expected from real rotation. 

Both data sets do show strong fluctuations in the kine- 
matics around this radial range, which may be part of 
the story, if there are very localized substructures with 
distinct kinematics whose detection depends on having 



; 


= 6.8-10.5' 

i? A T 
S 


o 


s. 


1 

• 


old : 


- 








o 

° 


• i J 


: • 


= 585 ±126 km s"' , <7 p = 474 ± 61 k 


- s 
m s" 


i i~ 


r 

J - 


V 

: r p 


1 1 1 1 1 1 1 1 1 — 

= 6.8-10.4' 


— I — 


H 1 — 




H 1 

new 




__«---- 




A 








o o • • o 


_D 




•• 


a 

o - 




• o o O O 
o o 

•o o 






oo 




D — 




= 80 ± 54 km s" 1 , a p = 307 ± 28 km 


S" 1 









100 200 



P.A. [deg] 

Fig. 16. — Radial velocities versus position angle in a radial bin at 
R p ~ 8' (~ 40 kpc), for old (top) and new (bottom) data. Symbol 
types are is in the previous two figures. The kinemetric solutions for 
these bins are shown as solid orange sinusoidal curves with dashed 
curves outlining the ± 2 a envelope; the solutions' rotation and 
dispersion parameters are listed in the panels. The PNe, UCDs 
(including the peculiar object S7023), and transition objects are 
not used in the fits. The black dot-dashed horizontal lines mark 
the systemic velocity. 

complete or fortunate target selection. In particular, the 
new data set has sparse coverage at R ~ 6.7' (~ 32 kpc) 
where many of the old extreme velocities were, and we 
highlight the need to provide more coverage in this area 
and to repeat more of the old velocity measurements. 
Another factor that will be discussed in §6.41 is that the 
brighter GCs targeted in the old surveys are systemat- 
ically different from the fainter ones that we can now 
reach. 

At present we are unable to conclusively assess the 
reliability of the old data. The differences between the 
old and new data may be due both to sample selection 
and to problems with the old data. Since a small number 
of large outliers can significantly skew the results, for the 
rest of the paper we will rely only on the new data set 
unless otherwise stated. 

6.3. Subpopulations: Parameter Surveys 

We next consider the kinematics of various subpopu- 
lations, starting with trends with color in c; q6.3.lH6.3.3l 
and with luminosity in £16.3.41 

6.3.1. Overview 

This extensive data set, with high quality photometry 
and velocities, allows us to consider the generic question 
of bimodality ("blue" and "red" GCs) in a new way: we 
investigate whether or not the GC kinematics show a 
strong transition or discontinuity with color that would 
imply two disconnected subpopulations. 

We could attempt to define subpopulations based on 
the six potentially-relevant parameters in our data set 
(color, magnitude, size, velocity, 2-D position). However, 
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Fig. 17. — Smoothed rms velocity fields of GCs in the plane of color versus log of the galactocentric distance. White points show the 
individual velocity measurement locations, and the color bars at the bottom indicate the v TulB scale (in km s — 1 ). On the left is a mock 
data set, and on the right is the true data set of new M87 GCs. A joint color-kinematical bimodality was input in the mock data, and 
correctly recovered as seen by the strong v TulB transition at small radii. The transition seen in the real data is more complicated, suggesting 
multiple subpopulations. 



such an exercise is beyond the scope of this paper (cf 
IChattopadhvav et~a l. 2009), and even our extensive data 
set will provide only sparse coverage of a six-dimensional 
parameter space, so for now we will carry out a more 
cursory treatment. 

We begin with a non-parametric study in the space of 
color, (g — i)o, versus galactocentric distance, R p , since 
both of these parameters are likely to be important to 
the GC kinematics. The simplest kinematical parameter 
that does not require detailed modeling is the rms veloc- 
ity Vrma- We will therefore examine the dependence of 
v Tms on color and radius jointly. 

6.3.2. Simulations 

Before analyzing the real data, we generate some mock 
data sets in order to optimize our techniques to pick up 
only real features in the data. We take the same positions 
R p for the mock data as in our real data set, and ran- 
domly generate a color and velocity for each object. The 
color is drawn from either a blue or red subpopulation, 
with its probability of belonging to each subpopulation 
being a function of radius motivated by our surface den- 
sity profile analysis f £|5.ip . This probability ranges from 
45% blue at R p ~ T to 85% blue at - 40'. For simplicity, 
we choose constant rms velocity profiles with radius, with 
v lms = 400 km s _1 for the blue GCs, and 250 km s _1 for 
the red GCs. The (g — i)o color distributions are assumed 
to be Gaussian with peaks and 1-ct widths of (0.79, 0.06) 
and (1.05,0.09) for the blue and red GCs, respectively. 

After generating the data set, we construct a grid of 
color versus log-radius, and at each grid point derive a 
smoothed u rms value based on a weighted average of the 
surrounding data points. The weighting uses a Gaus- 
sian kernel with "distances" of 0.22 dex in log-radius and 
0.07 dex in color. 

The results for one mock data set are shown in the left 



panel of Figure 1171 An obvious feature in the central re- 
gions is the strong transition in u rms correctly recovered 
at (g — i)o ~ 0.95, reflecting the presence of two dis- 
tinct underlying subpopulations. At larger radii where 
there are few red GCs, the color-location of the kine- 
matical transition drifts away from the input value, pre- 
sumably because of redward contamination for the blue 
subpopulation. There are also localized fluctuations in 
v Tms caused by statistical fluctuations — particularly in 
regions with very few nearby data points. 

6.3.3. Kinematical Evidence for Subpopulations 

We next consider the real data in the right panel of Fig- 
ure El where again, we are now using only the new data 
set, and are omitting the known "UCDs" (rh > 5.25) as 
well as the bright objects (io < 20); there are 410 objects 
in this revised sample. Some "hot" and "cold" spots can 
be seen which are in regions with few data constraints 
and so can be dismissed as noise. As in the simulation, a 
kinematical "boundary" is indeed seen at (g — i)o ~ 0.95. 
However, there are a number of dissimilarities with the 
simple model that appear subtle but may be important. 

One difference is that the kinematical transition with 
color is less gradual in the data than expected, with v Tms 
reaches its maximum value at very blue colors of (g — 
*)o J$ 0.75. Another is that the transition color shifts 
blueward at larger radii, i? p > 8' (> 40 kpc), while is 
seen in the simulation, we would expect a redward shift 
if anything. 

To investigate these features in more detail, we make 
standard one-dimensional profiles of v Tms versus color in 
select radial bins, and versus radius in select color bins. 
Figure US] shows the results. We first consider the radial 
bin R p = 3'-8' (15-40 kpc), which we may consider the 
main body of the GC system. We see that there may be 
a discontinuity at (g — i)o ~ 0.83, with GCs bluewards 
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and redwards having v rms — 370 ± 38 km s" 1 and 255 ± 
23 km s _1 respectively. 

This result would seem initially to not confirm the sim- 
ple bimodality model with a color boundary at (g — i)o ~ 
0.93, and may even be c onsistent with a non-lin ear color- 
metallicity relation (e.g. JYoon. Yi. fc Ledl200"6l ). since an 
assumed smooth u rms -metallicity anti-correlation might 
then lead to a strong increase of observed ?; rms at very 
blue colors. However, such a scenario would need addi- 
tional data and modeling to test, and in any case we have 
already seen that the M87 GC color distribution is more 
complicated than in a classical bimodality picture. We 
found evidence in i j5.4l for at least a third GC subpopu- 
lation with intermediate colors of (g — i)o ~ 0.86-1.01. 
It appears that these objects have relatively "cold" kine- 
matics, which is similar to the redder GCs, and compli- 
cates the kinematics tests for distinct subpopulations. 

In Figure [TH] we also plot v rms versus color in other 
radial bins, and versus radius in different color bins, and 
see that the picture becomes even more complicated. For 
example, at large radii, both the very bluest and the very 
reddest GCs have higher w rms than the bulk of the GCs 
at the same radii (see also Figure 117]) . These kinematical 
irregularities seem to be driven largely by a small number 
of extreme- velocity object^, and could be due to recent, 
discrete accretion events resulting in families of GCs that 
are clumped in chemo-dynamical phase-space. 

We therefore conclude that the outer regions of M87 
may host a significant fraction of unrelaxed material that 
makes kinematical analysis of subpopulations in these re- 
gions somewhat precarious. The outer subpopulations 
may be considered still in the process of assembly, and 
might not even bear much connection to the central sub- 
populations. These conclusio ns, based on basic color an d 
v lms trends, parallel those of IRomanowskv et al.l ([201 1) , 
who analyzed phase-space in more detail. 

A similar explanation may apply to a kinematical pecu- 
liarity of the reddest GCs at small radii (R p = 1.6'-2.2' 
or 8-10 kpc). These objects have a much higher v rms 
than at intermediate radii, which we trace to a peculiar 
offset in their overall velocities. As shown in Figure Q1J1 
the red GCs inside 2.2' (10 kpc) have a median velocity 
of 1464 ± 49 km s _1 , which is significantly offset from 
the overall v sys — 1307 km s , and is driven mostly by 
the far-red objects with (g ~ i)o <; 1.01. This effect does 
not seem to be driven by uneven sampling of a strongly 
rotating system, since the offset objects are distributed 
over a wide range of azimuths, and we cannot think of a 
plausible observational or instrumental effect to explain 
it. 

Inspecting the older data set, the offset also appears in 
the red GCs, and can be seen to extend to even smaller 
radii (Figure As mentioned in £ )5.4[ there are also 
indications of peculiar shifts in the red GC colors at pre- 

28 Note that among the far-red GCs, H20573 and H66419 have 
similar colors, luminosities, and phase-space locations, supporting 
a common origin. For the far-blues, the extreme-velocity objects 
are T16997, H23346, and H58443, which also lie in a similar re- 
gion of the color-color diagram (see earlier discussion of the latter 
object). None of these objects has a measured size but we do not 
consider them as likely stars (in the case of those with low ve- 
locities), but instead as probably parts of unrelaxed substructure 
aroun d M 87 (see also two PNe close to H58443 in phase-space in 
Figure □!}. 



cisely these radii, supporting a picture of the inner re- 
gions of M87 as incompletely mixed. 

6.3.4. Trends with Luminosity 

We next consider magnitude dependencies of the GC 
kinematics. Because of the strong variations in the mag- 
nitude selection with radius in our data (Figure [TU|) , it 
is difficult to consider this theme in much detail. We 
adopt the simple approach of examining the velocity dis- 
tribution as a function of radius for "bright" and "faint" 
subsamples, separately for blue and red GCs. 

For the red GCs we do not find any strong trends with 
magnitude, but for the blue GCs we find a possible tran- 
sition at in ~ 20 (Mi ~ — 11) which we illustrate in Fig- 
ure [201 (left panel). For the innermost regions where we 
have kinematics data on bright blue GCs, these objects 
appear to avoid the systemic velocity. Although we have 
very low number statistics for these objects, intriguingly, 
the bright UCDs (of all colors) appear to share in the 
same pattern. 

A key question here is whether these kinematical pe- 
culiarities are primarily linked to luminosity or to size, 
since four out of the five bright blue "GCs" have un- 
measured sizes and might be UCDcH To investigate 
further, we break down the sample in the radial range 
R p ~ 2.5'— 7' (~ 10-30 kpc) into four subsamples of size 
and magnitude, and show their velocity histograms in the 
right panel of Figure [201 Because some of the subsamples 
include very few objects, we have boosted the statistics 
by calling objects with intermediate sizes (rh ~ 5-10 pc) 
"UCDs" , and by including all of the old data in the hope 
that the occasional rogue velocity measurement will not 
muddy the waters. 

For each of the four subsamples, the velocity dispersion 
(from a Gaussian fit) and kurtosis are indicated in the 
Figure panels. We see that the bright and faint compact 
blue GCs have consistent dispersions and kurtoses, sug- 
gesting that these belong to the same population. Both 
UCD subsamples, on the other hand, have higher dis- 
persions than the GCs, which provides another line of 
evidence that UCDs (as defined by large sizes but not 
necessarily high luminosities) are a distinct population 
from normal GCs. 

This picture becomes more complicated when compar- 
ing the UCD subsamples, which although having simi- 
lar velocity dispersions, have discrepant kurtoses. The 
bright UCDs have a negative kurtosis, driving the pe- 
culiar "double-peaked" velocity distribution previously 
mentioned, while the faint UCDs have a higher kurtosis. 
There are a number of potential explanations for this dif- 
ference: multiple populations of UCDs; additional cases 
of catastrophic measurement errors; or a statistical fluc- 
tuation (the kurtoses are different at only 1.5 a). Firmer 
conclusions about UCD kinematics will require an en- 
larged set of new data. 

We will return to additional luminosity trends in £16.4.41 

29 These four GCs are S279, S348, S501, and VUCD10. The 
fifth, S1265, has a normal GC size, and is the closest of the five 
objects to v ByB , but intriguingly has a very similar distance, color, 
magnitude, and velocity to the bona fide UCD S1629. Also, two 
bright transition objects that fit the same pattern are S77 and 
S13 7. which have sim ilar positions, velocities (from Cohen 2000 
and Hancs et al. 2001), colors, magnitudes, and sizes to each other. 
These seem likely to share a common origin. 
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Fig. 18. — Root-mean-square velocity profiles of GC subsamples. The left panel shows i) rms versus GC color, in three radial bins (see 
legend in panel). The right panel shows t> rm s versus galactocentric distance, in three color bins (see legend). 
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Fig. 19. — Velocities versus galactocentric distance for redder 
GCs. Old data with (g — i)o > 1.01 are shown with gray X symbols, 
while the new data are shown as orange open pentagons for 0.86 < 
(<? — *)o < 1.01, and red filled circles for (g — i)o > 1.01. At small 
radii, the reddest GCs show a strong velocity asymmetry relative 
to v aya (the UCD candidates, which arc not plotted here, have the 
same behavior). 



and discuss some general implications in H8.3I To be safe, 
we have in general simply omitted all objects with ig < 
20 (whether blue or red) from our kinematical analyses 
in this paper. 



6.4. Subpopulations: Kinemetry 

We now move on to kinemetric models of the GC sub- 
populations (see methods in ^6.2. H and overall analysis in 
M6.2.3|) . Given our discussion of multiple GC subpopula- 
tions in preceding sections, we adopt a trimodal analysis 
as our new default. The three subpopulations are the 
blue GCs or "BGCs" (analyz ed in §6.4. the interme- 
diate col or GC s or "MGCs" (SgXU), and the red GCs or 
"RGCs"(|nX31). The color b oundaries are (g-i) = 0.86 
and 1.01 (derived in §5.4.2|) . We also exclude the large 
(r n > 5.25 pc) objects (or "UCDs") as well as all bright 
objects (i n < 20) because of the strong s uspicio n that 
they represent different subpopulations (§6.4.41) . Our 
kinemetry results will be summarized in £16.4.51 

An overview of the kinemetric results is shown in Fig- 
ure [2TJ For each subpopulation, radial profiles are shown 
for rotation position angle and amplitude (#0 and u ro t)i 
projected velocity dispersion (er p ), and projected veloc- 
ity kurtosis (k p ). For clarity, the radial binning is the 
same for all kinematic parameters of a given subpopu- 
lation, and so represents a compromise between optimal 
binning for, e.g., dispersion and kurtosis. 

We continue to use circular kinemetric models (q = 1) 
for the results and figures discussed here, in order to allow 
for kinematic position angle variations (see §6.2.1). We 
have also carried out parallel kinemetric analyses using 
flattened models (based on the photometric distributions 
of the stars or the GCs), but, at least for the present data, 
these models turn out to yield similar results. 

Although we have used Monte Carlo methods to esti- 
mate the uncertainties in our kinemetry fits, we suspect 
that these uncertainties are underestimated, particularly 
for cases where azimuthal coverage is sparse. To help 
evaluate the results in Figure [2U we therefore provide 
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Fig. 20. — Demonstration of magnitude dependence of GC kinematics, using io = 20 as the boundary between "bright" and "faint" 
objects, and rh=5.25 pc as the boundary between compact and extended objects. Left panel: Absolute value of the difference between 
object velocity and systemic velocity, as a function of galactocentric radius. Blue GCs and UCDs (of all colors) are shown, in bright and 
faint subsamples (top and bottom). As in previous plots, colored and gray symbols show new and old data, respectively. Symbol types 
denote sizes, with filled and open triangles the UCDs and transition objects, filled squares and circles the "normal" compact GCs, and 
crosses and open circles for unknown sizes. An approximate circular velocity profile is shown in the top panel (see iJTJ. Right panel: Velocity 
distributions, for data in the radial range R p = 2.5'— 7' (in histograms with bins of 100 km s" 1 ), and best-fit Gaussians (solid curves with 

dotted outlines of VN scatter). The top and bottom panels show the bright and faint objects, respectively. The left and right panels show 
"GCs" and "UCDs", distinguished by being compact or extended, respectively. The best-fit values for the Gaussian dispersions and the 
kurtoses are indicated in each panel, along with the 1 a uncertainties. 



Figure l22l to show some details of the fits. We also plot 
some velocity distributions in Figure [23] as a counterpart 
to the kurtosis metrics, and in Table ITU we report some 
numerical values for the kinematics in bins of sub-popu- 
lation and radius. 

6.4.1. Blue GCs 

Considering first the BGCs, we find some indications 
that they do not obey a smooth, simple rotation field. In 
the central regions (inside ~ 6.6' or ~ 30 kpc), there ap- 
pears to be significant rotation (v TO t/<r p ~ 0.3), but with 
a nearly 180° twist near ~ 2.3' (~ 11 kpc) from minor- 
axis rotation to the east to minor-axis rotation to the 
west (in this region, the VIRUS-P data suggest the stars 
have minor-axis rotation toward the west). However, the 
azimuthal coverage along the direction of the maximum 
rotation is poor, and we cannot yet be confident in this 
result. 

There are additional suggestions of BGC kinematical 
twists at larger radii. At two locations (10'-13' or ~ 50- 
65 kpc, and 23'-30' or - 110-140 kpc), high rotation 
along the minor axis to the east is inferred (v TO t /&p ~ 0.5- 
0.9). However, the azimuthal coverage and number of ve- 
locities measured in these bins are rather sparse, and the 
rotation does not seem to persist at radii in between, so 
we will regard this outer-rotation finding as provisional. 
Considering all of the BGCs inside 27' together, they 
have v r ot/0p < 0.22. 

The velocity dispersion and rms velocity profiles of the 
BGCs decline mildly with radius, with an overall power- 
law exponent for the latter of ~ —0.12 ± 0.05. However, 



rather than declining smoothly, these profiles may have 
a sharp transition at ~ 10' (~ 50 kpc), with constant 
profiles inside and outside this radius (at ~ 370 km s _1 
and ~ 300 km s _1 , respectively; see also top panels of 
Figure l22|) . This feature coincides with the possible onset 
of a kinematical twist previously mentioned. 

The kurtosis of the BGCs may also change at ~ 10', 
from k p = —0.1 ± 0.4 in the inner regions to 0.6 ± 0.5 in 
the outer regions (a 1 a difference). We will discuss this 
theme in more detail later, but the implication would be 
for a shift from isotropic to more radial orbits. As an al- 
ternative check on this result, we show the reconstructed 
velocity distributions of these two GC subsamples in the 
top panels of Figure [23l with best-fit Gaussians shown 
for comparison. Here it can be seen that the two veloc- 
ity distributions show similar non-Gaussian deviations 
(peaky profiles in the center, and excess objects at high 
relative velocities), even though their kurtoses are for- 
mally inconsistent. These deviations are suggestive of 
radially-biased orbits in both bins and are a reminder 
that the kurtosis metric can be a blunt tool. 

6.4.2. Intermediate- color GCs 

Turning next to the MGCs, we found earlier that their 
v Tms values were closer to those of the RGCs than those 
of the BGCs. Figure [21] suggests that the MGCs have 
a different rotation profile from the RGCs, particularly 
at large radii, where the two subpopulations may be 
counter-rotating. This would support the idea (already 
suggested by photometry) that there are more than the 
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Fig. 21. — Summary plot of kinematical profiles for M87 subpopulations, as functions of galactocentric radius. The first column shows 
blue GCs; the second shows intermediate-color GCs; the third shows red GCs and PNe (red and green color-schemes, respectively); and 
the fourth shows bright and large objects ("UCDs" with > 5.25 pc). The second and third columns also include the VIRUS-P stellar 
kinematics data. The top three rows show the kinemetric fit parameters Oq, v TO t, and <r p (position angle wrapped around past 360°, rotation 
amplitude, and projected velocity dispersion). The bottom row shows the velocity kurtosis. Points with error bars represent best fits in 
independ ent radial bins, along with the 1 a uncertainties. The dashed light blue curves in the top row show the photometric major axis as 
in Figure [15] 



two standard GC subpopulations. 

The MGCs at small radii (R p < 2.1 arcmin or 10 kpc) 
appear to be kincmatically very cold (cr p = 168 ± 
31 km s -1 ; see Figures \W\ and |2"2"|) . However, with only 
12 objects in this bin, this result might be a chance sam- 
pling fluctuation. 

The MGC kurtosis is generally consistent with zero. 
Figure [23] shows that ther is also an asymmetric excess 
of lower GC velocities at large radii, which may be part o f 
a substructure in the MGCs (jRomanowskv et al~ll2011[ ). 

6.4.3. Red GCs and PNe 



The RGCs also have generally weak rotation close to 
zero at all radii, except for a rise suggested outside ~ 10 ' 
(~ 50 kpc). Figures |2"T1 al s o sho ws results from the 
PN sample of IDohertv et al.l (120091 ) . adopting R p = 40' 
(190 kpc) as a cut-off radius, beyond which there may 
be severe contamination from intracluster PNe (see Fig- 
ure 112]). 

There is remarkable agreement between the PNe and 
the outer RGCs in all four kinematical parameters, which 
is not shared by the other subpopulations, or even by the 
red GCs when using the standard bimodal color division. 
Unfortunately, this conclusion is not certain because as 
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Fig. 22. — Velocities versus position angle in bins of radius and subpopulation (cf Figure lT6)l , The subpopulations (blue GCs, intermediate- 
color GCs, red GCs, UCDs, and bright GCs) are as labeled in each row. The columns show inner, intermediate, and outer radial bins, as 
labeled in the panels. The best-fit rotation position angle and amplitude values are reported in each panel. There are additionally PNe 
included along with the red GCs at large radii, and the kinemetric model there is a joint fit to both subpopulations. The lower panels plot 
the UCDs and bright GCs together where possible, but the kinemetric models shown are generally fits to the UCDs, except for the middle 
panel where a joint fit is applied. Vertical lines mark the semi-major axes of the stellar isophotes. 



shown in Figure [22l there are very few RGC data points 
at large radii, and the PN azimuthal coverage is very 
uneven (affecting the rotation conclusions in particular) . 

Also, it is possible that at these large radii, the very 
low velocities are not from objects bound to M87 but 
from free-floating intergal actic populations tha t are near 
M87 only in projection. iDohertv et all (|2009() adopted 
this interpretation for their three low-velocity PNe, but 
as discussed when introducing Figure I12[ we consider it 
also plausible that these are bona fide M87 objects on 
v ery radial orbit s. 

IDohertv et al.l ((2009) commented that their inferred 
halo rotation from the PNe wa s low, sugg e sting that 
the fast rotation of the GCs from I Cote et al.l (|2001h was 
driven by contaminants or else meant that the GCs do 
not trace the stars. We have found that there were ap- 
parent issues with the older data that produced the very 



high rotation signal, but also that some rotation is still 
possible in both the PNe and the RGCs, depending on 
how outliers are handled. If we did exclude the 5 lowest 
PN and RGC velocities, we would find an dramatically 
lower overall dispersion for these objects at large radius 
(~ 200 km s _1 rather than ~ 400 km s _1 ). One impor- 
tant avenue for future investigation would be to increase 
the number of PNe and RGCs with measured velocities 
at large radii. 

The velocity distributions of the outer RGCs and PNe 
do support the kurtosis agreement, with peaky profiles 
and an excess of very low velocities (Figure |2"51). Such be- 
havior suggests fairly radial orbits (again, assuming that 
these are not intergalactic contaminants). Note that the 
velocity distribution of the inner RGCs is asymmetric, 
reflecting the velocity offset we previously identified for 
these objects (see Figures [T^l and |2"2"|) . 
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Fig. 23. — Velocity distributions of GCs in bins of subpopulation 
and galactocentric radius. The solid curves show the observed re- 
sults, which use optimal Gaussian-kernel smoothing. The dashed 
curves show the best-fitting Gaussians, convolved with the smooth- 
ing kernel. Visible differences between these curves are generally 
significant. The rows from top to bottom show blue, intermediate 
color, and red GCs, and UCDs. The left column shows the inner 
radii, and the right column shows the outer radii, with the bound- 
ary radius Hp = 10' in all cases except 3.5' for the red GCs. The 
sixth panel includes PNe at the same radii for reference, and the 
lowest panels show the bright GCs. Optimal Gaussian smooth- 
ing has been used for the velocity distributions of the data. The 
normalizations are arbitrary. 

There is some overlap bet ween the GCs and th e stellar 
kinematics from VIRUS-P (|Murphv et all 120111) . Both 
the MGC and RGC subsamples show some kinematic 
similarities to the stars (see also Figure [14]), although the 
stars do not show the large mean velocity offset found in 
the RGCs in particular. 

6.4.4. Luminous GCs and UCDs 

We also attempt to examine the kinemetry of the 
bright GCs and the UCDs, which is difficult because of 
the fairly small data sets for these objects. First of all, 
Figure I2T1 suggests that the bright GCs and UCDs have 



consistent kinematics in general, except for remarkably 
high rotation at large radii for the UCDs, which is not 
shared by the bright GCs. However, as shown by Fig- 
ure El there is one low- velocity UCD (H38554) that is 
driving half of the rotation, and without this object the 
two subpopulations would be more consistent. The veloc- 
ity distributions further show some similar asymmetries 
(Figure l23l see also Figure [20)1 . 

More generally surveying the different color subpopu- 
lations' kinemetry, the UCDs seem to be generally con- 
sistent with the MGCs as well. Because of the lower GC 
velocity dispersion th at we have now foun d, we do not re- 
produce the finding of lFirth et al.l (|2008l ) that these have 
a higher dispersion than the UCDs around M87. 

One interesting result here is that not only do the 
bright GCs have a higher velocity dispersion overall than 
the faint ones ( £16.3]) . but their rotation and dispersion 
profiles both show spikes (of marginal significance) to 
relatively high values at ~ 8' (Figure l2"Tj) . This is the 
same radial region where older work found elevated ro- 
tation and dispersion values (see £16.2.31) , suggesting that 
those results were driven by behavior in the bright GCs 
that does not reflect the bulk of the GC system, and may 
trace primarily UCD kinematics (see also £J6.3.4p . 

6.4.5. Ktnemetry Summary 

To summarize the kinematic results, overall there is lit- 
tle dynamically significant rotation found in M87, with 
some localized suggestions of more substantial rotation. 
The rotational "blips" may be due to unmixed substruc- 
tures, and in general will require more complete spectro- 
scopic coverage to ascertain definitively. We see various 
differences between the BGCs, MGCs, and RGCs that 
tentatively support the distinct nature of these subpop- 
ulations. The UCDs and bright GCs also appear to have 
distinct kinematics which are most closely related to the 
MGCs. The central stellar kinematics has some similar- 
ities to both the MGCs and RGCs, while the PN kine- 
matics is most similar to the RGCs. See also Table [14] 
for a summary. 

7. DYNAMICS 

M87 has been one of the most intensively modeled 
galaxies in the Universe, and it is beyond the scope of 
this paper to carry out another detailed dynamical anal- 
ysis or even to adequately review the previous work. This 
will be the subject of future papers, and for now we will 
derive simple scale-free estimates of the mass profile and 
compare them to previous results. Our dynamical anal- 
ysis is presented in £17.11 comparisons to other studies in 
£17.21 and implications for the dark matter halo in £17.31 

7.1. Scale-free Analyses 

Our dynamical analysis begins with a method descrip- 
tions in £17.1.11 our results in £17.1. 2[ and checks of the 
methods with other data sets in £17.1.31 

7.1.1. Methods 

The first key approximation we will make is of dynam- 
ical equilibrium, which is certainly questionable in the 
light of the indications of substructure found in M87, 
and will be an issue for even the more advanced dynam- 
ical models. Nevertheless, we are restricted at this point 
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to simple assumptions, which will also include spherical 
symmetry. 

The basic Jeans equation describing the relation be- 
tween the gravitating mass profile and the kinematics 
and density of a tracer population can then be expressed 
in the simplified form: 



10 R [kpc] 100 



1000 



v 2 c (r) = Hr)-2P(r) +1 (r)}a 2 r (r), 



(6) 



where r is the three-dimensional galactocentric radius, 
a(r) = ~d\\\v / dhir is the slope of the tracer density 
profile v(r), f3(r) = 1 — (Tg/cr 2 . is the anisotropy parame- 
ter that describes the balance of the tangential and radial 
components of the velocity dispersion ag(r) and o~ r (r), 
7(r) = — d\na 2 /d\nr is the slope of the internal velocity 
dispersion, and the circular velocity is related to the cu- 
m ulative mass profile by v 2 ( r) = GM(r)/r (cf Eq. 4.215 
of lBinnev fe Tremaindf2008h . 

One advantage of the circular velocity expressed as a 
function of angular radius is that once derived, this quan- 
tity is independent of the distance to the galaxy and of 
any luminosity model and filter bandpass characterizing 
the stellar distribution. Note also that unlike the case of 
a thin-disk galaxy, the luminosity and basic velocity pro- 
files are not enough to specify the mass profile because 
of the additional factor (3. This is the "mass-anisotropy 
degeneracy" that is the bugaboo of studying elliptical 
galaxy dynamics and which may be alleviated through 
several techniques, some of which we will incorporate 
here. 

We next adopt the simplification that the galaxy's dy- 
namics are scale-free, i.e. the quantities a, /3, 7, and v c 
in equation [5] are all independent of the radius, which 
also implies that the projected velocity dispersion and 
kurtosis profiles a p (R p ) and k p (R p ) are scale-free. The 
conversion between observed dispersion and circular ve- 
locity can then be expressed as a constant: 

v c (r = R p ) = kap(Rp), where 

,2 _ (g + 7)(« + 7-2/3) T[(q + 7 )/2]r[(a - l)/2] 
a + 7 - (a + 7 - 1)0 r[(a + 7 - l)/2]r[a/2] 

, , , (7) 

(IDekel et al.l 120051: see also lEfstathiou et all Il980t 
IGerhardl 119931: IWatkins et al. 112010). Note that the ob- 
served dispersion er p we will use here actually has the 
(dynamically weak) rotation folded in as v vms . Also, real 
galaxies and GC systems are not exactly power-law, but 
we will adopt the approximation that the slope of a p at 
R p provides the line-of-sight averaged slope 7 for ay at 
r = R p , and similarly, the slope of the surface density 
profile yields the effective a after adding 1. 

Some interesting cases are an isotropic system {f3 = 0) 
where the pre-factor in k 2 reduces to (a + 7); a constant- 
dispersion system (7 = 0) where k 2 — a(a — 2(3)/(a+(3 — 
a/3); and an isotropic constant-dispersion system where 
k 2 = a. Note also the complicated dependence of k on 
the parameters. For example, the classic situation where 
radial anisotropy suppresses a p in galaxy halos (i.e., k 
increases with 0) applies only when {a + 7) > 3; the 
reverse is true for shallower density profiles with flat ve- 
locity dispersions such as the case of the observed M87 
GC system, where radial anisotropy should actually boost 
the dispersion (k decreases with /3). This also means that 
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Fig. 24. — Root-mean-square velocity profiles for subpopulations 
in Virgo, relative to the systemic velocity of M87. The different 
subpopulations are indicated in the panel legend; the stellar kine- 
matics are from VIRUS-P, using true second-moment estimates 
provided by J. Murphy & K. Gebhardt. The differences between 
GC subpopulations are relatively minor on this scale, and we lump 
all of the GCs together. The open and filled green squares show 
the PNe with high-velocity "outliers" alternatively excluded or in- 
cluded. 

in regions where (a + 7) ~ 3, the inferred mass will be 
fairly insensitive to the (often uncertain) value of (3; for 
M87, this should happen at ~ 90 kpc for the total or blue 
GC system, and at ~ 25 kpc for the red GC system. 

An additional modeling ingredient is the use of the ob- 
served velocity kurtosis k p to estimate the anisotropy /?. 
INapolitano et al.l (|2009l ) found that for a certain class of 
spherical systems with constant profiles of (3 and velocity 
dispersion (7 = 0),/? can be derived geometrically by a 
set of projection integrals (see their equations BIO B12). 
While all of these assumptions will be at some level in- 
accurate for M87, we can use this approach to derive a 
first plausible guess for (3. 

7.1.2. Results 

Now we carry out the dynamical analysis using all of 
the M87 GCs from the "new" data seQ We show the 
corresponding rms velocity profile in Figure [24] (which 
includes other details that will be discussed later). Before 
binning the data with radius, we illustrate the method 
with the data combined. The GCS dispersion is fairly 
constant overall, with 7 = 0.13 ± 0.07 and a p — 320 ± 
11 km s _1 at a median radius of 6.9' (33 kpc), while 
the kurtosis is k p — 0.32 ± 0.24. We then infer j3 ~ 
0.3 ± 0.2 and after substitution in equation [7] we find 
k ~ 1.64 ± 0.09 and therefore v c ~ 525 ± 28 km s" 1 
(allowing for the statistical uncertainties in a, /3, 7, and 
ct p)- 

30 401 objects, where we have omitted the IGCs, UCDs, bright 
objects, the possibly dwarf-bound GC H35970, and the weird, high- 
velocity object S923. 
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Fig. 25. — Mass profile of M87, expressed as circular velocity 
versus radius. Orange filled points with error bars show results 
from this paper based on simple dynamical modeling of the new 
GC data; the equivalent modeling using the outer stellar kinemat- 
ics data of VIRUS-P is shown as a blue star symbol. Curves with 
shaded regions show models from the literature based on X-ray 
data (Jpas ct al. 2010; very narrow uncertainty band) and old GC 
data I IMurphy et al.H2011l) . A constant mass-to-light ratio model 
is shown as a dashed curve for compariso n. Green dashed an d or- 
ange dotted curves show m a ss mo dels from[McLaughlin (1999) and 
IRoman owskv & Kochanck (2001), respectively. Th e red triangle 
with e rror bars shows another GC-based result, from lDeason et alj 

poTl) . 

The inclusion of S923 would have boosted both the 
dispersion and kurtosis, and implied a higher anisotropy 
/3 ~ 0.7 ± 0.1: a large change induced by just one ob- 
ject out of - 40CE3. The net effect on the v c estimate 
would be very small because the dispersion and kurtosis 
work in opposite directions in the central regions of the 
GCS. However, the same would not be true if such an 
object were found in the outer regions, highlighting the 
critical importance of contamination rejection in discrete 
velocity samples. 

We next follow the same technique, while breaking the 
dataset down into three radial bins with ~ 100-150 GCs 
each. The results are shown in Figure [35] the circular 
velocity shows some indication of a decline with radius 
but is consistent with being constant at v c ~ 530 km s _1 
over the range of radii probed by the GCs. The large un- 
certainties in each radial bin are driven by the uncertain 
local dispersion slope 7. 

We leave the next steps to future work, including sepa- 
rate examination of the dynamics of GC subpopulations 
in color, luminosit y, and size, and the use of more de- 
tailed models (e.g., IMurphv et al.ll20Tll) . 

31 Investigating the dynamics of the subpopulations in more de- 
tail is beyond the scope of this paper, but if we included the UCDs 
and bright objects in our analysis, then the results j3 and v c results 
above would not be affected significantly. We also infer the blue and 
red GCs to have overall /3 = 0.2 ± 0.3 and 0.6 ± 0.2, respectively — 
although these values should not be compared directly since they 
are measured at different characteristic radii. 



7.1.3. Sanity Checks 

Our final step here is to carry out two checks on the 
reliability of our methods. The first is to apply the same 
simple modeling techniques to the stellar kinematics data 
from VIRUS-P (introduced in gOjJ). Using the data 
from the last three points (where /14 is approximately 
converted to n p ), we derive an estimate of v c ~ 524 ± 
14 km s _1 at ~ 14 kpc, as shown in Figure [25] 

This v c v alue is almost identical to the full 
IMurphv et alj (|2011[ ) modeling results at the same radius 
(where the dispersions of the stars and the old and new 
GC data are all very similar). Our inferred /3 ~ — 0.5±0.1 
for the stars can be compared to the Murphy et al. find- 
ing of /3 = — O.llgg, which should be fairly insensitive to 
any errors in the mass profile caused by the old GC data, 
given the slope a ~ —2.8 of the luminosity profile in this 
region, which should mean that the anisotropy is driven 
by the velocity distribution shape. These anisotropy re- 
sults are compatible but the formal error bars in our 
method are clearly missing some sizable systematic un- 
certainties. 

The second check is to apply our simple model to the 
old GC data, to verify that our new finding of decreased 
mass is caused by the data, not the dynamical meth- 
ods. However, the steep outer dispersion rise in the old 
data turns out to be problematic for our scale- free model, 
apparently allowing no physical solution for k with near- 
isotropic orbits, while slightly more anisotropic assump- 
tions formally permit any value for k. Put another way, 
the "observation" 7^—3 implies a mass density profile 
that impossibly increases with radius, p(r) oc r. This 
is a warning of the limitations of scale-free models, and 
also suggests that sharp increases in galaxy velocity dis- 
persion may signify problems in the data (note though 
that the rising stellar dispersion profiles found in some 
cD galaxies are not as steep as in the old M87 GC data) . 

Some benchmark values to report from our models are 
a dynamical mass of M = (9.2 ± 4.0) x 1O U M and 
(4.3 ± 1.1) x 10 12 M Q at r = 11 and 84 kpc, respectively. 
The cumulative mass-to-light ratios are then = (20 ± 
9)Tb,q and = (46 ± 15)Ts i0 , respectively. 

7.2. Comparisons 

In Figure [25] we also show for comparison some mass 
profiles derived from the literature. The first ones we will 
discuss come from state-of-the-art anal ysis of the latest 
M87 data till now: X-ray gas emission (IDas et al.H2010D 
and c ombined stellar and GC kinematics (Murphy et al. 
I2011f h Both techniques found a circular velocity that 
increases outside ~ 10 kpc, but with inconsistent ampli- 
tudes. 

Our new results at large radii are incompatible with 
any reasonable extrapolation outwards of the Mur- 
phy et al. models, which would have three times as much 
dynamical mass within ~ 85 kpc as our new estimate. 
This is not surprising since those models were based on 
the old GC kinematics data set with a high outer disper- 
sion which we no longer find (Figures [TBI and [24]) . There 
may also be inconsistency of our mass results with the 
outer X-ray profile, but we do not regard this as a firm 
result given the systematics in our dynamical analysis. 

Three other sets of models are also shown in Fig- 
ure [25] One is a mass model constructed for Virgo by 
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iMcLaughlinl (|1999T) that was later found by iCote et all 
( 20011) to reproduce the GC dynamics if roughly isotropic 
orbits were adopted. The three dotted curves are sample 
models from the joint dynamical analys is of stars and 
GCs in lRomanowsky fc Kochanekl (pool . The red open 
triangle is from a constant-an isotropy powe r -law a nal- 
ysis of the GC dynamics from iDeason et ail (|201lD . It 
can be seen that all these models, which used the old GC 
data set, are generally close to our new results using the 
new GC data, and dissimilar to the Murphy et al. model 
based on the old data. 

These comparisons are puzzling. It is possible that 
some of the models got the "right" answer for the wrong 
reasons, but we did notice in a broader inventory of pre- 
vious models of M87 that a surprisingly wide range of 
mass results were obtained from different studies that 
used basically the same data sets. This suggests that 
some methods may be less sensitive to errors and out- 
liers in the data than others, which is a possibility that 
should be tested further through modeling of simulated 
data sets. 

7.3. Dark Matter Halo 

We next consider some general implications that can 
be drawn for the dark matter (DM) surrounding M87 — 
without doing any additional dynamical modeling. We 
begin with the assumption that the Murphy et al. results 
are robust within the region probed by the stellar dy- 
namics (r ~ 14 kpc) and that our new mass constraint 
at ~ 85 kpc is accurate. We then explore a range of 
cosmologically-motivated mass models that may be com- 
pared to these constraints. 

Our basic model consists of a stellar mass distri- 
bution and a ACDM halo. For the former, we use 
the tr iple power-law model of lRomanowskv fc Kochanekl 
(|2001| ) based on the B-band surface brightness profile: 
although this does not use the most modern photometric 
results (e.g., iKormendv et al.l 12009). it is good enough 
for our purpo s es. F or the latter, we adopt the classic 
iNavarro et all (|1997t l profile: 
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where p s and r s are the characteristic density and scale 
radius. The overall model is specified by the latter pa- 
rameters along with the (constant) stellar mass-to-light 
ratio T* : s. 

Given three free parameters, it is fairly easy to fit a 
wide range of constraints. However, there are strong 
prior probabilities on all of these parameters. In par- 
ticular, there is a statistical relation expected between 
p s and r s (or alternatively, between h alo virial mass M v j r 
and concentration c V i r = r vir /r s ; e.g. JPrada et aljfeoill ). 

As a first rough guess for M v j r , we may use a sim- 
ple relation between velocity dispersion and virial ra- 
dius, based o n the approximation o f a singular isothermal 
sphere ( after ICarlberg et "all 1 19971) : 

(9) 
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We consider two alternatives: that the DM halo of 
M87 is smoothly contiguous with a relaxed Virgo clus- 
ter halo; or that M87 hosts its own group-scale halo 



that is still decoupled from the larger Virgo environ- 
ment. In the first case, we use the cluster galaxies' 
overall cr p = 803 ± 29 km s" 1 to find r vir ~ 2.8 Mpc, 
with a corresponding Af v i r ~ 1.2 x 1O 15 M . For com- 
pariso n, more detail e d ana lyses of the V i rgo ki nemat- 
ics by iFouque et aD (pOOll) and by iTuliyl (|2005l ) imply 
M vir - 1O 15 M , and by lRines fc Diaferid $2QQB) implies 
r vir ~ 1.9 Mpc and M vir ~ 4 x 1O 14 M . 

For the second case, we use the GC velocity dispersion, 
which will provide a lower limit to the mass, since the 
GCs are probably colder than the DM halo itself. We find 
r vir £ 1.1 Mpc and M vir £ 8 x 10 13 M Q . The - 2.3 keV 
temperature of the X-ray gas around M87 also suggests 
r vir ~ 1.4 Mpc and M yir ~ 1.7 x lO 14 -^ from standard 
scaling relations fe.g.. lUrban et al.ll201lD . 

For comparison to these dynamical metrics, classi- 
cal number-density surveys of the Virgo cluster at op- 
tical and X-ray wavelengths characterize it as a region 
of ~ 2 Mp c projected radius with distinct substruc- 
tures (e.g.. iBinggeli et all 119871 iBohringer et al.l 119941: 
iSchindler et all Il999l ). For example the "BCG" M49 
(NGC 4472) hosts a subcluster "B" at 1.3 Mpc projected 
distance from M87, while the overdensity "A" associated 
with M87 itself may be characterized as a region of ~ 0.6- 
1.0 Mpc projected radius. Given the overall constraints, 
we consider r v - 1T ~ 1.2 Mpc and M V i r ~ 1O 14 M to be a 
reason able model for the subhalo around M87. 

The McLaughlin (1999) model already shown in Fig- 
ure [H has r vir ~ 2.2 Mpc, M vir ~ 6 x 1O 14 M , and 
c v ir — 4. This would be a fairly reasonable model for 
a relaxed, cluster-wide halo, but as the Figure shows, it 
does not accommodate as much mass in the ~ 20 kpc re- 
gion as currently inferred. In fact, it turns out that this is 
a generic problem, and no standard halo (of either group 
or cluster mass) fits our constraint^. 

There are several possibilities for a DM halo with a 
higher central density. One is a halo with a much higher 
concen tration than average; this i s show n conceptually 
by the Romanowsky & Kochanek (2001) model in Fig- 
ure [25] that comes close to matching the "data" . The im- 
plication would be that M87 is at the center of a group 
halo that collapsed at an earlier time than average for its 
mass. 

Another possibility is that a process like adiabatic 
contraction has raised the central halo density (e.g., 
iBlumenthal et alJ 119861). Alternatively, a "cored" 
isothermal halo (e.g., INapolitano et all 120111 ) would 
work, but this is not surprising given the extra free pa- 
rameter with no prior constraint, not to mention the lack 
of theoretical motivation for the profile itself. Finally, 
there is the possibility that our simple dynamical models 
are wrong, and that the mass at ~ 85 kpc is somewhat 
higher, allowing for a normal M V - 1T ~ 1O 14 M halo. 

Unfortunately we cannot resolve these possibilities 
without further extensive modeling. For now, we con- 
clude that the properties of the DM halo(s) of M87 and 
Virgo are poorly known — an unfortunate circumstance 
for the nearest galaxy cluster which should be reme- 

32 In lSpitler et alJ Q201H 1 we focused on modeling the "blue" GC 
subpopulation and found that if the central DM density cusp has a 
log slope of —1.2 rather than —1 as in the NFW profile above, then 
a halo with M v i r ~ 5 X 1O 13 M0 can fit the data, using a plausible 
/3(r) anisotropy profile. 
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died as soon as possible through more detailed model- 
ing. In particular, the DM density is of great inter- 
est in potentially providing observable gamma-ray sig- 
natures from the D M pa r ticles (e.g.. iGorchtein et alJ 
20101: iCuesta et all 120111 : iSanchez-Conde et all 120111 : 
Saxena et al.H201lT ). 



8. DISCUSSION 

We now bring together the preceding compendium of 
new observational results on M87 into discussions of sev- 
eral thematic areas. In §8.1l we address the discrepancies 
between our new GC rotation and velocity dispersion 
profiles with those from the literature. <JH2] discusses the 
evidence for multiple GC subpopulations and their re- 
lationships to individual field star components. In i j8.3l 
we delve further into the kinematical properties of can- 
didate UCDs and weigh the evidence for different for- 
mation scenarios. §8.41 is an analysis of the reality of 
proposed structural transitions in M87, especially in the 
halo, and the related issue of intergalactic GCs. In M8.5I 
we consider various possibilities for recent mergers. In 
£18.61 we present a detailed comparison of M87 with the 
cluster-central Fornax elliptical NGC 1399. Finally, in 
§8.71 we discuss the larger implications of this work for 
understanding the formation of massive galaxies. 

8.1. Rotation and Dispersion Revisited 

One of our key results is that we do not reproduce 
the higher outer rotati on and dispersion previously found 
for GCs in M87 (e.g., iCote et all [20011 ). Neither do we 
find any evidence for a steeply falling dispersion profile 
indicative of a truncation of the stellar halo at a radius of 
~ 150 k pc, as suggested by k inematics of a small sample 
of PNe (|Dohertv et al.|[200l ). Instead, we find a fairly 
constant velocity dispersion profile out to a projected 
radius of nearly 200 kpc. 

In previous sections we have briefly discussed the ori- 
gin of the discrepancies with previous GC kinematics 
findings. Part of the problem appears to be the pres- 
ence of what we term catastrophic outliers; e.g., two 
GCs (S878 and S1074) had published radial velocities 
of ~ 2200 and ~ 2500 km s _1 that were erroneous by 
~ 900 to 1000 km s _1 . Both are located in the radial 
range 30-35 kpc, and thus contributed to the high disper- 
sion inferred (another relevant case is S7023, discussed in 
§4.4. 3p . However, other GCs with extreme velocities have 
been confirmed by multiple measurements (e.g., S66 and 
S176). Therefore, we cannot conclude that such outliers 
are the sole cause of the disagreement; additional dupli- 
cate measurements of older radial velocities are needed. 

Another factor that may be relevant is the luminos- 
ity difference between the samples. A subset of the 
more massive GCs may be "contaminating" the sample 
as UCDs with distinct kinematics (see below), and the 
median magnitude of the Cote et al. GC sample is ~ 0.8 
mag brighter than in our sample. 

Improved photometry is also important for more ac- 
curate kinematics because of the often underappreciated 
complication of high-velocity Galactic star contamina- 
tion in GC data sets. We have used both high-precision 
color-color diagnostics, and HST- measured sizes, to help 
w eed out stars and (t o a les ser extent) galaxies. 

Romano wskv et all ([20111 ) found that there is substan- 
tial kinematical substructure among GCs in the halo of 



M87, especially at radii > 50 kpc. Therefore we consider 
it plausible that another important reason for our dis- 
crepancy with previous work (see especially Figure II 6|) 
is that each study is sampling different regions of sub- 
structure. If true, this could indicate the dominance of 
substructure starting at even smaller radii (perhaps 30 
kpc) than currently demonstrated. Additional velocities 
in this "transition" region, with good azimuthal cover- 
age, are needed to assess the importance of substructure. 

8.2. Subpopulations: Bimodality and Beyond 

The prototype of GC bimodality is the Milky Way, 
where the disk/bulge and halo subpopulations are dis- 
tinguished by strongly different distributions in metallic- 
ity, rotation, and positions relative to the center of the 
galaxy. Similar studies of external galaxies are gener- 
ally confined to GC color, which is used as a metallicity 
proxy bec ause of the generally la rge ages found for the 
GCs (e.g- lBrodie fc Strad"er1l200l . 

Bimodal GC color distributions have been estab- 
lished for many galaxies including M87 (e.g., iPeng et all 
l2006f ) , but there are lingering controversies over whether 
color bimodality reflects a true underlying metal- 
licity bimodality, or a nonlinearity of the color - 
meta l licity relation (e.g.,lRichtlerll2006trYoon. Yi. fc Le< 
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more general analysis of GC chemo-dynamical phase- 
space would thus be helpful in resolving this matter 
definitively. 

In NGC 1399, ISchuberth et~aT1 dMlOft showed that 
there is a distinct offset in the GC velocity dispersion 
at the same color as the bimodal separation between the 
blue and red GC subpopulations. We have attempted 
a similar analysis in the case of M87 ( §6.3p . but our 
findings are less clear. Certainly we have found strong 
kinematical differences between different color subpopu- 
lations (e.g., Figures [2"T1 and |2"2"]) . but the picture is com- 
plicated by the evidence that we have unco vered for at 
least a third distinct GC sub-population (cf iBlom et alJ 
[2012] ). 

The third population has intermediate spatial distri- 
bution and colors with respect to the blue and red sub- 
populations, a stronger tail to high luminosities, and 
somewhat different kinematics. Curiously, this popula- 
tion appears to provide the best kinematical match with 
the VIRUS-P stellar kinematics, while the far-red GCs 
match well with the PNe. The radial density profile of 
the full subpopulation of red GCs is also a very good 
match to the F-band stellar surface brightness of M87. 

Such comparisons are important for understanding 
the origins and interrelations of different subpopulations 
within galaxies. The conventional wisdom about GCs 
is that the metal-poor subpopulation corresponds to the 
faint metal-poor stellar halo, while the metal-rich sub- 
population is strongly associated with the formation of 
the hot spheroid ("bulge"). Our results in M87 generally 
support this picture. 

The parent stellar populations of the bright PNe 
in early-type galaxie s are somewhat mysterious (e.g., 
iCiardullo et al.l I2005D . Empirically, the PN number 
densities and kinematics seem to generally agree well 
with the properties of the surrounding field stars (e.g., 
iCoccato et al.|[2009T) . In M87 we find a preliminary indi- 
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cation that the PNe are associated with the most metal- 
rich GCs. This can be tested in the future with a kine- 
matical sample of PNe at small radii, where the far-red 
GCs show a peculiar velocity offset relative to systemic. 

Detailed, orbit-based dynamical models will be needed 
to understand more clearly the relations between these 
different subpopulations. 

8.3. Bright GCs and UCDs 

When using GCs as discrete kinematical tracers, a 
rarely-considered point is the luminosity of the clusters. 
Brighter GCs are more frequently targeted, for obvious 
reasons, and this has probably caused significant biases in 
many of the GC kinematics studies to date. Not only can 
tide-driven evolution produce co rrelations between GC 
luminosity and kinematics (e.g., iVesperini et al.l F2003) . 
but there is now abundant photometric evidence for a 
transition between normal GCs and extended UCDs at 
the bright end of the "GC" luminosity function. 

Understanding the nature of UCDs and the degree to 
which they "contaminate" the normal GC population is 
a ongoing area of exploration, and here we add just a few 
lessons learned from our analysis of M87. The first is that 
a joint examination of size, luminosity, and color (£ 34.4. f \ 
suggests that there is a genuinely distinct population of 
extended objects that are not simply a continuation of 
th e bright GC populat ion. We develop this theme further 
in lBrodie et~atl (|2011l ). 

Next, we have found that the most luminous objects in 
M87 (io ^ 20) have kinematics distinct from fainter GCs 
( ^6.31 and ^6.4|l . including elevated velocity dispersions 
and peculiar radial velocity distributions that avoid the 
systemic velocity. Most of the bright objects are also ex- 
tended, and there is some evidence that size rather than 
luminosity may be the key parameter in these peculiar 
kinematics trends (although this conclusion hinges on the 
reliability of the older velocity data). This reinforces our 
conclusions from photometry that the UCDs and GCs 
are distinct classes of objects. 

These kinematical diff erences are part o f an em erging 
pattern discussed by iRomanowskv et al~l (|2009f ) where 
the brighter "GCs" (with sizes unknown) in massive el- 
lipticals show peculiar kinematics. This is typically man- 
ifested as a higher dispersion and/or double-peaked ve- 
locity distribution, with the transition at Mi ~ —10 to 
— f f . Additional, recent support for this picture has come 
from NGC 5128 dWoodlev et al.l l2010h and NGC 1399 
(jSchuberth et al.1 l2010f ). In the latter study, the blue 
GCs seem more strongly affected, which a lso appears to 
be the case in M87. On the other hand, iMisgeld et al.l 
poll found for NGC 3311 that velocity dispersion de- 
creases for the brighter objects. 

To interpret these findings, we consider two simple sce- 
narios. One is that the UCDs are the remnant cores 
of larger galaxies which have been stri pped through 
a process such as tidal threshing fe.g.. iBassino et all 
11994 iBekki et al.l I2001D . The other is that they 
are bona fide star clusters with large sizes from 
a range of possible causes, e.g., from birth in a 

dwarf galaxy environment or fro m m ergers of star 

clust e rs (e.g.. [Fellhauer fc Kroupal l2002t iBurkert et al 
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These two scenarios may be difficult to distinguish if 
they both lead to final sizes that are established by tidal 
limitations from the surrounding M8 7 gravity well (see 
further discussion in lBrodie et al.ll2011l ). However, there 
should be residual kinematical and dynamical signatures 
of the origins of UCDs. If they began as normal nucleated 
dwarf galaxies, then they have become stripped down 
to UCD size by making close passages to the center of 
M87. In this scenario, one would na turally expect them 
to be on preferent ially radial orbits (|Bassino et al.lll994l : 
IBekki et all I2003L although iGoerdt et al.l 120081 found 
more complicated orbital trends) . The UCDs as a popu- 
lation should then reside in a centrally concentrated dis- 
tribution, with radial velocities th at decline ste eply with 
increasing galactocentric distance (Bckki 2007]). 

In the second scenario, the current UCD sizes are sim- 
ilar to their original sizes when formed as extended star 
clusters, and they represent the surviving objects that 
have avoided plunging close to center of M87. The UCDs 
would then be expected to reside on more tangential or- 
bits, showing dispersions that are low in the center and 
increase outwards. The number density profile should 
also have a shallow core, with large objects found near 
the center only in projection. 

The implications for velocity distribution shape are less 
clear. Although radial and tangential orbits are classi- 
cally expected to produce "peaky" and "double-horned" 
velocity distributions, respectively, this is for objects fol- 
lowing a s mooth power-law density distr ibution with ra- 
dius (e.g., Ivan der Marcl fc Franxl 119931 ) . The velocity 
distributions might differ in a situation with a cored den- 
sity law caused by ongoing tidal disruption. For example, 
the shell of objects on radial infall around M87 has a ve- 
locity distribution that changes from peaky near its edge 
to double-horn ed at smaller radii (a "ch evron" pattern 
in phase-space: IRomanowskv et al.ll2l)TTl ). 

Very qualitatively, we can outline some simple 
geometrically-based expectations for a well-mixed, quasi- 
equilibrium population following a cored density law and 
with radially-biased orbits (the stripped dwarf scenario) . 
This population can be approximated as a superposi- 
tion of many shell structures, with near-zero projected 
velocities for objects near apocenter, and high veloci- 
ties exceeding the circul ar velocity v c for objects n ear 
pericenter (cf Figure 3 of [Mcrrificld & Kuiiken 1991). A 
broad range of velocities is expected, except in the core 
region — here there are few objects near apocenter, and 
thus only the high-velocity objects near pericenter are 
seen. Therefore at small radii, one might expect to see a 
shell-like diverging chevron pattern of velocities. 

For near-circular orbits (the star cluster scenario), we 
would again expect a broad distribution of projected ve- 
locities, but with the range extending no higher than 
~ v c . At small projected radii, the circular orbits would 
generally be moving in the plane of the sky, and so the 
observed velocities would decrease towards zero. 

In summary, the origins of UCDs could be revealed 
by consideration of several kinematic aspects: the slope 
of their projected velocity dispersion profile, the detailed 
shape of their radius-velocity phase-space distribution at 
small radii, and their maximum observed velocities. The 
final aspect may be difficult to test in practice, since a 
sufficiently accurate v c (r) profile may not generally be 
available for comparison. 
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Now considering the overall properties of the M87 sys- 
tem of UCDs, we find a fair fraction of these to be at 
large radii (~ 20% at ~ 100 kpc), whi ch may argue 
again st the galaxy-threshing scenario (as iThomas et al.1 
120081 concluded in Fornax). However, we have not yet 
attempted to construct the M87 UCD density profile, 
which would require a careful accounting of selection ef- 
fects on both size and velocity measurements. The chal- 
lenge of estimating sizes over a wide field of view is such 
that no existing wide-field UCD study of any system has 
yet achieved this goal. 

Turning to the M87 UCD kinematics, their projected 
velocity dispersion profile is high compared to the nor- 
mal GCs, but is nearly constant with radius (Figure l2Tj) . 
which does not seem consistent with either UCD sce- 
nario. As discussed in i)6.3.4[ the bright UCDs show a 
remarkable double-peaked velocity distribution at small 
radii, but their phase-space distribution shows both di- 
verging and converging behavior with radius (top left 
panel of Figure [20]), which suggests a combination of 
the behaviors predicted above for the two UC D scenarios 
(cf. t he top and bottom panels of Figure 1 in IRix et akl 
fl997h . There is one UCD with a velocity well in excess of 
v c , but otherwise the highest velocities are close enough 
to v c to require more detailed modeling for interpreta- 
tion. 

The overall kinematic picture of the M87 UCDs may 
argue for a mixture of objects with different origins, with 
additional data and orbital analyses needed for clarifi- 
cation. In NGC 1399, the UCD velocity dispersion was 
found to increase with r adius (once the central UCDs are 
inclu ded in the analysis; I Thomas et a l. 2008; Greg g et"aLl 
2009) , which could imply that the Fornax UCDs are pre- 
dominantly star clusters. 

8.4. Structural Transitions and Intergalactic GCs 

As discussed in SJTJ a key goal of this study was to 
search for transitions in the halo of M87 between the 
galaxy and the surrounding cluster. One would naturally 
look at the stellar surface photometry, which in many el- 
liptical galaxies does show sharp changes in isophote po- 
sition angle, ellipticity, or higher-order shape that surely 
reflect internal transitions in the orbital structure. 

Despite the extensiv e photometric work on M87 
over the years (e.g. , Ide Vaucouleurs fc Nietol 119781 : 
iKormendy et alJ 120091 : Uanowiecki et al.l 12010t h it has 
been difficult even to establish whether or not it hosts a 
cD envelope. We have made an inventory of the M87 data 
and found no obvious, agreed-upon photometric features, 
other than the strong increase of ellipticity with radius. 
This behavior ha s previously been interpreted as rota- 
tional flattening (Ki ssler-Patig fc GebhardU [1998) . but, 
given our revised rotational results, we suggest instead 
that it may be due to preferential accretion of material 
along a dark matter filament. 

We consider next the possibility of transitions in kine- 
matics and GC metallicities, reviewing our findings so 
far (recall that 1' corresponds to 4.8 kpc). There is a 
twist in the stellar kinematics at ~ 1.5 kpc, and changes 
in the colors and mean velocities of the red GCs inside 
~ 10 kpc, extending down to at least ~ 2 kpc. The 
blue GCs show a peak-color transition somewhere around 
~ 15-25 kpc. The velocity distribution of the UCDs and 
bright blue GCs has an unusual "double-peaked" shape 



over the range ~ 10-30 kpc. The far-red and blue GCs 
have dispersion drops at ~ 17 kpc and ~ 40-50 kpc, 
respectively. The 50 kpc radius seems to be a zone of 
multiple transitions, including the blue GC, far-red GC, 
and PN rotation and dispersion, and the UCD rotation 
and velocity distribution shape. 

This is a portrait of complex behavior that does not 
match up clearly with any transitions in the galaxy pho- 
tomet ry. We do not reproduce the result of ICote et al.1 
(|2001h that a twist in the BGC rotation coincides with an 
onset of the cD envelope at a radius of ~ 21 kpc. How- 
ever the ~ 50 kpc kinematical transition zone does cor- 
res pond with the onset of a halo substructure identified 
by iRomanowskv et al.1 (|2011l) with different techniques. 

As a final exercise we look for tra nsitions in t he far 
outer halo of M87. Given that (i) ICote et al.l (|2001[ ) 
found a strongly increasing outer GC dispersion that 
extrapolated smoothly to the hig h dispersion of cluster 
galaxies at larger radii; and (ii) iDohertv et al.l (2009) 
found a low outer dispersion in the PNe, along with a 
paucity of PN detections that they interpreted as a trun- 
cation of the stellar halo at ~ 160 kpc, it is of great 
interest to revisit this issue with our new data. 

As shown with the rms velocity profiles in Figure l24l 
we do not confirm either of the above claims. All of 
the subpopulations that we have studied around M87 
stay at w rms ~ 300-500 km s _1 , out to the largest 
radii we have measured — in contrast to the cluster 
galaxies ([Rines fc Gellerl [2008D which have u rms ~ 700- 
1000 km s _1 , even at radii where they overlap with the 
GCs and PNe. 

We have already seen that the high GC dispersion pre- 
viously reported may have been caused by problems with 
the older data, and we have discussed how the low PN 
dispersion was predicated on "outlier" removal. We find 
no evidence for a steeply declining dispersion in the GCs 
near 160 kpc. Similarly, the photometric number den- 
sity of metal-poor GCs shows no notable features beyond 
100 kpc (Figure O; a single Sersic profile is a good fit all 
the way to the edge of t he photometric sampl e (see also 
similar constraints from iTamura et ail [20 06b). We can 
therefore categorically state that there is no edge to the 
stellar halo of M87 around these radii. 

As discussed in Sj7l there is already ample evidence that 
the Virgo cluster consists of multiple subsystems which 
have yet to merge and relax. Thus it is plausible that 
the dark matter halo and GC system surrounding M87 
are basically decoupled from the greater Virgo environ- 
ment, which would explain the relative coldness of the 
GC kinematics. However, kinematical measurements at 
larger radii should pick up a transition eventually, either 
due to the dynamical interface between M87 and Virgo, 
or because of an intergalactic (IGC) population seen in 
projection. 

I GCs are known to exist in Virgo from photome- 
try (|Williams et alj|2007t iDurrell et al.11201 11 ) becoming 
domin ant at perhaps ~ 190 kpc from M87 (|Lee et al.l 
2010b). Although we have not probed this far out yet 
with kinematics, projection effects should make many 
IGCs appear at smaller radii, with some of them stand- 
ing out by their extreme velocities. 

To investigate this issue quantitatively, we estimate 
the number of IGCs expected among our "new" spec- 
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Fig. 26. — Distribution of radial velocities of objects around M87, in two radial bins (as labeled in the panels). Light blue, orange, red, 
and green histograms show data for stars, GCs, Virgo galaxies, and PNe, respectively; all of the histograms, except the stars and the GCs 
in the left panel, are rescaled by a factor of three for visibility. The PNe in the right panel include a few that extend to R p = 84'; also 
shown are three IGCs at larger radii. 



troscopic dataset of 485 objects. For each of these ob- 
jects, we calculate its a priori probability for being an 
IGC via the local ratio of IGC surface density (esti- 
mated to be 0.1-0.5 arcmin~ 2 in iTamura et aTll2006bl) . 
and the IGC+GC density using our Sersic model for the 
GCS density profile from £15.11 Summing over all objects 
and allowing for statistical fluctuations, we find that we 
should have obtained between 17 and 87 IGC spectra, 
with 90% of these occurring at R p Z 8' ( > 40 kpc). 

To estimate how many of these IGCs should have 
"Virgo-like" velocities that clearly distinguish them from 
GCs bound to M87, in Figure we plot histograms 
of observed velocities for various populations. First we 
show GCs, stars, PNe, and Virgo galaxies in the outer re- 
gions of M87 where we have GC spectroscopy (left panel, 
R p =10'-41'). For a larger sample of galaxies, we also 
consider their velocity distribution from a larger-radius 
bin (right panel), assuming this is similar to the region at 
smaller (projected) radius. We also show the PN veloci- 
ties at larger radii, which as mentioned previously, show 
a transition to Virgo-like velocities outside ~ 40'. 

If we define extreme velocities to be more than 
1000 km s _1 relative to v sys for M87 (i.e., v p < 
307 km s _1 or v p > 2307 km s _1 ), then the various sub- 
samples of galaxy and intergalactic PNe imply an "ex- 
treme" fraction of between 16% and 67%. Therefore our 
spectroscopic data set should include between 3 and 58 
IGCs with extreme velocities. 

We do find 4 objects (including the weird S923) with 
extreme velocities (see Figure [T2l , but these are all at 
small radii (8-30 kpc), and except for S923, are only 
slightly in the "extreme" range. Therefore we deem them 
very likely to be objects bound to M87, seen around the 
pericenters of near-radial orbits. There are also several 
objects (GCs and PNe) at large radii (90-130 kpc) with 
relative velocities close to —1000 km s — 1 . As previously 
discussed, these are somewhat ambiguous: whether or 
not they are bound to M87 is highly dependent on the 
uncertain distribution of mass outside these radii. 

There could be additional lower-velocity GCs lurk- 



ing in our sample of foreground "stars" , many of which 
do not have size measurements to verify their classifi- 
cations. Indeed, there are three confirmed GCs in the 
older velocity data with reported relative velocities of 
around —1200 km s _1 (see ^4.4.21) . However, we also 
expect a significant fraction of IGCs to be found with 
relative velocities below ~ —1500 km s _1 and above 
~ 1000 km s _1 , which are not yet observed. 

Therefore we conclude that the IGCs around the core 
of Virgo have either a surface density of ~ 0.1 arcmin -2 
or lower, or a surprisingly cold radial velocity dispersion. 
The latter possibility is supported by the three bona fide 
IGCs with spectroscop y (at ~ 800 kpc distances from 
M87; lFirth et al.ll2008f) . Their velocities are close to v sys 
for M87 (see Figures [T^] and , so we speculate that 
much of the intergalactic material in the core of the Virgo 
cluster may be relatively cold and flowing in roughly the 
plane of the sky. One might suppose such an arrange- 
ment to be a co nsequence of filame ntary accretion into 
the cluster (e.g.. iKnebe et al.ll200l) . but the Virgo fila- 
ment is t hought to run al ong rather than across the line- 
of-sight (|Mei et al J 12007( 1. i.e., probably not coincident 
with the major axis of M87. 

8.5. Recent Merger (s)? 

The large substructure outside ~ 5 kpc is discussed 
in detail in iRom anows kv et al.l (|201lH . where the over- 
all conclusion is that a massive galaxy was accreted 
less than ~ 1 Gyr ago. Here we have also found ex- 
amples of small pairs and groups of objects ( §6.3jl that 
could be relics of disrupted dwarf gala xies. Similar fea- 
tures have been found in o ther systems dCote et alj|2003t 
IRomanowskv et al.l 120091 : iWoodlev fc Harris! l20Tlli . and 
it should be kept in mind that localized peculiarities in 
the GC kinematics might be caused by small substruc- 
tures with coherent kinematics. 

We now review additional signs of more significant in- 
teractions or mergers in M87, starting with the central re- 
gions. There is a known velocity of fset in the M87 stellar 
kinematics at a radius of ~ 0.2 kpc ( Dre ssier fc Richstonc 



M87 Globular Cluster Kinematics 



37 



19901 IJarvis fe Melnickl U99lt IJarvis fe Peletierl HMD: 
Carter fe Jenkinsl fT992). This feature has been suggested 



as the nuclear remnant of a smaller accreted galaxy that 
settled down to the center of M87 by dynamical friction. 
We have also found (8 36.3P that the far-red GCs show pe- 
culiar color and velocity shifts inside ~ 10 kpc, suggesting 
this region is unrclaxed. 

There is a well known central jet that is estimated to 
have been active for ~ 0.1 Gyr (|Owen. Eilek. fe Kassiml 
2000), and which could have been fed by gas in a merger. 
There is furthermore a complex of dusty and filamentary 
warm and hot gas extending out to radii of ~ 15 kpc that 
has b een suggested as a possible byproduct of a merger 
(e.g.. iSparks et all 12004 ) . Additional "cold fronts" are 
found in the hot gas at ~ 30 and 90 kpc which have 
been explained as gas sloshing provoke d by a "fly-by" 
of a m assive galaxy group ~ 1 Gyr ago (jRoediger et al.1 

HSU). 

It was also recently discovered that the supermassive 
black hole in M87 is slightly off-center, with one possi- 
ble expla nation being a galaxy merger within the past 
- 1 Gyr (|Batcheldor et all 12010) . 

Previous work on GCs in M87 derived high outer ro- 
tation, with a major merger as one possible explanation. 
However, we have found that this rotation detection was 
probably spurious. We find no strong signatures of rota- 
tion. 

lArp fc Bertolal (fl97ll) and 

IWeil. Bland-Hawthorn, fc Malinl ([19971 ) identified a 
broad fan of material at a radius of ~ 50 kpc. Weil et al. 
modeled this as the product of a smaller galaxy accr etion 
within the past ^0.5 Gyr. Uanowiecki et al.1 (|2010l ) cast 
doubt on the existence of this structure, but did identify 
other small ste llar substructures i n the outer halo o f 
M87 (see also IMihos et al.1 [2005t IRudick et all l20Toh . 
IRomanowskv et al.l ( 20111 ) identified a small moving 
group of GCs at ~ 150 kpc that may be associated with 
one of these stellar features, and could in principle be 
part of a substructure extending inwards. 

As previously mentioned, several dE galaxies are found 
at ~ 40 kpc distance from M87 (NGC 4476, NGC 4478, 
NGC 4486A, NGC 4486B, and IC 3443), most of them 
with signs of severe disturbance as might be expected 
in these regions. These include peculia rities in the 
stellar isophotes, colors, and kinemati c s (iRoodl 119651: 
Faberl [19731: ISandage fc Binggelil [l98l IPrugniel et al l 
1987t Ivan den Bosch et al.l 11994 lHallidav et all 120011: 
Lucero et all 120051: iFerrarese et all 120061) . Any or all 



of these galaxies could have disturbed the kinemat- 
ics of the GCs in this ~ 40 kpc region, both from 
gravitational scattering of M87's pre-existing GCs, and 
(more likely) from depositing their own stripped GCs 
on unmixed, coherent orbits (see, e.g., discussion in 
IKissler-Patig fc Gebhardtl [19911 . 

The overall picture in M87 to date provides various 
indications for interactions with galaxies from the sur- 
rounding environment, probably including multiple inde- 
pendent events. Until now, these indicators have all been 
relatively subtle, with the overall visible picture of M87 
suggesting a placid galaxy where any recent interactions 
were presumably weak, involving accretion of low-mass 
galaxies. 

Focusing on the inner substructure that we have de- 
tected at radii of ~ 4-9 kpc, the crossing time here is 



~ 10-20 Myr, so assuming that the merger happened a 
few crossing times ago in order to still be very visible 
in phase space (Figure I19p but not in real space, we in- 
fer an event time no more than ~ 0.1 Gyr ago. This 
agrees with the timescales from the black-hole offset and 
central jet studies, and it may be that all of these fea- 
tures were caused by a significant gas-rich galaxy accre- 
tion ~ 0.1 Gyr ago that is still in the process of settling 
to the center of M87. 

There is a snag with this interpretation. The far-red, 
probably metal-rich GCs that show the velocity "slosh- 
ing" are unlikely to have been brought in by the accreted 
galaxy (which should typically host bluer clusters), but 
would instead be central GCs of M87 itself. Even if the 
center of M87 were somehow moved around enough by 
the merger to show a velocity offset, the problem is that 
the same offset is not seen in the other central subpopu- 
lations of M87: the blue GCs and (particularly) the field 
stars. More data and simulations for this central region 
should help clarify the situation. 

8.6. Comparisons to NGC 1399 and Beyond 

We now have the opportunity to draw illuminating 
comparisons and contrasts between the halo kinemat- 
ics of M87 and NGC 1399, the central galaxy in the 
second nearest cluster (Fornax). An overview of the 
two systems is shown in Figure 1271 Our new combined 
data set of GC radial velocities around M87 has 737 ob- 
jects including 487 with high-resolution measurements 
(18 km s _1 median uncertainty), out to galactocentric 
radii of nearly 200 kpc. The combined NGC 1399 dataset 
is comparable, with 729 objects total and 527 "class 
A" measurements ( 30 km s" 1 median uncertainty), also 
out to ~ 200 kpc (iRTchtler et al.l 120041: iBergond et all 
I2007t iSchuberth et al.ll2010l ). The median luminosities 
for the high-quality GC measurements are also similar. 
NGC 1399 has a m uch larger PN dataset than M87 
(|McNeil et al.l[20Tol) . 

One of the first striking differences between these stud- 
ies is that in M87 we had to eliminate only one GC from 
our sample as potentially bound to an adjacent dwarf 
galaxy, while NGC 1399 saw a great deal of "interloper" 
removal with dramatic ramifications for the kinematical 
and dynamical results. As we have already commented, 
overjudicious outlier trimming can be treacherous in the 
halos of galaxies, which may very well contain radially- 
biased orbits with velocities extending well past the nor- 
mal wings of a Gaussian velocity distribution. 

Some of the differences between the galaxies that 
we discuss below might be caused just by disparate 
treatments of o utlier s. Also, as emphasized by 
ISchuberth et~aT1 {2008), a very small handful of erro- 
neous velocity measurements (with presumably marginal 
S/N) can dramatically change the kinematics results, 
and so we consider only conclusions that are based on 
class-A measurements around NGC 1399. Even so, there 
were two cases of high-quality repeat velocity measure- 
ments with catastrophic (~ 4 a) differences. 

Like M87, NGC 1399 had earlier reports of very high 
outer rotation and velocity dispersion (in its PN kinemat- 
ics) which were overturned by the later higher-quality 
data. Both galaxies are now known to have extremely 
low rotation amplitudes overall, ~ 20 km s _1 . 

NGC 1399 has a peak central stellar velocity dispersion 
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Fig. 27. — Halo kinematics data in M87 (left) and NGC 1399 (right). Orange circles show GCs and green open squares show PNe (in the 
case of M87). The known "UCDs" are included for M87 but not for NGC 1399. A quality criterion has been used for both datasets: post- 
2003 measurements for M87, and "class A" for NGC 1399. A brown open circle in the NGC 1399 panel shows the approximate phase-space 
position of the satellite galaxy NGC 1404, whose accompanying GCs we have made no attempt to exclude. The equ ivalent radius R m f or 
the NGC 1399 data was calculated approximately by using an overall position angle of 90° and cllipticity of e = 0.2 (Dirsch ct al. 2003). 



1 which declines steadily with radius and 
200 km s _1 (as traced by the PNe and 



of - 370 km s" 
by ~ 50 kpc is 

RGCs). This is very similar to the behavior in M87, 
where we have also found that the ~ 50 kpc coldness 
may be related to a massive substructure, with the dis- 
persions rising somewhat at larger radii. The NGC 1399 
BGCs have an initially higher dispersion that also drops 
to ~ 200 km s _1 at ~ 150 kpc. 

Thus NGC 1399, like M87, shows no transition to 
hotter cluster ki nematics (which would be ~ 300- 
400 km s _1 ; e.g., iDrinkwater et al1l2001f ). In particu- 
lar, there is little sign of a transition at ~ 60 kpc where 
llkebe et all (|1996l ) claimed to see a galaxy-cluster inter- 
face from X-ray observations. 

The amplitudes of the GCS projected velocity disper- 
sions in M87 and NGC 1399 are not straightforward to 
interpret without understanding the mass, number den- 
sity, and orbital anisotropy profiles in detail. However, 
it seems initially remarkable that the dispersions of the 
two systems are so similar while the overall mass of Virgo 
may be larger than Fornax by a factor of ~ 10 fe.g.. lTullvl 
2005) . This situation further bolsters our suggestion that 
the M87 GCs are associated with a group-mass subhalo 
rather than with the entire Virgo cluster. 

The NGC 1399 data were previously interpreted as 
showing a dynamically relaxed central galaxy accom- 
panied by a subpopulation of vagrant' GCs (see refer- 
ences above), which are some combination of IGCs (i.e., 
associated with the entire Fornax cluster) and GCs on 
tidal streams stripped from other galaxies by NGC 1399. 
Such conclusions were motivated both by fluctuations in 
the velocity dispersion and by extreme radial velocities. 
The latter point is apparent in Figure [27l where there is 
one low velocity at ~ 170 kpc, and an asymmetric low- 



velocity tail at ~ 10-80 kpc (which is also seen in the 
PNe). 

M87 is generally more symmetric than NGC 1399 in 
phase-space but harbors signatures of cold substructure 
and numerous extreme velocity objects that could mark 
highly eccentric orbits. Both systems also show strong 
fluctuations in the mean velocities of their RGCs at small 
radii that imply they are out of equilibrium. The outer 
kinematics in both cases also suggest significant infall 
and accretion, even though no major disturban ces are 
apparent in optical imaging (e.g.. lTal et al.lr2009D . More 
work is needed to characterize these processes quantita- 
tively. 

One area where the systems may differ is in their ve- 
locity distribution shapes (considered after removing the 
bright GCs with their potentially different kinematics). 
NGC 1399 was found to have a near-zero kurtosis for the 
GCs overall, becoming slightly negative for the BGCs. In 
M87, we found generally positive kurtosis for both BGCs 
and RGCs. These differences are preliminary, pending 
uniform analysis of both datasets with the same outlier 
treatments, and considering also the variations of kurto- 
sis with radius. 

Keeping these caveats in mind, the kurtosis implica- 
tion would be that the GC orbits are radially-biased in 
M87, but isotropic or tangential in NGC 1399, perhaps 
implying formational differences between the two galax- 
ies. Radially-biased orbits are a fairly generic expecta- 
tion for galaxy halos, whether because of cosmological 
infall or fall-back of tidal tails from major mergers, but 
a transition to isotropy may occur within the central re- 
gions of massive groups or clusters during their initia l 
"fast accretion" phase (e.g.. iBiviano fc Poggiantil l2009). 
Therefore NGC 1399 may be in a more advanced state of 
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assembly than M87, consistent with the general idea that 
the Fornax cluster is more dynamically evolved than the 
Virgo (see also below). 

Other comparisons between the systems, concerning 
their luminosity-dependent GC kinematics, their UCDs, 
and kinematical e vidence for bimodality, w ere discussed 
in earlier sections. iSchuberth et al.l (|2010ft also found in 
NGC 1399 that the canonical connection between the 
RGCs and the visible starlight is supported by their sim- 
ilar density profiles, which in M87 we were able to estab- 
lish further through ellipticity comparisons (NGC 1399 
is nearly round so this test would be weaker). 

The Fornax cluster is optically more dense and sym- 
metric than Virgo, which is thought to refle ct a more ad- 
vance d stage of dynamical evolution (e.g., Uordan et al.l 
l2007f ). Two other nearby massive systems have also been 
recently studied via the dynamics of both their group 
galaxies and their GCs: the Eridanus A group (with 
central galaxy NGC 1407; e.g. . iRomanowskv et al1l2009ft 
and the Hydra I cluster ( with central galaxy NGC 3311; 
e.g., iRichtler et ail l201lh . The GC dispersion profile 
around NGC 1407 remains constant at ~ 250 km s _1 
out to at least 60 kpc, while in NGC 3311 there is a 
sharp dispersion increase (in both stars and GCs) that 
reaches ~ 700 km s _1 by 50 kpc. 

We propose that these four systems represent a quali- 
tative evolutionary sequence for clusters and their central 
galaxies. Eridanus A is a massive galaxy group that may 
be the progenitor of a cluster like Virgo which is in an 
active assembly phase. Fornax is a more evolved sys- 
tem, but like Virgo is still not internally relaxed enough 
to heat the central GC population or to accumulate a 
substantial cD envelope. Hydra I is closer to a "fully re- 
alized" cluster with a massive cD envelope in place that 
bridges the central galaxy and its host halo (although this 
region still does not ap pear well mixed in phase-space; 
Vcntimiglia et al 
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8.7. Implications for BCG and GCS Formation 

M87 is likely to remain a touchstone for formational 
tests of both BCGs and GCs, with the multi-dimensional 
properties of its field stars and GCs providing invaluable 
constraints. Here we will focus on a few specific aspects 
of the GCs: rotation, velocity dispersion, spatial distri- 
bution, and orbital anisotropy. 

8.7.1. Rotation 

The high outer rotation measured in the past for the 
GCs in M87, as well as in early studies with GCs and 
PNe in other galaxies such as NGC 1316, NGC 1399, 
and NGC 5128, has been interpreted as evidence for a 
major merger that transferred initial orbital a ngular mo- 
ment u m to the outer parts of the remnant dHui et al 



1995] lArnaboldi et all 119981; Ikissler-Patig fc Gebhardt 
1998t iCote et al.l 120011: iVitvitska et al.l 12002ft . At one 



point, this appeared to be shaping up as a new trend: 
"rotation may in fact be a common by-pro duct of the 
formation of GC systems" (|C6te et al.l 120031 ) . However, 
as we have seen in the case of M87, those earlier stud- 
ies seem to have been prone to systematically overes- 
timating the outer rotation, which with newer data on 
some of the same galaxies has turned out to be consider- 
ably lower (jWoodlev et all 120101: ISchuberth et all 120101: 



iMcNeil et alj |2~010). Other, more recent rotation studies 
could be affected by bias in t he analysis techniques as 
we have discuss ed here (e.g., IRomanowskv et ail 120091 : 
ILee et al|[20T0al ). 

Recent work on extended kinematic tracers in larger 
samples of early-type galaxies suggests that high outer 
rotation may be t he exception rather t han the rule (Proc- 
tor et al. 2009; iCoccato et al.l [2009h ■ and even in this 
context, M87 has gone from being one of the most dra- 
matic outer rotators to one of the weakest. Although 
more theoretical work is needed on the halo rotation of 
BCGs in different formational scenarios, low rotation as 
in M87 and NGC 1399 is qualitatively suggestive of the 
accretion of multiple, small systems with uncorrelated 
angular momenta rather than of a single massive merger. 

8.7.2. Velocity Dispersion and Assembly Constraints 

We have already discussed in ^8. 41 and ^8.6l the remark- 
able lack of a velocity dispersion increase with radius 
in M87, and the possible implication of a decoupling 
with the overall Virgo cluster. More generally, the ve- 
locity dispersion of a population of objects in a cluster is 
thought to reflect their time of infa ll, assuming they were 
born outside of t he cluster (e.g., iDiemand et alj [2~005t 
iMoore et al.ll2006l ). For example, spiral galaxies are gen- 
erally found with more extended radial distributions and 
correspondingly hotter velocity dispersions than ellipti- 
cals in cl usters, probably refle cting their more recent in- 
fall fe.g.. lTullv fc Shavaill98l . 

Qualitatively, if the GCs were accreted by M87 at an 
"early" stage, this idea could explain the disparity ob- 
served between the velocity dispersions of GCs and clus- 
ter galaxies even if there is no decoupling between M87 
and Virgo. An intriguing similarity can be seen between 
the GCs and the relatively cold kinematics of the sub- 
set of nucleated dEs around M87 with round isophotes, 
which has been in terpreted as a sig nature of an early 
formational epoch (jLisker et alll2009l ). 

If we assume first an extreme scenario where the GC 
system of M87 is entirely accreted, and that to a first 
approximation it is built up monot onically from the 
inside-out like a series o f tree rings (|Wang et all 1201 it 
iSalvador-Sole et al.ll20ll . then we can adopt a simple 
model to put limits on its assembly age. The basic idea 
is that the GCs observed now at some radius r cannot 
have arrived at a time when r <; 2r v j r (t wice the virial ra- 
dius of M87: IDiemand fc Kuhlenl 12008ft . Then using our 
mass model for M87 from S}7] (i.e., w c i rc ~ 500 km s -1 ) 
and cosmological form ulae for virial quantities (e.g., 
iBrvan &: Norm an 1998), we can map between r and some 
upper- limit on the assembly redshift zt(r). This is a con- 
servative upper-limit beca use it assumes the m ass within 
r did not grow with time (jCuesta et al.ll2008ft . 

Using this schematic, we find that the GC popula- 
tion (of all metallicities) at the ~ 100 and 200 kpc radii 
around M87 must have been accreted sometime after 
Zf ~ 8 and 5, respectively. These are not very strin- 
gent constraints, but they do illustrate conceptually the 
kind of inferences that one might make about the accre- 
tion history with more detailed models, e.g., based on 
the ellipticity, kinematics, and radial density profile of 
the GC system (IDiemand et al.ll200l IMoore etalll2006t 



Griffen et a l. 2010). We pursue an approach like this in 
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Another inference that one could make from standard 
ACDM cosmology is the characteristic radius of accret- 
ing material at the present day. Mate rial on its first 
infall typically penetrates to r ~ 0.3 r v ; r (jDiemand et al.1 
2007), so given our estimate in £)7.3I of r vir ~ 1.2 Mpc, 
we expect accretion around M87 to now be occurring 
typically at radii of r ~ 350 kpc to 2.5 Mpc (where the 
infalling galaxies should not be disrupted but might lose 
some of their outermost GCs to tidal stripping) . 

This expectation, along with the lack of a velocity dis- 
persion increase at large radii (suggesting little recent ac- 
cretion from the cluster environment), would appear to 
create some tension with the d iscovery of a massive accre - 
tion event at r ~ 50-100 kpc (|Romanowskv et al.ll201l[ ). 
The nominal dynamical friction timescales at 350 kpc are 
too long for even a 10 12 M Q subhalo to sink far inward, 
but it is plausible tha t an entire group of galaxies fell in 
and broke apart (e.g.. lRudick et al~ll2006| ). delivering the 
shell progenitor to small radii. 

We may apply similar arguments to the scenario at 
the other extreme, where the GCs are all formed in situ. 
Considering that the associated starbursts would proba- 
bly occur inside the virial radius, we can conclude that 
the formation of the M87 GC system as a whole, out 
to ~ 200 kpc, continued until after z ~ 2 (otherwise at 
higher z the GCs would have been forming outside the 
galaxy). Such timing would not sit comfortably with di- 
rect estimates of the GC ages, and would contradict the 
current observational picture of BCGs experiencing the 
bulk of their in situ starburst activity at z J> 4, with only 
fairly quiescent merging since then. 

This argument thus appears to rule out a wholly in situ 
formation scenario for the M87 GC system (and possibly 
for the cD envelope by association), with scope remain- 
ing for in situ formation within the inner ~ 100 kpc ra- 
dius. There is in principle a loophole where the in situ 
GCs could have formed early on in a compact configu- 
ration, and been propelled to larger radii much later by 
bombardment from infalling substructure. However, it is 
likely that the energy requirement of this "puffing up" 
would involve a large amount of infalling mass which 
would add its own considerable population of accreted 
GCs. 

8.7.3. Orbital Anisotropy 

A last clue to the formation of BCG halos is orbital 
anisotropy (also di scuss ed i n ^8.61) . As summariz ed by 
iHwang et all (pOfJl and lRomanowskv et al.l (f2009l ). M87 
and other high-mass ellipticals were previously found 
to host near-isotropic or even tangentially-biased or- 
bits that contradict general expe ctations from theory for 
radially-biased orb its (see also iNapolitano et all 120111 : 
iDeason et al.l 1201 1| ). An implication of these results 
might have been the existence of a process that drives the 
orbits toward isotropy in the central regi ons of groups or 
clusters (e.g.. iBiviano fc Pog gianti 2009). 

However, our new observational results overturn the 
previous findings of isotropy in M87, and we instead find 
indications of radially-biased orbits overall (©. We have 
not tracked down the exact reason for this difference, but 
it may be caused in large part by "contamination" of the 
previous sample by UCDs and bright GCs with their dis- 
tinct kinematics f !j6.3l) . The kurtosis trends for the blue 
GCs also suggest that the system may be near-isotropic 



inside ~ 50 kpc and more radial at larger distances, which 
might reflect the two-phase assembly discussed earlier, 
with the inner and outer regions forming in fast and slow 
accretion phases, respectively. 

8.7.4. Metallicity Gradients and General Conclusions 

Returning to the original questions raised in <JT] about 
BCG assembly, our overall conclusion is that the outer 
regions of M87 do not show the hot kinematics expected 
if they were formed from infalling cluster material. It is 
possible that here the bona fide cD envelope is still in 
a very early stage of development, with the halo around 
M87 only now beginning to merge with the larger clus- 
ter environment. The multiple transitions at ~ 50 kpc 
(including the appearance of strong substructure) may 
mark the region now experiencing cluster infall. 

Another important feature in this context is the metal- 
licity increase of the M87 BGCs that we have found inside 
~ 20 kpc ( §5.4.1|) . Similar transitions have been found 
in the Milky Way and NGC 1407, at ~ 10 an d - 60 kpc, 
respectively (|Harrisll200H: iForbes et al.l 120111 ). A metal- 
licity transition with radius may be seen as a sign of a 
two-phase assembly s cenario, with the inner component 
havin g formed in situ (jZolotov et al.ll2009t iCoccato et al.1 
[2010l front eFal1l20ll . However, this scenario generally 
describes the overall stellar population (which is mostly 
metal-rich in the inner regions), while here we are con- 
sidering only the metal-poor component. 

Understanding the gradient goes to heart of the GC 
bimodality mystery: when and where were the metal- 
licity subpopulations formed? In each giant galaxy, an 
interruption in the early, intense star formation period 
may have produce d a break in GC metallicity (e.g., 
iBeaslev et all 12002*). with a radial gradient in each sub- 
population reflecting the dissipative nature of this period. 
Alternatively, virtually all of the metal-poor GCs might 
have been acquired through the accretion of smaller 
galaxies, with the inner gradient caused by a small num- 
ber of massive accretors (with more metal-rich stars) 
preferentially settling at smaller radii because of dynam- 
ical friction. However, this interpretation may simply 
shuffle the bimodality puzzle off to the smaller galax- 
ies, and indeed less mas sive galaxies also show evidence 
for G C subpopulations (jPeng et alJl2006t iStrader et al.l 
12006]) . 

There is a critical need for theoretical work to help 
interpret the wealth of kinematical and dynamical data 
now being obtained on M87 and other BCGs. Although 
recent years have seen a few high-re solution simulations 
of BCGs in a cosmolog ical c ontext dRudick et all 120061 
iRuszkowski fc Springell l200l iDolag et all 120101 ) . as far 
as we know, the re has be e n no study of their predicted 
kinematics since Dubinski (1998). Neither has there been 
any study of the long-term dynamical evolution of GCs 
within galaxies assembling in a cosmol ogical context, ex- 
cept for the case of the Milky Way (jPrieto fc Gnedinl 
l200l . 

9. SUMMARY 

We have presented a new photometric and spectro- 
scopic catalog of 737 confirmed GCs around M87, the 
giant elliptical at the center of the Virgo cluster. For 
451 of these GCs, we have new, precise radial veloci- 
ties from observations with the MMT instrument Hec- 
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tospec and the Keck instruments DEIMOS and LRIS. 
Our catalog also includes half-light radii, measured from 
high-resolution HST images, for nearly half of the GCs. 
The majority of these sizes are new measurements from 
archival data. 

The high precision of our new radial velocities, with a 
median uncertainty of 18 km s _1 , has allowed us to iden- 
tify phase-space structures with low velocity dispersions. 
These features are discussed in detail in a companion 
paper (jRomanowsky et aTll2011h and are likely due to a 
recent accretion event. 

The construction of our catalog required a critical as- 
sessment of literature data. A large portion of this effort 
was the classification of individual objects as GCs or fore- 
ground stars. This task is difficult using radial velocities 
alone, partially because of the existence of the substruc- 
tures mentioned above. We used a combination of multi- 
color photometry, velocity, and HST imaging to reliably 
identify nearly all of our objects as stars or GCs; few 
ambiguous objects remain. Conversely, the firm identi- 
fication of many such foreground stars and background 
galaxies has allowed us to accurately delineate the gri 
colors of GCs, enabling better selection for future photo- 
metric and spectroscopic studies. 

We have found that a subset of published radial ve- 
locities, including some of the most extreme values, are 
erroneous, although it is unclear whether or not simi- 
lar issues affect a substantial fraction of literature data. 
We generally (see below) come to kinematical conclusions 
that differ from those of previous work, even using ob- 
servations taken at the same galactocentric radius. The 
causes of these discrepancies are not entirely clear; the er- 
roneous literature measurements may play a role. Other 
factors include the possibility of pervasive substructure 
and the confounding effects of sample selection, such as 
a correlation between GC luminosity and kinematics. 

The half-light radii of confirmed M87 objects in our 
catalog allow us to help clarify the relationship between 
true GCs and UCDs with a more unbiased sample than 
before (since UCDs were not specifically targeted in our 
spectroscopic observations). Defining the certain UCDs 
as luminous objects with sizes > 10 pc, we have dis- 
covered 18 new ones, and found that overall the UCDs 
follow a narrow track in color-magnitude space. We infer 
that the UCDs are a population distinct from a mere ex- 
tension of the bright GCs. Among luminous objects, we 
estimate that few true GCs have sizes > 5 pc, while some 
UCDs range to smaller sizes that overlap with the GCs. 
The U CDs are discussed further in another companion 
paper (|Brodie et al.ll2011ft . 

We have presented new Sersic profile fits to the GC sur- 
face densities of both individual metallicity sub-popula- 
tions (metal-poor and metal-rich) and the full set of M87 
GCs. These fits are more accurate than single power- 
laws and can be readily deprojected to yield the three- 
dimensional number densities required for mass model- 
ing. Blind fitting for the surface density of background 
objects in the Sersic profiles gives similar results to the 
spectroscopic fraction of contaminants, lending credence 
to Sersic fits for GC systems in general. 

The surface density profile of the metal-rich GCs is a 
good match to that of the stellar light of M87 itself over 
an enormous radial range, from ~ 8-100 kpc. Separately, 
we have also estimated the radial profiles of ellipticity e 



for the GC subpopulations. The metal-poor GCs have 
a flat radial profile with e ~ 0.3. The metal-rich GCs 
have e that increases with radius and is generally consis- 
tent with the integrated light of M87. The latter finding, 
along with the radial profile comparison, supports an as- 
sociation between the metal-rich GCs and the bulk of the 
field starlight. 

We have critically examined the gri color distribution 
of GCs as a function of luminosity and galactocentric 
radius, using spectroscopically confirmed GCs to guide 
the analysis. There is evidence for a third old subpop- 
ulation of objects, in addition to the classic blue and 
red groups, with intermediate colors and brighter than 
typical luminosities. Only a fraction of these objects ap- 
pear to be UCDs. This additional subpopulation may 
be partially responsible for the identification of a sig- 
nificant, extended color gradient among the metal-poor 
GCs in previous work. Our photometric analysis yielded 
only marginal evidence for such a gradient outside of 
~ 3'. By contrast, while we have found no evidence 
for a monotonic color gradient among metal-rich GCs, 
the mean color of this subpopulation within 5' (~ 25 
kpc) has much larger radial variations than expected 
from counting statistics. This suggests the presence of 
accreted metal-rich GCs that are still radially unmixed. 

The kinematics of the GC system of M87 have been 
analyzed, while remaining cognizant of the presence of 
substructure in the distribution of radial velocities. Pre- 
vious work found high amounts of rotation (u ro t up to 
~ 400 km s _1 ); our new analysis show little evidence for 
rotation (u ro t < 150 km s _1 at all radii, and typically 
- 20 km s" 1 ). 

We find a much lower value for the outer velocity dis- 
persion than in the literature (with v rms ~ 300 km s _1 ), 
and no indication of a previously claimed transition be- 
tween the outer halo of M87 and the potential of the 
Virgo cluster (which has w rms ~ 800 km s _1 ). In partic- 
ular, the current data suggest that M87, to radii > 100 
kpc, may be dynamically decoupled from Virgo itself. 
We similarly cannot confirm the assertion that M87 has 
an "edge" at ~ 150 kpc. These apparently mistaken con- 
clusions seem to principally be due to the confluence of 
three separate effects: (i) the presence of radial orbits 
and phase-space substructure, (ii) small number statis- 
tics, and (iii) a few catastrophic errors in determining 
GC radial velocities. We observe fewer IGCs than ex- 
pected; the surface density of these objects may be lower 
than reported in the literature, or their kinematics may 
be cold along our line of sight so that they are difficult 
to identify. 

Considering the blue and red GCs separately, we have 
attempted to find evidence for a clear kinematical di- 
chotomy between these two classical subpopulations, but 
have been stymied by the complicating effects of an ad- 
ditional, intermediate-color subpopulation. These sub- 
populations have complex kinematics and are not well- 
described by a single value of rotation or position angle. 
There are indications of unusual kinematics at even finer 
gradations — for example, the reddest GCs inside of 10 
kpc (projected) have a mean velocity that is offset from 
systemic by nearly 200 km s . Conclusions based on 
smaller samples of GCs are necessarily less certain than 
those discussed thus far, and additional data are needed 
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for confirmation of tentative features. 

Using a scale-free analysis we have estimated the mass 
enclosed within ~ 85 kpc of M87 to be (5 ± 2) x 1O 12 M . 
For the GC system as a whole, we infer mildly radial 
orbits with /3 ~ 0.4, with the derived mass only weakly 
dependent on the anisotropy. The mass is lower than 
the most recent X-ray based determination at the same 
radii, although the difference may not be significant given 
the systematic uncertainties in our analysis. Our mass 
results are also much lower than the previous estimate 
using the old GC data set. This suggests that some of 
the discrepancies found between optic al and X-ray based 
mass analysis in various galaxies (e.g.. lShen fc Gebhardti 
I2010D might be driven by problems with the optical data. 

Looking toward the future, we identify the following 
priorities for future improvements. The first is a need 
for complete radial and azimuthal velocity coverage, with 
special attention given to regions with unusual features. 
This would require data in the very center (< 8 kpc) 
where the far-red GCs have an offset v sys ; the region of 
maximal disagreement between new and literature data 
(25-50 kpc); and in the distant outer halo (> 150 kpc), 
hopefully extending at least to 250 kpc where a popula- 
tion of intergalactic PNe has been found, and where the 
transition to the Virgo Cluster may finally be revealed. 
Improvements in the supporting data are equally im- 
portant, including photometry and HST sizes, to enable 
the identification of subpopulations, exclude foreground 
stars, and continue to clarify the relationship between 
GCs and UCDs. 

The extensive chemo-dynamical dataset of GCs around 
M87 provides a unique opportunity for insight into the 
assembly of dark matter and stellar halos around galax- 
ies. The results so far, while bedeviled with complexi- 
ties at all scales, are broadly consistent with the overall 
conclusion that M87 is in active assembly, but still de- 
coupled dynamically from the Virgo cluster. Despite our 
wish list above, in this new era of ultra wide-field, high- 
precision kinematics, data have ceased to be the main 
limiting factor. As near-field cosmology theory starts to 
reach beyond the Local Group, M87 and Virgo would be 
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TABLE 1 

Spectroscopic Observing Log 



Instrument 



Date 



Exposure Time 



Keck/DEIMOS 
MMT/Hectospec 
Keck/DEIMOS 
Keck/LRIS 



2007 Mar 20/21 3-3.5 hr 

2010 Feb 17 2 hr 

2010 Mar 12 0.6 hr 

2010 Apr 8-11 2.5-3.5 hr 



TABLE 2 

Data for 
Keck/DEIMOS 
Spectroscopic 
Candidates 



TABLE 3 
Data for 
Supplementary 
Keck/DEIMOS 
Spectroscopy 



TABLE 4 

Data for 
MMT/Hectospec 
Spectroscopic 
Candidates 



TABLE 5 
Data for Keck/LRIS 
Spectroscopic 
Candidates 



TABLE 6 

Globular Cluster/UCD Radial 
Velocities with Multiple 
Measurements 



Zolotov, A., Willman, B., Brooks, A. M., Governato, F., Brook, 
C. B., Hogg, D. W., Quinn, T., & Stinson, G. 2009, ApJ, 702, 
1058 
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TABLE 7 
New 
Half- 
Light 
Radii 



TABLE 8 

Spectroscopic Data for 
All Foreground Stars 



TABLE 9 

Spectroscopic and 
Photometric Data for 
All M87 Globular 
Clusters and UCDs 



TABLE 10 

Supplementary 
Photometric Data 
for M87 Globular 
Clusters and 
UCDs 



TABLE 11 

Sersic Profile Fits 



Sample 


Population 


N a 


Rs 




m 






(arcmin -2 ) 


(arcmin) 






Phot. 


Full 


2.72 X 10 4 


5.11 X 10" 


-5 


5.50 


Phot. 


Blue 


6.10 X 10 2 


1.62 X 10" 


-2 


3.69 


Phot. 


Red 


5.08 x 10 4 


2.66 x 10" 


-5 


5.33 


Spec. 


Full 


7.36 x 10 4 


4.96 x 10" 


-fi 


6.27 


Spec. 


Blue 


3.22 x 10 2 


4.99 x 10" 


-2 


3.22 


Spec. 


Red 


3.49 x 10 4 


5.51 x 10" 


-5 


5.09 



These normalizations only include GCs with 19 < in < 22.5. 



TABLE 12 
Interloper Surface Density 
Estimates 

Sample Photometric Spectroscopic 
(arcmin -2 ) (arcmin -2 ) 

Full 0.51 ±0.16 0.44 ±0.15 

Blue 0.18 ±0.17 0.25 ±0.11 
Red 0.36 ±0.09 0.37 ±0.19 
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TABLE 13 

Ellipticity and Position Angle Estimates 







l'< R < 3' 


3 '< R < 5' 


5'< R < 8.5' 


8.5' < R < 12' 


Full 




0.24 ± 0.05 


0.25 ± 0.06 


0.31 ± 0.07 


0.32 ± 0.07 


Full 


P. A. 


-39 ± 7 


-21 ± 8 


-41 ± 8 


-29 ±8 


Blue 




0.27 ± 0.07 


0.28 ± 0.07 


0.30 ± 0.07 


0.31 ± 0.08 


Blue 


P.A. 


-38 ±8 


-19 ±8 


-41 ± 8 


-25 ±9 


Red 


e 


0.12 ±0.03 


0.19 ± 0.06 


0.34 ± 0.10 


0.31 ± 0.12 


Red 


P.A. 


-30 ±7 


-32 ±8 


-40 ±9 


-38 ±9 



TABLE 14 

Kinematic al Results 



Subpopulation 


Radius range (R p ) 


N 


00 


Vrot 




"rms 


W p 










(km s- 1 ) 


(km s- 1 ) 


(km s- 1 ) 






All GCs 


1.5' 


-38' 


410 


135° ± 11° 


21 ± 19 


320 ± 11 


320 ± 11 


0.32 ±0.24 


UCDs/bright GCs 


0.7' 


-4T 


73 


185° ± 66° 


17 ±47 


368 ± 30 


368 ± 16 


0.83 ±0.56 


BGCs 


1.5' 


-38' 


242 


100° ± 40° 


32 ± 26 


335 ± 15 


335 ± 15 


0.25 ±0.31 


MGCs 


1.5' 


-32' 


92 


243° ± 29° 


79 ±44 


297 ± 22 


298 ± 22 


0.04 ± 0.50 


RGCs 


1.5' 


-26' 


76 


134° ± 38° 


62 ±42 


295 ± 23 


295 ± 24 


0.97 ± 0.55 


All GCs 


1.5' 


-10' 


250 


186° ± 40° 


38 ± 26 


335 ± 15 


335 ± 15 


0.26 ±0.31 


UCDs/bright GCs 


0.7' 


-10' 


46 


123° ± 72° 


17 ±61 


385 ± 41 


385 ± 41 


0.73 ±0.69 


BGCs 


1.5' 


-10' 


120 


178° ± 49° 


43 ±41 


373 ± 25 


373 ± 24 


-0.08 ± 0.44 


MGCs 


1.5' 


-10' 


64 


225° ± 37° 


71 ± 53 


306 ± 26 


306 ± 28 


0.25 ±0.59 


RGCs 


1.5' 


-10' 


66 


126° ± 54° 


33 ±41 


286 ± 25 


286 ± 25 


1.16 ±0.58 


All GCs 


10' 


-38' 


160 


62° ± 31° 


54 ± 30 


293 ± 16 


294 ± 17 


0.36 ±0.38 


UCDs/bright GCs 


10' 


-4T 


27 


230° ± 27° 


178 ± 104 


338 ± 44 


338 ± 48 


1.46 ± 0.87 


BGCs 


10' 


-38' 


122 


60° ± 25° 


81 ± 35 


289 ± 18 


292 ± 19 


0.57 ±0.44 


MGCs 


10' 


-32' 


28 


276° ± 37° 


97 ± 66 


280 ± 36 


280 ± 39 


-0.53 ±0.86 


RGCs 


IT 


-26' 


10 


121° ±37° 


182 ± 123 


348 ± 66 


348 ± 87 


0.98 ± 1.33 


PNe 


10' 


-37' 


16 


182° ± 46° 


166 ± 124 


387 ± 63 


401 ± 76 


1.91 ± 1.09 



Note. — The UCDs/bright GCs sample incorporates all objects with > 5.25 pc, or io < 20. These objects 
are excluded from all the other subsamplcs, along with the IGCs, H35970, and S923. The blue GCs (BGCs) have 
(g — i)o < 0.86, the red GCs (RGCs) have (g — i)o > 1.01, and the MGCs have intermediate colors. 



